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Antiferromagnetism and superconductivity are both fundamental and common states of matter. 
In many strongly correlated systems, including the high Tc cuprates, the heavy fermion compounds 
and the organic superconductors, they occur next to each other in the phase diagram and influence 
each other's physical properties. The 50(5) theory unifies these two basic states of matter by a 
symmetry principle and describes their rich phenomenology through a single low energy effective 
model. In this paper, we review the framework of the 50(5) theory, and its detailed comparison 
with numerical and experimental results. 
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I. INTRODUCTION 

The phenomenon of superconductivity (SC) is one of the most profound manif estations of quantum mechanics in 
the macroscopic world. The celebrated Bardeen-Cooper-Schrieffer (BCS) theorv ljBardeen et all Il957|) of supercon- 
ductivity provides a basic theoretical framework to understand this remarkable phenomenon in terms of the pairing 
of electrons with opposite spin and momenta to form a collective condensate state. Not only does this theory quan- 
titatively explain the experimental data of conventional superconductors, the basic concepts developed from this 
theory, including the concept of spontaneously broken symmetry, the Nambu-Goldstone modes and the Anderson- 
Higgs mechanism, provide the essential buildin g blocks for the unified theor y of fundame ntal forces. The discovery of 
high temperature superconductivity CHTSCI lfBednorz and Miilleitll986HWu et a^ill987(l in the copper oxide material 
poses a profound challenge to the theoretical understanding of the phenomenon of superconductivity in the extreme 
limit of strong correlations. While the basic idea of electron pairing in the BCS theory carries over to HTSC, other 
aspects, like the weak coupling mean field approximation and the phonon mediated pairing mechanism, may not ap- 
ply without modifications. Therefore, the HTSC systems provide an exciting opportunity to develop new theoretical 
frameworks and concepts for strongly correlated electronic systems. 

Since the discovery of HTSC, a tremendous amount of experimental data has been accumulated on this material. 
In this theoretical review it is not possible to give a detailed review of all the experimental findings. Instead, we 
refer the readers to a number of excellent review articles on the subjectljCamp uzano et al\. |2002|: [Damascelli et all 
gnn.l't'Tm ada et nil . Il 9981: iKastner et oIl . ll998HMa,n]eHl 998tl(>^nstein and Mi]li4l2nnflHTimusk and Stattlll flflflHYeTl 
2002) . Below, we summarize the phase diagram of the HTSC cuprates and discuss some of the basic and (more or 
less) universal properties in each phase. 

To date, a number of different HTSC materials have been discovered. The most studied of these include the 
hole doped La2-xSr^CuOi+s (LSCO), YBa2CuzOfi+5 (YBCO), Bi2Sr2CaCu20s+s (BSCO), and Tl2Ba2Cu06+s 
(TBCO) materials and the electron doped Nd2-xCexCu04 (NCCO) material. All these materials have two di- 
mensional (2D) Cu02 planes and display an antiferromagnetic (AF) insulating phase at half-filling. The magnetic 
properties of this insulating phase are well approximated by the AF Heisenberg model with spin S* = 1/2 and an 
AF exchange constant J ^ lOOme^. The Neel temperature for the three dimensional AF ordering is approximately 
given by Tn ^ 300 — 500K. The HTSC material can be doped either by holes or by electrons. In the doping range 
of 5% ^ a; < 15%, there is an SC phase, which has a dome-like shape in the temperature versus doping plane. The 
maximal SC transition temperature, Tc, is of the order lOOK . The three doping regimes are divided by the maximum 
of the dome and are called the underdoped, optimally doped, and overdoped regimes, respectively. The generic phase 
diagram of HTSC is shown in Fig. ^ 

One of the main questions concerning the HTSC phase diagram is the transition region between the AF and SC 
phases. Partly because of the complicated material chemistry in this regime, there is no universal agreement among 
different experiments. Different experim ents indicate several different po ssibilities, including phase separation with 
an inho mogeneous density distribution llHowald et all f200lt iLang et all l2002|) . uniform mixed phase between AF 
and SC l|Brewer et aZ.LIl988HMiller et allhOO^ and periodically ordered spin and charge distributions in the form of 



stripes or checkerboards (Tranouad a et oil 1X993) . 

The phase diagram of the HTSC cuprates also contains a regime with anomalous behavior conventionally called 
the pseudogap phase. This region of the phase diagram is indicated by the dashed line in Fig. ^ In conventional 
superconductors, a pairing gap opens up at Tc- In a large class of HTSC cuprates, however, a gap, which can be 
observed in a variety of spectroscopic experiments, starts to open up at a temperature T* , much higher than Tc- Many 
experiments indicate that the pseudogap "phase" is not a true thermodynamical phase but rather a precursor toward 
a crossover b ehavior. The ph enomenology of the pseudogap behavior is extensively reviewed in ljTallon and LoramL 
EOOI; Timusk and Stattl . Il999() . 

The SC phase of the HTSC has a number of striking properties not shared by conventional superconductors. 
First of a ll, phase sensitive experiments ind icate that the SC phase of most cuprates has d-wave like pairing 
symmetry ijHarlineenl 119951: iTsuei and Kirtlevl [JOOO). This is also sup ported by the photoemi s sion experiments , 
which show the existence of the nodal points in the quasiparticle gap ijCampuzano et oil . 120021 : iDamascelli et al\ . 
l2003i) . Neutron scattering experiments find a new type of collective mode, carrying spin one, lattice momentum close 
to (tt, tt) and a resolution limited sharp resonance energy around 20 ~ AOmeV. Most remarkably, this resonance mode 
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appears only below Tc in the optimally doped c uprates. It has been found in a number of materials, inc l uding the 
YBCOBSCO and the TBCO classe s of materials llDai et oi.Lll99fillT993lFong et all ll999H200d 119951 ll99fil:lHe et all 
l2nn2ll2nmtfMook e.t oiJ . ri9Mll993trR.ossat,-Mignod p.t oil 1199111 . Another property uniquely different from the con- 
ventional superconductors is the vortex state. Most HTSC are type II superconductors in which the magnetic field 
can penetrate into the SC state in the form of a vortex lattice, with the SC order being destroyed at the center of 
the vortex core. In conventional superconductors, the vortex core is filled by normal metallic electrons. However, a 
number of different experimental probes, including neutron scattering, muon spin resonance (/xsR), and nuclear mag- 



normal metallic! Fuiita et al|,| 


2003 


ikakuvanagi et a/.|.l2003; Kang et a/1 12003: Katano et aL. 200G:|Khaykovich et al\. 


120021 iLake et aL. 


,2001^ 


|200^ 


Lev: 


,r2002£ .Miller et al. 


.l200aiMitrovic et a/J.l200L.,2003(). This phenomenon has been 



observed in almost all HTSC materials, including LSCO, YBCO, TBCO and NSCO; thus, it appears to be a universal 
property of the HTSC cuprates. 

The HTSC materials also have highly unusual transport properties. While conventional metals have a dependence 
of resistivity, in accordance with the predictions of the Fermi liquid theory, the HTSC materials display a linear T 
dependence of the resistivity near optimal doping. This linear T dependence extends over a wide temperature window 
and seems to be universal among most of the cuprates. When the underdoped or sometimes optimally doped SC 
state is destroyed by applying a high magneti c field, the resulti ng "n ormal state" is not a conventional conducting 
state (|Ando et all 119951 119961: Boebinger et all Il996t iHill et all l200l|) but exhibits insulating like behavior, at least 
along the c axis, i.e. the axis perpendicular to the Cu02 planes. This phenomenon may be related to the insulating 
AF vortices mentioned in the previous paragraph. 

The HTSC materials attracted great attention because of the high SC transition temperature. However, many of the 
striking properties discussed above are also shared by other materials, which have a similar phase diagram but typically 
with much reduced temperature and energy scales. The 2D organic superconductor k— {BEDT —TTF)2X (X=anion) 
display a similar phase diagram in the temperature versus pressure plane, where a direct f irst order transition between 
the AF and SC phases can be tuned by pressure or magnetic field l)Lefebvre et all 1200(1 ISingleton and Miel kc. 2002 
iTaniguchi et alV 120031) . In this system, the AF transition temperature is approximately Tiy ^ 30j^, while th e 
SC transition temper ature is Tr, ~ I'iK. In heavy fermion compounds CeCu'^iSi^-^Ge■,■^■:l ^\ta.6ka, et ai\ . l200l|) . 
CePd2Si2 and Ce Jn-; ijMathur et"aZI . Il998|) . the SC phase also appears near the boundary to the AF phase. In all 
these systems, even though the underlying solid state chemistries are rather different, the resulting phase diagrams are 
strikingly similar and robust. This similarity suggests that the overall feature of all these phase diagrams is controlled 
by a single energy scale. Different classes of materials differ only by this overall energy scale. Another interesting 
example of competing AF and SC can be found in quasi-one-dimensional Bechgaard salts. The most well studied 
material from this famil y, (TMTSF)2PFf,, i s an AF insulator at ambient pressure and becomes a triplet SC above a 
certain critical pressure ijJerome et al. l ll98nHLee ^7^11. l2003[ Il997t IVuletic et ai.Ll2002^ . 

The discovery of HTSC has greatly stimulated the theoretical understanding of superconductivity in strongly 
correlated systems. Since the theoret ical literature is extensive, ^the readers are refer r ed to a number of excel- 
lent reviews and re present ative artic les fAbrikosov, 2000; Anderson, '1997'; Anderson et al., 2003; Balcnts et a/., 199g 



Carlson et aH 12002: .Chakravartv et al 2001: Chubukoy et al 2002: Dagotto. 1994; Franz et al 2002b: Fu et al. 

2004 [ [Inui et aJl 119881: lloffe and Milli^ 120021: iLaughlinL l2002t iNorman and Pen ni l2003t ISachdevl l2002at IScalapincl 

lis 



1995^ 'Sc hrieffer a/J . 119891: ISenthil and Fisheil 120011: IShen et all 120021: IVarmai 119991 IWen and Leei Il99(it IZaancE . 



il999b). The present review article focuses on a part icular theory, which unifies the AF and SC phases of the HTSC 
cuprates based on an approximate 50(5) svmmetrv llZhand Il997j) . The SO{5) theory draws its inspiration from the 
successful application of symmetry principles in theoretical physics. All fundamental laws of Nature are statements 
about symmetry. Conservation of energy, momentum and charge are direct consequences of global symmetries. The 
form of fundamental interactions are dictated by local gauge symmetries. Symmetry unifies apparently different 
physical phenomena into a common framework. For example, electricity and magnetism were discovered indepen- 
dently and viewed as completely different phenomena before the 19th century. Maxwell's theory and the underlying 
relativistic symmetry between space and time unified the electric field, E, and the magnetic field, B, into a common 
electromagnetic field tensor, i^^^. This unification shows that electricity and magnetism share a common microscopic 
origin and can be transformed into each other by going to different inertial frames. As discussed in the introduction, 
the two robust and universal ordered phases of the HTSC are the AF and SC phases. The central question of HTSC 
concerns the transition from one phase to the other as the doping level is varied. 

The SO{5) theory unifies the three dimensional AF order parameter {N^, Ny, N^) and the two dimensional SC 
order parameter (i?eA, JmA) into a single, five dimensional order parameter called the superspin, in a way similar to 
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the unification of electricity and magnetism in Maxwell's theory: 
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This unification relies on the postulate that a common microscopic interaction is responsible for both AF and SC 
in the HTSC cuprates and related materials. A well-defined 5*0(5) transformation rotates one form of the order 
into another. Within this framework, the mysterious transition from the AF to SC phase as a function of doping is 
explained in terms of a rotation in the five dimensional order parameters space. Symmetry principles are not only 
fundamental and beautiful, but they are also practically useful in extracting information from a strongly interacting 
system, which can be tested quantitatively. As seen in the examples applying the isospin SU(2) and the SU{3) 
symmetries to the strong interaction, some quantitative predictions can be made and tested even when the symmetry 
is broken. The approximate SO{5) symmetry between the AF and SC phases has many direct consequences, which 
can be tested both numerically and experimentally. We shall discuss a number of these tests in this review article. 

Historically, the 5*0(5) theory concentrated on the competition between AF and SC orders in the high Tc cuprates. 
The idea of some order competing with supercon ductivity is common in several theories. The staggered flux or the 
d-density wave phase has been suggested in Refs. (Affleck and Marston.'l988VChakravartv et aL'.'200lVWen an d Leel 
Kme,) . the spin-Peierls order has been discussed in Refs. (Park and Sach dev, 2001; Vo ita and Sachde y, 1990,), and 
spin and charge density wave orders have been considered in Refs. (jKivelson et a?J . f20blllZaanenlll999aHzhang all 
|2002() . The 5*0(5) theory extends simple consideration of the competition between AF and SC in the cuprates by 
unifying the two order parameters using a larger symmetry and examining consequences of such symmetry. 

The microscopic interactions in the HTSC materials are highly complex, and the resulting phenomenology is 
extremely rich. The 50(5) theory is motivated by a confluence of the phenomenological top-down approach with the 
microscopic bottom-up approach, as discussed below. 

The top-down approach: Upon first glance at the phase diagram of the HTSC cuprates, one is immediately impressed 
by its striking simplicity; there are only three universal phases in the phase diagram of all HTSC cuprates: the AF, 
the SC and the metallic phases, all with homogeneous charge distributions. However, closer inspection reveals a 
bewildering complexity of other possible phases, which may not be universally present in all HTSC cuprates, and 
which may have inhomogencous charge distributions. Because of this complexity, formulating a universal theory of 
HTSC is a formidable challenge. The strategy of the 5*0(5) theory can be best explained with an analogy: we see 
a colorful world around us, but the entire rainbow of colors is composed of only three primary colors. In the 5*0(5) 
theory, the superspin plays the role of the primary colors. A central macroscopic hypothesis of the 50(5) theory is that 
the ground state and the dynamics of collective excitations in various phases of the HTSC cuprates can be described in 
terms of the spatial and temporal variations of the superspin. This is a highly constraining and experimentally testable 
hypothesis, since it excludes many possible phases. It does include a homogeneous state in which AF and SC coexist 
microscopically. It includes states with spin and charge density wave orders, such as stripe phases, checkerboards 
and AF vortex cores, which can be obtained from spatial modulations of the superspin. It also includes quantum 
disordered ground states and Cooper pair density wave, which can be obtained from the temporal modulation of the 
superspin. The metallic Fermi liquid state on the overdoped side of the HTSC phase diagram seems to share the same 
symmetry as the high temperature phase of the underdoped cuprates. Therefore, they can also be identified with the 
disordered state of the superspin, although extra care must be given to treat the gapless fermionic excitations in that 
case. If this hypothesis is experimentally proven to be correct, a great simplicity emerges from the complexity: a full 
dynamical theory of the superspin field can be the universal theory of the HTSC cuprates. Part of this review article 
is devoted to describing and classifying phases which can be obtained from this theory. This top-down approach 
focuses on the low energy collective degrees of freedom and takes as the starting point a theory expressed exclusively 
in terms of these collective degrees of freedom. This is to be contrasted with the conventional approach based on weak 
coupling Fermi liquid theory, of which the BCS theory is a highly successful example. For an extensive discussion 
on the relative merits of both approaches for the HTSC problem, we refer the readers to an excellent, recent review 
article in Ref. l|Carlson et a/.l 120021. 

The 50(5) theory is philosophically inspired by the Landau-Ginzburg (LG) theory. The LG theory is a highly 
successful phenomenological theory, in which one first makes observations of the phase diagram, introduces one order 
parameter for each broken symmetry phase and constructs a free energy functiona l by expanding in terms of different 
order parameters (a review of earlier work based on this approach is given in Ref. l|Vonsovskv et alUldS^ ). However, 
given the complexity of interactions and phases in the cuprates, introducing one order parameter for each phase with 
unconstrained parameters would greatly compromise the predictive power of theory. The 5*0(5) theory extends the LG 
theory in several important directions. First, it postulates an approximately symmetric interaction potential between 
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the AF and the SC phases in the underdoped regime of the cuprates, thereby greatly constraining theoretical model 
building. Second, it includes a full set of dynamic variables canonically conjugate to the superspin order parameters, 
including the total spin, the total charge and the so called tt operators. Therefore, unlike the classical LG theory, 
which only contains the classical order parameter fields without their dynamically conjugate variables, the 15*0(5) 
theory is capable of describing quantum disordered phases and the quantum phase transitions between these phases. 
Because the quantum disordered phases are described by the degrees of freedom canonically conjugate to the classical 
order parameters, a definite relationship, the so-called SO{5) orthogonality relation, exists between them, which can 
give highly constrained theoretical predictions. Therefore, in this sense, the S0{5) theory makes great use of the 
LG theory but also goes far beyond in making more constrained and more powerful predictions which are subject to 
experimental falsifications. 

The bottom-up approach: Soon after the discovery of the HTSC cuprates, Anderson l) Andersonl Il987j) introduced 
the repulsive Hubbard model to describe the e l ectron ic degrees of freedom in the Cu02 plane. Its low energy limit, 
the t — J model, is defined bv l)Zhang and Ric^ - llOSSj) 



H = -t 



E 

{x.x') 



c^„{x)c„{x ) + J ^ S{x) ■ S{x'), 

{x.,x') 



(2) 



where the t term describes the hopping of an electron with spin a from a site x' to its nearest neighbor x, with double 
occupancy removed, and the J term describes the nearest neighbor spin exchange interaction. The main merit of 
these models does not lie in the microscopic accuracy and realism but rather in the conceptual simplicity. However, 
despite their simplicity, these models are still very hard to solve, and their phase diagram cannot be compared directly 
with experiments. 

The t — J model certainly contains AF at half-filli ng. While it is s till not fully settled whether it has d-wave SC 
ground state with a high tr ansition temperatu re l,Prvadko et a^J . |2003j) . it is reasonably convincing that it has strong 
d-wave pairing fluctuations l|Sorella et a/.l . [2002f) . Therefore, it is plausible that a small modification could give a robust 
SC ground state. The basic microscopic hypothesis of the SO{5) theory is that AF and SC states arise from the same 
interaction with a common energy scale of J. This common energy scale justifies the treatment of AF and SC on 
equal footing and is also the origin of an approximate 50(5) symmetry between these two phases. By postulating an 
approximate symmetry between the AF and SC phases, and by systematically testing this hypothesis experimentally 
and numerically, the question of the microscopic mechanism of HTSC can also be resolved. In this context, early 
numerical diagonalizations sho wed that the low-l ying states of the t — J model fit into irreducible representations 
of the 50(5) symmetry group ijEder et all Il99^ . If the SO{5) symmetry is valid, then HTSC shares a common 
microscopic origin with the AF, which is a well understood phenomenon. 

The basic idea is to solve these models by two steps. The first step is a renormalization group transformation, 
which maps these microscopic models to an effective superspin model on a plaquette, typically of the size of 2a x 2a 
or larger. This step determines the form and the parameters of the effective models. The next step is to solve the 
effective model either through accurate numerical calculations or tractable analytical calculations. 

There is a systematic method to ca rry out the first step. Using the contractor- renormalizati on-group 
rCORE UMorningst ar and WeinsteinL Il996|) approach, Altman and Auerbach l|Altman and Auerbachl 1200^ derived 
the projected SO{5) model from the Hubbard and the t — J model. Within the approximations studied to date, a 
simple and consistent picture emerges. There are only five low energy states on a coarse grained lattice site, namely 
a spin singlet state and a spin triplet state at half-filling and a d-wave hole pair state with two holes. These states 
correspond exactly to the local and dynamical superspin degrees of freedom hypothesized in the top-down approach. 
The resulting effective S0{5) superspin model, valid near the underdoped regime, only contains bosonic degrees of 
freedom. This model can be studied by quantum Monte-Carlo simulations up to very large sizes, and the accurate 
determination of the phase diagram is possible (in contrast to the Hubbard and t — J models) because of the absence 
of the minus sign problems. Once the global phase diagrams are determined, fermionic excitations in each phase can 
also be studied by approximate analytic methods. Within this approach, the effective 50(5) superspin model derived 
from the microscopic physics can give a realistic descriptio n of the phenomenology and phase diagram of the HTSC 
cuprates and account for many of their physical properties l)Dorneich et a^J . l2002albtl . This agreement can be further 
tested, refined and possibly falsified. This approach can be best summarized by the following block diagram: 



Hubbard and t-J 
model 
of the electron 



coarse graining 
transformation 



Quantum SO (5) 

model 
of the superspin 



analytical and numerical 
calculations 



Phase diagram 
Collective modes 
other experiments.. 



While the practical execution of the first step already introduces errors and uncertainties, we need to remember 
that the Hubbard and the t — J models are effective models themselves, and they contain errors and uncertainties 
compared with the real materials. The error involved in our coarse grain process is not inherently larger than the 
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uncertainties involved in deriving the Hubbard and the t — J models from more realistic models. Therefore, as long 
as we study a reasonable range of the parameters in the second step and compare them directly with experiments, 
we could determine these parameters. 

This review is intended as a self-contained introduction to a particular theory of the HTSC cuprates and related 
materials and is organized as follows. Section ^] describes three toy models which introduce some important concepts 
used in the rest of the article. Section lTTll introduces the concept of the 5*0(5) superspin and its symmetry transforma- 
tion, as well as effective dynamical models of the superspin. The global phase diagram of the SO{h) model is discussed 
and solved numerically in section llVl SectionlVl introduces exact 50(5) symmetric microscopic models, the numerical 
tests of the 50(5) symmetry in the t — J and Hubbard models, and the Altman-Auerbach procedure of deriving the 
50(5) model from microscopic models of the HTSC cuprates. Section rVTI discusses the tt resonance model and the 
microscopic mechanism of HTSC. Finally, in section lVlII we discuss experimental predictions of the SO{5) theory and 
comparisons with the tests performed so far. The readers are assumed to have general knowledge of quantum many 
body physics an d are referred to severa l excellent t extbooks for pedagogical introductions to the basic concepts and 
theoretical tools llAbrikosov et all Il99.ll: lAndersonL Il997l: lAiierbachL Il994l: iPines and NozieresL hmd ISachdevl l200nt 
ISchrieffeii Il964t iTinkhamL Il995l) . 



II. THE SPIN-FLOP AND THE MOTT INSULATOR TO SUPERCONDUCTOR TRANSITION 

Before presenting the full 50(5) theory, let us first discuss a much simpler class of toy models, namely the anisotropic 
Heisenberg model in a magnetic field, the hard-core lattice boson model and the negative U Hubbard model. The low 
energy limits of this class of models are all equivalent to each other and can be described by a universal quantum field 
theory, the 50(3) quantum non-linear sigma model. Although these models are simple to solve, they exhibit some of 
the key properties of the HTSC cuprates, including strong correlation, competition of different orders, low superfluid 
density near the insulating phase, maximum of T^, and the pseudogap behavior. 

The spin 1/2 anisotropic AF Heisenberg model on a square lattice is described by the following Hamiltonian: 

H= J2 S'{x)Vix,x')S'ix') + J Y,(.S''ix)S''{x) + Syix)Syix))~Bj2S"ix). (3) 

{x,'x') {x,x') X 

Here, 5" = ^r" is the Heisenberg spin operators and t" is the Pauli matrix. J describes the nearest-neighbor 
exchange of the xy components of the spin, while V{x, x') describes the z component of the spin interaction. We shall 
begin by considering only the nearest neighbor (denoted by {x,x')) spin interaction V. B is an external magnetic 
field. At the point of _B = and J = V , this model enjoys an 50(3) symmetry generated by the total spin operators: 

5"-^5"(a;) , =*e„/3^5^ , [iJ,5"]=0.. (4) 

X 

The order parameter in this problem is the Neel operator, which transforms according to the vector representation of 
the 50(3) group 

N'- ^Y^i^^yS^ix) , [S'',NP]=ie^p^N\ (5) 

X 

Here {~Y = 1 if a; is on an even site and (— )^ — —1 if x is on an odd site. The symmetry generators and the order 
parameters are canonically conjugate degrees of freedom, and the second part of Eq. ((SJ is similar to the Heisenberg 
commutation relation [x,p\ — ih between the canonically conjugate position and momentum. Just like p can be 
expressed as jd/dx, one can express 

where the second part of the equation, called the 50(3) orthogonality relation, follows directly from the first. Both the 
symmetry algebra, the canonical conjugation and the orthogonality constraint are fundamental concepts important 
to the understanding of the dynamics and the phase diagram of the model. 

Let us first consider the classical, mean field approximation to the ground state of the anisotropic Heisenberg model 
defined in Eq.Q. For V > J , the spins like to align antiferromagnetically along the z direction. In the Ising phase, 
5^(a;) ~ (—1)^5, the ground state energy per site is given by eising{B) — —-^5^, where z is the coordination number, 
which is 4 for the square lattice. Note that the energy is independent of the B field in the Ising phase. For larger values 
of S, the spins "flop" into the XY plane, and tilt uniformly toward the Z axis. (See Fig. [2^). Such a spin-flop state is 
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given by S^{x) = ScosO and S^{x) = (—1)^ SsinO. The minimal energy configuration is given by cosO = B/zS{V + J), 
and the energy per site for this spin-flop state is exviB) = —zJS^/2 — B^/2z{V + J). Comparing the energies of 
both states, we obtain the critical value of B where the spin-flop transition occurs: Bd — zS^/V^ — .P. On the other 
hand, we require |cos0| < 1, which implies a critical field Bc2 = zS{V + J) at which |cos0| = 1, and the staggered 
order parameter vanishes. Combining these phase transition lines, we obtain the "class i?" transition in the ground 
state phase diagram depicted in the B-J/V plane (see Fig. Oli). Here and later in the article, the "class B" transition 
refers to the transition induced by the chemical potential or the magnetic field. At the SO{2>) symmetric point, V — J 
and Bel = 0. For F < J, the ground state has XY order even at i? = 0, and there is no spin-flop transition as a 
function of the magnetic field B. The Ising to XY transition can also be tuned by varying J /V aX B — 0, and the 
phase transition occurs at the special S'0(3) symmetric Heisenberg point. This type of transition is also depicted in 
Fig. 12^ and will be called "class A" transition in this paper. 

The spin 1/2 Heisenberg model can be mapped to a hard-core boson model, defined by the following Hamiltonian: 

H=Y^ n{x)V{x,x')n{x')~'^J ^ {h\x)h{x') + h.c.) - ^x'Y,n{x). (7) 

{x,x') {x.x') X 

Here h{x) and h^{x) arc the hard-core boson annihilation and creation operators and n{x) = h^{x)b{x) is the boson 
density operator. In this context, V , J and ^ describe the interaction, hopping and the chemical potential energies, 
respectively. There are two states per site; |1) and |0) denote the filled and empty boson states, respectively. They 
can be identified with the spin up | t) and the spin down | [) states of the Heisenberg model. The operators in the 
two theories can be identified as follows: 

bix)"^ ^ {S'-'{x)+iSy{x)) b{x) ^ {~Y{S''{x)-iSy{x)) n{x) ^ S'{x) + ^. (8) 

We see that these two models are identical to each other when n — B + zV/2. From this mapping, we see that the 
spin-flop phase diagram has another interpretation: the Ising phase is equivalent to the Mott insulating phase of 
bosons with a charge-density-wave (CDW) order in the ground state. The XY phase is equivalent to the superfluid 
phase of the bosons. The two paramagnetic states correspond to the full and empty states of the bosons. While 
Heisenberg spins are intuitively associated with the 50(3) spin rotational symmetry, lattice boson models generically 
have only a U{1) symmetry generated by the total number operator N = '^,j.n(x), which transforms the boson 
operators by a phase factor: b^{x) e*"6^(a;) and b{x) e~*"6(x). From this point of view, it is rather interesting 
and non-trivial that the boson model can also have an additional SO (3) symmetry at the special point J — V because 
of its equivalence to the Heisenberg model. 

Having discussed the Heisenberg spin model and the lattice boson models, let us now consider a fermion model, 
namely the negative U Hubbard model, defined by the Hamiltonian 

H^-tJ2 (4(a;)c.(x') + /i.c.) + C/^(nT(a;)-i)(ni(a;)-i)-/i5I"-W' (9) 

{x,x') X X 

where Ca-{x) is the fermion operator and n^{x) = c^^{x)ca{x) is the electron density operator at site x with spin a. t, 
U and fi are the hopping, interaction and the chemical potential parameters respectively. The Hubbard model has a 
pseudospin SU (2) symmetry generated by the operators 

X a 

where rj^ = -rf rk ir\y and a = x,y, z, as before. Yang and Zhang ljYand . Il989t lYang and Zhanei llQQClt IZhand Il99(]l) 

pointed out that these operators commute with the Hubbard Hamiltonian when ii = { i.e. [H, rj"] = 0) ; therefore, 
they form the symmetry generators of the model. Combined with the standard SU{2) spin rotational symmetry, the 
Hubbard model enjoys a 50(4) = SU{2) (g) SU{2)/Z2 symmetry. This symmetry has important consequences in the 
phase diagram and the collective modes in the system. In particular, it implies that the SC and CDW orders are 
degenerate at half-filling. The SC and the CDW order parameters are defined by 

A-=^c.TC4 , A+-(A-)t , A^ = i^(-l)'n,, , [rj'^ , A^] ^ te^p^A\ (11) 

i ia 

where A^ = A^ ±iA^. The last equation above shows that the rj operators perform the rotation between the SC and 
CDW order parameters. Thus, ry" is the pseudospin generator and A" is the pseudospin order parameter. Just like the 
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total spin and the Neel order parameter in the AF Heisenberg model, they are canonically conjugate variables. Since 
[H, 77"] = at /i = 0, this exact pseudospin symmetry implies the degeneracy of SC and CDW orders at half-filling. 

The phase diagram of the U < Q Hubbard model corresponds to a ID slice of the 2D phase diagram, as depicted 
in Fig|21i. The exact pseudospin symmetry implies that the "class _B" transition line for the U < Q Hubbard model 
exactly touches the tip of the Mott lobe, as shown by the B' line in Fig|3^. At /i = 0, SC and CDW are exactly 
degenerate, and they can be freely rotated into each other. For /i 7^ 0, the system is immediately rotated into the 
SC state. One can add additional interactions in the Hubbard model, such as a nearest neighbor repulsion, which 
breaks the SU{2) pseudospin rotation symmetry even at ^ = 0. In this case, the pseudospin anisotropy either picks 
the CDW Mott insulating phase or the SC phase at half-filling. By adjusting the nearest neighbor interaction, one 
can move the height of the "class i?" transition line. 

We have seen that the hard-core boson model is equivalent to the Heisenberg model because of the mapping (jS)). 
The U < Q model, on the other hand, is only equivalent to the Heisenberg model in the low energy limit. In fact, it 
is equivalent to a [/ > Hubbard at half-filling in the presence of a Zeeman magnetic field. The ground state of the 
half-filled Hubbard model is always AF; therefore, its low energy limit is the same as that of the Heisenberg model in 
a magnetic field. All three models are constructed from very different microscopic origins. However, they all share the 
same phase diagram, symmetry group and low energy dynamics. In fact, these universal features can all be captured 
by a single effective quantum field theory model, namely the 50(3) quantum non-linear a model. This model can 
be derived as an effective model from the microscopic models introduced earlier or it can be constructed purely from 
symmetry principles and the associated operator algebra as defined in Eq. Q and (|SJ). The fact that both derivations 
yield the same model is hardly surprising, since the universal features of all these models are direct consequences of 
the symmetry. 

The 50(3) non- linear cr model is defined by the following Lagrangian density for a unit vector field Uq with n\ — \: 

C = ^ujI^P ~ ^{diUaf - V {n) , Uaf} ^ UaidtUfi - iBp^n^) - [a ^ 13), (12) 

where the Zeeman magnetic field is given by Ba — 5 £0/37-8/37- Without loss of generality, we pick the magnetic field 
B to be along the z direction, x p are the susceptibility and stiffness parameters and V{n) is the anisotropy 
potential, which can be taken as V{n) — — fn^. Exact 50(3) symmetry is obtained when g ^ B = Q. g > 
corresponds to easy axis anisotropy or J/V < 1 in the Heisenberg model, g < corresponds to easy plane anisotropy 
or J/V > 1 in the phase diagram of FigO In the case of 5 > 0, there is a phase transition as a function of B. To see 
this, let us expand the first term in (|12|l in the presence of the B field. The time independent part contributes to an 
effective potential K// = V{n) — -§-(n^ -|- n^), from which we see that there is a phase transition at Bd — \fgjx- 
For B < Bel, the system is in the Ising phase, while for B > Bd the system is in the XY phase. Therefore, tuning B 
for a fixed g > traces out the "class B" transition line, as depicted in Fig|3^. On the other hand, fixing B = Q and 
varying g traces out the "class A" transition line, as depicted in Fig|3^. Therefore, we see that the SO{3) non-linear 
a model has a similar phase diagram as the microscopic models discussed earlier. For a more detailed discussion 
of phase transitions in S O (3) non-linear a models we refer the readers to an excellent review paper by Auerbach et 
al l|Auerbach et all l200(ll) . 

In D = 2, both the XY and the Ising phase can have a finite temperature phase transition into the disordered state. 
However, because of the Mermin- Wagner theorem, a finite temperature phase transition is forbidden at the point 
B = g = 0, where the system has an enhanced 50(3) symmetry. The finite temperature phase diagram is shown 
in Fig. ^p. Approaching from the SC side, the Kosterlitz-Thouless transition temperature Txy is driven to zero at 
the Mott to superfluid transition point J/V — 1. In the 2D XY mod el, the superfluid dens i ty an d the transition 
temperature Txy are related to each other by a universal relationship ijNelson and Kosterlitd Il977|) : therefore, the 
vanishing of Txy also i mplies the vanishing o f the superfiuid densit y as one approaches the M ott to superfiuid 
transition. Scalettar et al l)Scalettar et aZlll989f) . Moreo and Scalapino jMoreo and Scalapind . llOQlf) have performed 
extensive quantum Monte Carlo simulation in the negative U Hubbard model and have indeed concluded that the 
superfluid density vanishes at the symmetric point. The 50(3) symmetric point leads to a large regime below the mean 
field transition temperature where fluctuations dominate. The single pa rticle spectral function of the 2D attractive 
Hubbard model has been studied extensively by Allen et al (Allen et ai\ . Il999|) near half-filling. They identified the 
pseudogap behavior in the single particle density of states within this fluctuation regime. Based on this study, they 
argued that the pseudogap behavior is not only a consequence of the SC phase fluctuations toon iach and Inuilll99Ct 
lEmerv and Kivclso n. 1995; Ucmura, 200'^) but also a consequence of the full 50(3) symmetric fluctuations, which 
also include the fluctuations between the SC and the CDW phases. Fig. ^ shows the generic flnite temperature 
phase diagram of these 50(3) models. In this case, the Ising and the XY transition temperatures meet at a single 
bi-critical point T{,c, which has the enhanced 50(3) symmetry. At the "class A" transition point g = B = Q, the 
quantum dynamics is fully 50(3) symmetric. On the other hand, at the "class i?" transition point T = The, only the 
static potential is 50(3) symmetric. We shall return to a detailed discussion of this distinction in section llll.CI 
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The pseudospin SU(2) symmetry of the negative U Hubbard model has another important consequence. Away from 
half-fiUing, the rj operators no longer commute with the Hamiltonian, but they are eigen-operators of the Hamiltonian, 
in the sense that 

[H,r]^]^T2^^V^. (13) 

Thus, the ry operators create well defined collective modes in the system. Since they carry charge ±2, they usually 
do not couple to any physical probes. However, in a SC state, th e SC order parameter mixes the ri oper ators with 
the CDW operator A^, via Eq. l(TT|l . From this reasoning, Zhang ljPemler et ol 1 Il99fit IZhand. Il99nt ll99lD predicted 
a pseudo-Goldstone mode in the density response function at wave vector (tt, tt) and energy — 2/i, which appears only 
below the SC transition temperature Tc- This prediction anticipated the neutron resonance mode later discovered in 
the HTSC cuprates; a detailed discussion shall be given in section IVTl 

From the toy models discussed in this section, we learned a few very important concepts. Competition between 
different orders can sometimes lead to enhanced symmetries at the multi-critical point. Universal properties of very 
different microscopic models can be described by a single quantum field theory constructed from the canonically 
conjugate symmetry generators and order parameters. The enhanced symmetry naturally leads to a small superfluid 
density near the Mott transition. The pseudogap behavior in the single particle spectrum can be attributed to the 
enhanced symmetry near half-filling and new types of collective Goldstone modes can be predicted from the symmetry 
argument. All these behaviors are reminiscent of the experimental observations in the HTSC cuprates. The simplicity 
of these models on the one hand and the richness of the phenomenology on the other inspired the 50(5) theory, which 
we shall discuss in the following sections. 



III. THE S0(5) GROUP AND EFFECTIVE THEORIES 

A. Order parameters and S0(5) group properties 

The 50(3) models discussed in the previous section give a nice description of the quantum phase transition from 
the Mott insulating phase with CDW order to the SC phase. However, these simple models do not have enough 
complexity to describe the AF insulator at half-filling and the d-wave SC order away from half-filling. Therefore, a 
natural step is to generalize these models so that the Mott insulating phase with the scalar CDW order parameter 
is replaced by a Mott insulating phase with the vector AF order parameter. The pseudospin 50(3) symmetry group 
considered previously arises from the combination of one real scalar component of the CDW order parameter with 
one complex or two real components of the SC order parameter. After replacing the scalar CDW order parameter 
by the three components of the AF order parameter and combining them with the two components of the SC order 
parameters, we are naturally led to consider a five component order parameter vector and the SO(5) symmetry group 
which transforms it. 

In section^ we introduced the crucial concept of order parameter and symmetry generator. Both of these concepts 
can be defined locally. In the case of the Heisenberg AF, at least two sites, for instance, 5*1 and ^2, are needed to 
define the total spin S — Si + S2 and the Neel vector N = Si — 82- Similarly, it is simplest to define the concept 
of the 50(5) symmetry generator and order parameter on two sites with fermion operators and da-, respectively, 
where tr = 1, 2 is the usual spinor index. The AF order parameter operator can be defined naturally in terms of the 
difference between the spins of the c and d fermions as follows: 

iV" - i(cV"c-dV"d) , n2 = Ni , 713 = N2 , ni = N3. (14) 

In view of the strong on-site repulsion in the cuprate problem, the SC order parameter should be defined on a bond 
connecting the c and d fermions. We introduce 

At ^ ^ct.«dt = 1 (_,t,t + ,t,t) , , , , (15) 

We can group these five components together to form a single vector tiq = (ni, n2, ns, 714, 715), called the superspin 
since it contains both swpereonducting and antiferromagnetic spin components. The individual components of the 
superspin are explicitly defined in the last parts of Eqs. f(T^ and l(T^ . 

The concept of the superspin is useful only if there is a natural symmetry group acting on it. In this case, since the 
order parameter is five dimensional, it is natural to consider the most general rotation in the five dimensional order 
parameter space spanned by Ua- In three dimensions, three Euler angles are needed to specify a general rotation. In 
higher dimensions, a rotation is specified by selecting a plane and an angle of rotation within this plane. Since there 



10 



are n{n— l)/2 independent planes in n dimensions, the group SO{n) is generated by n(n— l)/2 elements, specified in 
general by antisymmetric matrices Lab — ~Lba, with a = 1, ..,n. In particular, the 50(5) group has ten generators. 
The total spin and the total charge operators, 



(16) 



perform the rotation of the AF and SC order par ameters within each subsp ace. In addition, there are six so-called tt 
operators, first introduced by Demler and Zhang jPemler and Zhan j . Il99fj) . defined by 



(17) 



They perform the rotation from AF to SC and vice versa. These infinitesimal rotations are defined by the commutation 
relations 



(18) 



The total spin components Sa^ the total charge Q, and the six tt operators form the ten generators of the SO{b) 
group. 

The superspin order parameters n^, the associated 50(5) generators Lab, and their commutation relations can be 
expressed compactly and elegantly in terms of the SO{b) spinor and the five Dirac F matrices. The four component 
SO{5) spinor is defined by 



di 



They satisfy the usual anti-commutation relations 

{vI/t,vl/,} = 5^, , {vl/^,vl/j = {vl/t^vl,t} = 0. 
Using the ^' spinor and the five Dirac F matrices (see appendix A), we can express and Lab sts 



2 tiv'^^ ' -^af 



2 M ^ 



The Lab operators forms the 50(5) Lie algebra and satisfy the commutation relation 

[Lab, Lcd\ — —i{SacLbd + SbdLac — SadLbc — SbcLad)- 



(19) 



(20) 



(21) 



(22) 



The Ha and the operators form the vector and the spinor representations of the 50(5) group, satisfying the 
equations 



[Lab, ric 



and 



'ii^acn-b - SbcUa) 



1 



[Lab,^,]^-TZ^. 



(23) 



(24) 



If we arrange the ten operators 5q, Q and tTq into LabS by the following matrix form 



Lnh — 



( 



V Q i(4-^.) \{^l-^y) OJ 











(25) 



and group Ua as in Eqs. l(Tl| and (fT5|l . we see that Eqs. and (|^ compactly reproduces all the commutation 
relations presented previously. These equations show that Lab and ria are the symmetry generators and the order 
parameter vectors of the 50(5) theory. The commutation relation Eq. H23|l is the 50(5) generalization of the 50(3) 
communication relation as given in Eqs. © and Hll|l . 
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In systems where the unit cell naturally contains two sites, such as the ladder and the bi-layer systems, the complete 
set of operators Lab, n-a and can be used to construct model Hamiltonians with the exact 5*0(5) symmetry, as we 
will show in section IV.AI In these models, local operators are coupled to each other so that only the total symmetry 
generators, obtained as the sum of local symmetry generators, commute with the Hamiltonian. For two dimensional 
models containing only a single layer, grouping the lattice into clusters of two sites would break lattice translational 
and rotational symmetry. In this case, it is better to use a cluster of four sites forming a squar e, which does not break 
rotational symmetry and can lead naturally to the definition of a d-wave pairing operator iAltma n and Auerbach| . 
|2002; Zhang el a;,, 1999). In this case, the Lab, "-a and ^P^t operators are interpreted as the effective low energy 
operators defined on a plaquette, which form the basis for an effective low energy S0{5) theory, rather than the basis 
of a microscopic SO{5) model. 

Having introduced the concept of local symmetry generators and order parameters based in real space, we will 
now discuss definitions of these operators in momentum space. The AF and SC order parameters can be naturally 
expressed in terms of the microscopic fermion operators as 

N"^J2i+n^"^P ' = Y^'^^^^^^^''^-^' ' d{p)^co.s{p,)-cos{py), (26) 
p p 

where 11 = (tt, tt) and d{p) is the form factor for the d wave pairing operator in two dimensions. They can be combined 
into the five component superspin vector Ua by using the same convention as before. The total spin and total charge 
operator are defined microscopically as 

5„ = ^ctr% , Q-i5](ctcp-l), (27) 
p p 

and the tt operators can be defined as 

- J2gip)cl^nr''Tycl^. (28) 
p 

Here the form factor g{p) needs to be chosen ap propriately to satisfy the 50(5) commutation relation (|22|l . In the 
original formulation of the SO{5) theory, Zhang |)Zhangl Il99^ chose g{p) = d(p) . In thi s case, the SO{5) symmetry 
algebra 1)22(1 only closes approximately near the Fermi surface. Later, Henlev l|Henle vl ^23) proposed the choice 
g{p) = sgn{d{p)){this construction requires introducing form factors for the AF order parameter as well). When the 
momentum space operators Sa, Q and tt^^,, as expressed in Eq. (|27|l and (|28|l . are grouped into Lab according to Eq. 
(|25|l . the symmetry algebra (|22(l closes exactly. However, the 7r-operators are no longer short ranged. 

The 50(5) symmetry generators perform the most general rotation among the five order parameters. The quantum 
numbers of the tt operators exactly patch up the difference in quantum numbers between the AF and SC order 
parameters, as shown in the Table I. 





charge 


spin 


momentum 


internal angular momentum 


A, or ni, ns 


±2 








d wave 


N" or 712,3,4 





1 


{■K,n) 


s wave 




±2 


1 


{TT,n) 


d wave 



TABLE I Quantum numbers of the AF, the d-wave SC order parameters, and the tt operator. Since the tt operator rotates 
the AF and SC order parameters into each other, its quantum numbers patch up the difference between the AF and SC order 
parameters. 

With the proper choice of the internal form factors, the tt operators rotate between the AF and SC order parameters 
according to H18|l . Analogous to the electro-magnetic unification presented in the introduction, the tt operators 
generate an infinitesimal rotation between the AF and SC order parameters similar to the infinitesimal rotation 
between the electric and the magnetic fields generated by the Lorentz transformation. These commutation relations 
play a central role in the 50(5) theory and have profound implications on the relationship between the AF and SC 
order - they provide a basis to unify these two different types of order in a single framework. In the AF phase, the 
operator iV" acquires a nonzero expectation value, and the tt and SC operators become canonically conjugate variables 
in the sense of Hamiltonian dynamics. Conversely, in the SC phase the operator A acquires a nonzero expectation 
value, and the tt and AF operators become canonically conjugate variables. This canonical relationship is the key to 
understanding the collective modes in the 50(5) theory and in HTSC. 
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The SO{5) group is the minimal group to contain both AF and SC, the two dominant phases in the HTSC cuprates. 
However, it is p o ssible to generalize this construction so that it includes other forms of order. For example, in Ref. 
i Podolskv et all 120041). it was demonstrated how one can combine AF and triplet SC using an 5*0(4) symmetry 
i Rozhkov and MillisLl200l . Such a construction is useful for quasi-one-dimensional Bechga ard salts, which undergo 
a tra i isition from an AF insulat ing state to a triplet SC state as a function of pressure Jerome et aLlllOSOtlLee et all 
l200.tll<)97tlVu]etic et »4l2002^ . 

To define an SO{A) symmetry for a one-dimensional electron system, we introduce the total spin, total charge, and 
Q operators 

Q = i^(4,,c+,fc + clj^^c_,fc-2) 

= T ^ - <.rycl_,) . (29) 

fc 

Here Cj_ creates right/left moving electrons of momentum ±fc/ + fc. The spin operators Sa form an 5*0 (3) algebra of 
spin rotations given by the second formula of equation We can also introduce isospin SO{3) algebra by combining 
the charge with the Q operators 

/. = i(et + e), 4 = l(et-e), = q 

[Ia,Ib] = ieabJc- (30) 

Spin and isospin operators together generate an S'0(4) « 50(3) x 5*0(3) symmetry, which unifies triplet supercon- 
ductivity and antiferromagnetism. We define the Neel vector and the TSC order parameter, 

= iT^i'l-M^^^-.k + ^-^kr^c+.K) (31) 
fc 

= - Vc+,fc(T^T")c_,_fe , 
k 



and combine them into a single tensor order parameter 

/ (Re^). 
Q - (Re^)j, 
V (Re*), 



(Im*)^ \ 

(Im^)j, Ny . (32) 
(Im*), I 



One can easily verify that Qaa transforms as a vector under both spin and isospin rotations 

\ScnQhfi\ — Kal3fQb-i [Ia,Qbp] = i^abcQcp- (33) 

One dimensional electron systems have been studied extensively using bosonization and renormalization group anal- 
ysis. They have a line of phase transitions between an antiferromagnetic and a triplet superconducting phase at a 
special ratio of the forward and backward scattering amplitudes. Podolsky et al pointed out that anywhere on this 
line the O operator commutes with the Hamiltonian of the system. Hence, one finds the 5*0(4) symmetry at the 
AF/triplet SC phase boundary witho ut any fine tuning of the parameters. Consequences of this symmetry for the 

Bechgaard salts are reviewed in Ref. l|Podolskv et all\2m4\ . 

Other extensions an d gen eralizations of 5*0(5) are discussed in Ref. ijLin et al\ . Il998t iMarkiewicz and VaughnL 
hflflStlMnrakami et »/.LlT9MlNava,kll20f)nllSchii]>lll998tlWii ef'ai\ . \2m?Aif . 

B. The S0(5) quantum nonlinear a model 

In the previous section, we presented the concepts of local 50(5) order parameters and symmetry generators. 
These relationships are purely kinematic and do not refer to any particular Hamiltonian. In section IV.AI we shall 
discuss microscopic models with exact 5*0(5) symmetry, constructed out of these operators. A large class of models. 
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however, may not have SO{5) symmetry at the microscopic level, but their long distance, low energy properties may 
be described in terms of an effective 5*0(5) model. In section^ we saw that many different microscopic models 
indeed have the SO{3) non-linear a model as their universal low energy description. Therefore, in order to present a 
general theory of the AF and SC in the HTSC, we first introduce the SO{5) quantum non-linear cr model. 

The 50(5) quantum non- linear a model describes the kinetic and potential energies of coupled superspin degrees 
of freedom. In the case of HTSC cuprates, the superspin degrees of freedom are most conveniently defined on a coarse 
grained lattice, with 2a x 2a lattice spacing in units of the original cuprate lattice spacing, where every super-site 
denotes a (non-overlapping) plaquette of the original lattice (see Fig. I29f) . There are 4"* = 256 states on a plaquette in 
the original Hubbard model, but we shall retain only the 6 lowest energy states, including a spin singlet state and three 
spin triplet states at half-filling, and two paired states with two holes or two particles away from half-filling (see Fig. 
EJ. In sections IVTTl and rvTDl we will show, with numerical calculations, that these are indeed the lowest energy states 
in each charge sector. Additionally, we will show explicitly that the local superspin degree of freedom discussed in 
this section can be constructed from these six low energy states. Proposing the 50(5) quantum non-linear a model as 
the low energy effective model of the HTSC cuprates requires the following physical assumptions: 1) AF and SC and 
their quantum disordered states are the only competing degrees of freedom in the underdoped regime. 2) Fermionic 
degrees of freedom are mostly gapped below the pseudogap temperature. For a c?-wave superconductor, there are 
also gapless fermion degrees of freedom at the gap nodes. However, they do not play a significant role in determining 
the phase diagram and collective modes of the system. Our approach is to solve the bosonic part of the model first, 
and th en include gapless fermions self-consistently at a later stage (Altniau and Auerbach, 2002; Dcmlcr and ZhanJ 
ll999aD . 

From Eqs. and the discussions in sections Hit. Al we see that Lab and are conjugate degrees of freedom, very 
much similar to [q-,p] = ih in quantum mechanics. This suggests that we can construct a Hamiltonian from these 
conjugate degrees of freedom. The Hamiltonian of the 50(5) quantum non-linear a model takes the following form 

H=^J2 Llbi^) + f E ^a{xW{x') + Bab{x)Lab{x) + ^ V{n{x)l (34) 

a:,a<& {x.x') x,a<b x 

where the superspin Ua vector field is subjected to the constraint 

nl = 1. (35) 

This Hamiltonian is quantized by the canonical commutation relations 122|) and 123|l . Here, the first term is the kinetic 
energy of the 50(5) rotors, where x has the physical interpretation of the moment of inertia of the 50(5) rotors. 
The second term describes the coupling of the 50(5) rotors on different sites through the generalized stiffness p. The 
third term introduces the coupling of external fields to the symmetry generators, while V{n) includes anisotropic 
terms which break the 50(5) symmetry to the conventional 50(3) x U{1) symmetry. The 50(5) quantum non-linear 
a model is a natural combination of the 50(3) non-linear a model describing the AF Heisenberg model and the 
quantum XY model describing the SC to insulator transition. If we restrict the superspin to have only components 
a — 2,3,4, then the first two terms describe the symmetric Heisenberg model, the third term describes the coupling 
to a uniform external magnetic field, while the last term can represent easy plane or easy axis anisotropy of the Neel 
vector. On the other hand, for a = 1,5, the first term describes Coulomb or capacitance energy, the second term is 
the Josephson coupling energy, while the third term describes coupling to an external chemical potential. 

The first two terms of the 50(5) model describe the competition between the quantum disorder and classical 
order. In the ordered state, the last two terms describe the competition between the AF and SC order. Let us 
first consider the quantum competition. The first term prefers sharp eigenstates of the angular momentum. On an 
isolated site, C = J2 ^ab Casimir operator of the 50(5) group in the sense that it commutes with all the 

50(5) generators. The eigenvalues of this operator can be determined completely from group theory - they are 0, 4, 6 
and 10, respectively, for the 1 dimensional 50(5) singlet, 5 dimensional 50(5) vector, 10 dimensional antisymmetric 
tensor and 14 dimensional symmetric, traceless tensors, respectively. Therefore, we see that this term always prefers 
a quantum disordered 50(5) singlet ground state, which is also a total spin singlet. In the case where the effective 
quantum non-linear a model is constructed by grouping the sites into plaquettes, the quantum disordered ground state 
corresponds to a plaquette "RVB" state, as depicted in Figs. andll2b,. This ground state is separated from the first 
excited state, the five fold 50(5) vector state, by an energy gap of 2/x- This gap will be reduced when the different 
50(5) rotors are coupled to each other by the second term. This term represents the effect of stiffness, which prefers 
a fixed direction of the Ua vector to a fixed angular momentum. This competition is an appropriate generalization of 
the competition between the number sharp and phase sharp states in a superconductor and the competition between 
the classical Neel state and the bond or plaquette singlet state in the Heisenberg AF. The quantum phase transition 
occurs near — 1- 
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In the classically ordered state, the last two anisotropy terms compete to select a ground state. To simplify the 
discussion, we first consider the following simple form of the static anisotropy potential: 

V{n) ^ -g{jil+Til + jil). (36) 

At the particle-hole symmetric point with vanishing chemical potential -B15 = fi = 0, the AF ground state is selected 
by (7 > 0, while the SC ground state is selected by 5 < 0. 5 = is the quantum phase transition point separating 
the two ordered phases. This phase transition belongs to "class ^" in the classification scheme of section Ull and is 
depicted as the "Al" transition line in Fig. ^1 This point has the full quantum SO{5) symmetry in the model 
described above. 

However, it is unlikely that the HTSC cuprates can be close to this quantum phase transition point. In fact, we 
expect the anisotropy term g to be large and positive, making the AF phase strongly favored over the SC phase at 
half-filling. However, the chemical potential term has the opposite, competing effect and favors SC. We can observe 
this by transforming the Hamiltonian into the Lagrangian density in the continuum limit: 

C = ^tOab{x,tf + ^(dku'^ix^t))^ - (37) 

where 

uJab ^ riaidtrib - iBbcHc) ~ {a b) (38) 

is the angular velocity. We see that the chemical potential enters the Lagrangian as a gauge coupling in the time 
direction. Expanding the first term in the presence of the chemical potential fi — i?i5, we obtain an effective potential 

K//(n) = Vin) - ^^(n? + nj), (39) 

from which we see that the bare V term competes with the chemical potential term. For ^ < — \f9Jx1 AF 
ground state is selected, while for [i > [ic, the SC ground state is realized. At the transition point - even though 
each term strongly breaks S'0(5) symmetry - the combined term gives an effective static potential which is 50(5) 
symmetric, as we can see from (|39|) . This quantum phase transition belongs to "class i?" in the classification scheme 
of section |H] A typical transition of this type is depicted as the "i31" transition line in Fig. El Even though the 
static potential is SO(h) symmetric, the full quantum dynamics is not. This can be seen most easily from the time 
dependent term in the Lagrangian. When we expand out the square, the term quadratic in [i enters the effective 
static potential in Eq. (|39|l . However, there is also a /i-dependent term involving a first order time derivative. This 
term breaks the particle hole symmetry and dominates over the second order time derivative term in the n\ and 715 
variables. In the absence of an external magnetic field, only second order time derivative terms of ^2.3,4 enter the 
Lagrangian. Therefore, while the chemical potential term compensates the anisotropy potential in Eq. (|39|l to arrive 
at an 50(5) symmetric static potential, its time dependent part breaks the full quantum 50(5) symmetry. This 
observation leads to the concept of the projected or static 50(5) symmetry. A model with projected or static 50(5) 
symmetry is described by a quantum effective Lagrangian of the form 

^ = \ X! (dtna f ~ xKiT-idins - n^dtui) -Veff{n), (40) 

a=2,3,4 

where the static potential Veff is 50(5) symmetric. 

We see that "class A" transition from AF to SC occurs at a particle hole symmetric point, and it can have a full 
quantum 50(5) symmetry. The "class B" transition from AF to SC is induced by a chemical potential; only static 
50(5) symmetry can be realized at the transition point. The "class A" transition can occur at half-filling in organic 
superconductors, where the charge gap at half-filling is comparable to the spin exchange energy. In this system, the 
AF to SC transition is tuned by pressure, where the doping level and the chemical potential stay fixed. The transition 
from the half-filled AF state to the SC state in the HTSC cuprates is far from the "class A" transition point, but static 
50(5) symmetry can be realized at the chemical potential induced transition. However, as we shall see in section lTV.BI 
there are also Mott insulating states with AF order at fractional filling factors, for instance, at doping level x = 1/8. 
The insulating gap is much smaller at these fractional Mott phases, and there is an effective particle-hole symmetry 
near the tip of the Mott lobes. For these reasons, "class A" transition with the full quantum 50(5) symmetry can be 
realized again near the tip of fractional Mott phases, as in organic superconductors. Transitions near the fractional 
Mott insulating lobes are depicted as the "A2" and "B2" transitions in the global phase diagram (see Fig. I13|) . In 
this case, a transition from a fractional Mott insulating phase with AF order to the SC state can again be tuned by 
pressure without changing the density or the chemical potential. 
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The 5*0(5) quantum nonlinear a model is constructed from two canonically conjugate field operators Lab and Ua- 
In fact, there is a kinematic constraint among these field operators. In the case of the Heisenberg model, the total 
spin operator and the AF Neel order parameter satisfy an orthogonality constraint, as expressed in Eq. 0. The 
50(5) generalization of this constraint can be expressed as follows: 

LabTlc + LbcUa + L^ani, = 0. (41) 

This identity is valid for any triples a, b and c, and can be easily proven by expressing Lab — naPb — nbPa, where pa 
is the conjugate momentum of n^. Geometrically, this identity expresses the fact that Lab generates a rotation of the 
Ua vector. The infinitesimal rotation vector lies on the tangent plane of the four sphere 5 ^, as defined b y Eq. 135|l . 
and is therefore orthogonal to the Ua vector itself. Extending this geometric proof, Wegner llWegneill200(i has shown 
that the 50(5) orthogonality relation also follows physically from maximizing the entropy. Taking the triple a, b, c 
to be 2,3,4, and recognizing that La/s = eap-yS-y, this identity reduces to the 50(3) orthogonality relation in Eq. 
©. This 50(5) identity places a powerful constraint on the expectation values of various operators. In particular, 
it quantitatively predicts the value of the tt order parameter in a mixed state between AF or SC. For example, let's 
take the a, b, c triple to be 1, 2, 5. Eq. (|^ predicts that 

ii5"2 + ^52"! + L2i"-5 = ^ {L25) = {Iimr^) = ^^^-^ , (42) 

where we chose the SC phase such that (n^) — 0. Here, Q — (Li^) measures the hole density. Since these four 
expectation values can easily be measured numerically and, i n principle, experirn entally, this relationship can be 
tested quantitatively. Recently, Ghosal, Kallin and Berlinskvll Ghosal et a/.lT2002(l tested this relationship within 
microscopic models of the AF vortex core. In this case, AF and SC coexist in a finite region near the vortex core, so 
that both (ni) and (71-5) are non- vanishing. They found that the 50(5) orthogonality constraint is accurately satisfied 
in microscopic models. 

In this section, we presented the 50(5) quantum non-linear a model as a heuristic and phenomenological model. 
The key ingredients of the model are introduced by observing the robust features of the phase diagram and the low 
energy collective modes of the HTSC cuprate system. This is the "top-down" approach discussed in the introduction. 
In this sense, the model has a general validity beyond the underlying microscopic physics. However, it is also useful 
to derive such a model directly from microscopic electronic models. Fortunately, this "bottom-up" approach agrees 
with the phenomenological approach to a large extent. A rigorous derivation of this quantum non- linear a model 
from an 50(5) symmetric microscopic model on a bi-layer system will be given in section FV.AI while an approximate 
derivation from the "realistic" microscopic t — J and Hubbard model will be given in section IV.DI 



C. The projected S0(5) model with lattice bosons 

In the previous section, we presented the formulation of the 50(5) quantum nonlinear a model. This model is formu- 
lated in terms of two sets of canonically conjugate variables - the superspin vector Ua and the angular momentum Lab- 
The two terms which break the full quantum 50(5) symmetry are the anisotropy term, g, and the chemical potential 
term, /i. Therefore, this model contains high energy modes, particularly excitations of the order of the Mott insulatirig 
gap at half-filling. For this reason, Greiter( Greiter, 1997) and Baskaran and Anderson! Baskaran and Anderson, 1993) 
questioned whether the effective 50(5) symmetry can be implemented in the low energy theory. In the previous 
section, it was shown that these two symmetry breaking terms could cancel each other in the static potential and 
the resulting effective potential could still be 50(5) symmetric. It was also pointed out that the chemical potential 
term breaks the 50(5) symmetry in the dynamic or time- dependent part of the effective La grangian. In response 
to these observations, Zhang et al. constructed the projected 50(5) models (|Zhang et all\l99^ . which fully project 
out the high energy modes, and obtained a low energy effective quantum Hamiltonian, with an approximately 50(5) 
symmetric static potential. 

The first step is to perform a transformation from the Ua and Lab coordinates to a set of bosonic operators. We 
first express the angular momentum operator as 

Lab = naPb - ntpa, (43) 
where Pa is the canonical momentum conjugate to Ua, satisfying the Heisenberg commutation relation: 



[na,Pb] = iSab- 



(44) 
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Furthermore, we can express the canonical coordinates and momenta in terms of the boson operators as 

""""75^*""^*^^ ^""^71^^""*^^' ^^^^ 

where the boson operators satisfy the commutation relation 

[ta,tl]=5ab. (46) 

and the (half-filled) ground state is defined by ta\^) = 0. There are five boson operators, ta — ^2,^3,^4 are the boson 
operators for the magnetic triplet excitations, also called the magnons, while 

^^ = ^^*^ + ^P^ ^^^^ 

are the linear combinations of the particle pair (tp) and hole pair (th) annihilation operators. In the 5*0(5) quantum 
non-linear a model formulation, there is an infinite number of bosonic states per site. However, due to the first term in 
Eq. (|34|) (the angular momentum term), states with higher angular momenta or, equivalently, higher boson number, 
are separated by higher energies. Therefore, as far as the low energy physics is concerned, we can restrict ourselves 
to the manifold of six states per site, namely the ground state \n) and the five bosonic states t\\^)- This restriction 
is called the hard-core boson constraint and can be implemented by the condition = 0. Within the Hilbert 

space of hard-core bosons, the original SO{5) quantum non-linear a model is mapped onto the following hard-core 
boson model: 

x,Q=2,3,4 2;, 1=1, 5 X 

- n-a{x)na{x) - nt{x)ni{x), (48) 

{xx') {xx'} 

where Ac — 2/x ~ and = 2/x — g are the creation energies for the charge pairs and the triplet magnons, /i 

is the chemical potential, and Jc and Js are the exchange energies for SC and AF, respectively. In the isotropic case, 
they are taken to be p in the second term of Eq. H34|) . Expressing rii and Ua in terms of the bosonic operators, we 
see that the Jc and Jg terms describe not only the hopping, but also the spontaneous creation and annihilation of the 
charge pairs and the magnons, as depicted in Fig0] 

When Ac = A^, Jc = Js and /i = 0, the model (|48ll has an exact quantum 5*0(5) symmetry. In this case, the energy 
to create charge excitations is the same as the energy to create spin excitations. This situation could be realized in 
organic and heavy fermion superconductors near the AF phase boundary or the HTSC near commensurate doping 
fractions such as a: = 1/8, as we shall see in section HV.BI However, for HTSC systems near half-filling, the energy to 
create charge excitations is much greater than the energy to create spin excitations, i.e. Ac 2> A^. In this case, the 
full quantum 50(5) symmetry is broken but, remarkably, the effective static potential can still be SO{5) symmetric. 
This was seen in the previous section by the cancellation of the anisotropy potential g term by the chemical potential 
fi term. In a hard-core boson model (|48l) with Ac 3> A^, a low energy effective model can be derived by retaining 
only the hole pair state while projecting out the particle pair state. This can be done by imposing the constraint 

tl{x)\n)^0 (49) 

at every site x. The projected Hamiltonian takes the form 

X X 

- JsYj '^a{x)na{x') - JcYj ni{x)ni{x'), (50) 

{xx') {xx') 

where Ac — Ac — /i. The Hamiltonian (|50|l has no parameters of the order of U. To achieve the static 50(5) symmetry, 
we need A^ ~ Ac and Js ^ Jc. The first condition can always be met by changing the chemical potential, whereas the 
second one requires certain fine tuning. However, as we discuss in Section VD (see Fig. I31II . this condition emerges 
naturally when one derives the model (|50|l from the Hubbard model in the relevant regime of parameters. 
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The form of the projected SO{5) Hamihonian hardly changes from the unprojected model 1)4811 . but the definition 
of ni and is changed from 

"1 = -^{h + t\) = ^{th + tp + tl + tl) 



to 



"i = ^(</. + 4) ^5 = ^(th-4)- (52) 

From Eq. (|51l) . we see that ni and 715 commute with each other before the projection. However, after the projection, 
they acquire a nontrivial commutation relation, as can be seen from Eq. (|51: 



[ni,n5] = i/2. (53) 

Therefore, projecting out doubly occupied sites, commonly referred to as the Gutzwiller projection, can be analytically 
implemented in the S0{5) theory by retaining the form of the Hamiltonian and changing only the commutation 
relations. 

The Gutzwiller projection implemented through the modified commutation relations between ni and is formally 
similar to the projection onto the lowest Landau level in the physics of the quantum Hall effect. For electrons moving 
in a 2D plane, the canonical description involves two coordinates, X and Y, and two momenta, Px and Py- However, 
if the motion of the electron is fully confined in the lowest Landau level, the projected coordinate operators become 
non-commuting and are given by [X, Y] — iZg, where Iq is the magnetic length. In the context of the projected 5*0(5) 
Hamiltonian, the original rotors at a given site can be viewed as particles moving on a four dimensional sphere S'^, as 
defined by Eq. (|35|l . embedded in a five dimensional Euclidean space. The angular momentum term ■^L'^t describes 
the kinetic motion of the particle on the sphere. The chemical potential acts as a fictitious magnetic field in the 
(^1,715) plane. In the Gutzwiller-Hubbard limit, where Ac 3> A^, a large chemical potential term is required to reach 
the limit Ac ~ A^. The particle motion in the (771,715) plane becomes quantized in this limit, as in the case of the 
quantum Hall effect, and the non-commutativity of the coordinates (771,775) given by Eq. H53|l arises as a result of 
the projection. The projection does not affect the symmetry of the sphere on which the particle is moving; however, 
it restricts the sense of the kinetic motion to be chiral, i.e., only along one direction in the (771,775) plane. (See Fig. 
ISJ. In this sense, the particle is moving on a chiral SO{5) symmetric sphere. The non-commutativity of the (771,775) 
coordinates is equivalent to the effective Lagrangian (see Eq. (|40|) of section lTll.Ij|l containing only the first order time 
derivative. In fact, from Eq. (|40|l . we see that in this case the canonical momenta associated with the coordinates 771 
and 775 are given by 

Pi = ^ XfJ-ns , P5 = ^ -XM"-i- 54) 

0771 0775 

Applying the standard Heisenberg commutation relation for the conjugate pairs (771, pi), or (775,^5) gives exactly the 
quantization condition H53|l . Note that in Eq. (|54|l xf^ plays the role of the Planck's constant in quantum mechanics. 
We see that the projected SO{5) Hamiltonian (|50l) subjected to the quantization condition (|53|l is fully equivalent to 
the effective Lagrangian Eq. H4U|) , discussed in the last section. 

Despite its apparent simplicity, the projected 50(5) lattice model can describe many complex phases, most of 
which are seen in the HTSC cuprates. These different phases can be described in terms of different limits of a single 
variational wave function of the following product form: 

l*> = n(cosf(a;) + sme{x){m^{x)ti{x) + A{x)tl{x)))\n) . (55) 

X 

where the variational parameters 777ct(a;) should be real, while A(x) is generally complex. The normalization of the 
wave function, (^'l^') = 1, requires the variational parameters to satisfy 

J2Mx)\' + \Aixr = l. (56) 

a 

Therefore, we can parameterize them as |r77c(a;)p = cos'^4){x) and |A(a;)p = sin'^(j){x), which is similar to the 50(5) 
constraint introduced in Eq. (|35|l . The expectation values of the order parameters and the symmetry generators in 
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this variational state are given by 



and 



(*|na(a;)|*) = —=sin20{x)Re{ma{x)) 
v2 



-sin29{x)Re{A{x)) 



(^'|n5(a;)|5') = -sin20{x)Im{A{x)) 



(ml{x)th{x)\'^) ^ smH{x) |A(a;)|2 



(57) 



mQix)\^) 

(^IS-^I*) = -{■^\ie°'^^tl{x)t^{x)\'f) ^ -ie^f^'^sin'^e ml{x)m^{x) 



(*|7r„(x)|^') 



sin^ 9 TO* (a;)A(x) 



(58) 



Initially, we restrict our discussions to the case where the variational parameters are uniform, describing a trans- 
lationally invariant state. Evaluating different physical operators in this state gives the result summarized in the 
following table: 







charge Q 


spin S 


AF order (ua) 


SC order (ui) 


TT order (ttq) 


(a) 


"RVB" state: sinO = 

















(b) 


Magnon state: cos9 = and sincf} = 





1 











(c) 


"Hole pair" state: cos6 = and cos(j) = 


-2 














(d) 


AF state: < sin9 < 1 and sm<j!> = 





indefinite 


/o 








(e) 


SC state: < sin9 < 1 and coscfi = 


indefinite 













(f) 


Mixed AF/SC state: < sinO < 1 and < sm<j!> < 


indefinite 


indefinite 










TT state: cos6 = and < cos<jf> < 1 


indefinite 


indefinite 











TABLE II Physical properties of various plaquette states classified according to the SO{5) order parameters and symmetry 
generators. 



As we can see, this wave function not only describes classically ordered states with spontaneously broken symmetry, 
but also quantum disordered states which are eigenstates of spin and charge. Generally, and A^ favor quantum 
disordered states, while Jc and Js favor classically ordered states. Depending on the relative strength of these 
parameters, a rich phase diagram can be obtained. 

The phase diagram of the projected SO{5) model with Jc = 2Js = J is shown in Fig. [T] Changing the chemical 
potential modifies Ac and traces out a one-dimensional path on the phase diagram. Along this path the system goes 
from the AF state to the uniform AF/SC mixed phase and then to the SC state. The mixed phase only corresponds 
to one point on this trajectory (i.e. a single value of the chemical potential Hc), although it covers a whole range of 
densities < p < pc- This suggests that the transition may be thought of as a first order transition between the AF 
and SC phases, with a jump in the density at p,c. The spectrum of collective excitations shown in Fig. |S1 however, 
shows that this phase diagram also has important features of two second order phase transitions. The energy gap 
to S' = 1 excitations inside the SC phase vanishes when the chemical potential reaches the critical value p,c- Such a 
softening should not occur for the first order transition but is required for the continuous transition into a state with 
broken spin symmetry. This shows that models with the iS'0(5) symmetry have intrinsic fine-tuning to be exactly at 
the border between a single first order transition and two second order transitions; in subsequent sections this type 
of transition shall be classified as type 1.5 transition. Further discussion of the phase diagram of the projected SO{5) 
model is given in section VC. Note that effective bosonic H amiltonians similar to (|50|l have also been considered in 
Refs. llvan Duin and Zaanenll2000l:lPark and Sachdevll200lD . 



IV. THE GLOBAL PHASE DIAGRAM OF S0(5) MODELS 

A. Phase diagram of the classical model 

The two robust ordered phases in the HTSC cuprates are the AF phase at half-filling and the SC phase away 
from half-filling. It is important to ask how these two phases are connected in the global phase diagram as different 
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tuning parameters such as the temperature, the doping level, the external magnetic field, etc, one varied. Analyzing 
the 5*0(5) quantum no nlinear a mod el, Zhang has classified four generic types of phase diagrams, presented as Fig. 
(lA)-(lD) in reference ((ZhanJ • the next section we are going to investigate the zero temperature phase 

diagram where the AF and the SC phases are connected by various quantum disordered states, often possessing charge 
order. In this section, we first focus on the simplest possibility, where AF and SC are the only two competing phases 
in the problem, and discuss the phase diagram in the plane of temperature and chemical potential, or doping level. 

Let us first discuss the general properties of phase transition between two phases, each characterized by its own 
order parameter. In particular, we shall focus on the phenomenon of the enhanced symmetry at the multicritical point 
at which physically different static correlation functions show identical asymptotic behavior. In the case of CDW to 
SC transition discussed in section IHl the CDW is characterized by an Ising like Z-i order parameter, while the SC is 
characterized by a V{V) order parameter. In the case of AF to SC transition, the order parameter symmetries are 
50(3) and C/(l), respectively. Generically, the phase transition between two ordered phases can be either a single 
direct first order transition or two second order phase transitions with a uniform mixed phase in between, where both 
order parameters are non-zero. This situation ca n be understood easily by describing the competition in terms of a 
LG functional of two competing order Darameters l)Kosterlitz et all . Il97fi[l . which is given by 

F = ri</.2 + r2(\)\ + u^<\)\ + U2<t>\ + 2ui2</'?0^ (59) 

Here, 01 and 02 are vector order parameters with A^i and Ni components, respectively. In the case of current interest, 
iVi = 2 and iV2 = 3 and we can view = nl+nl as the SC component of the superspin vector, and 02 = n^ + n^ + nl 
as the AF component of the superspin vector. These order parameters are determined by minimizing the free energy 
F, and are given by the solutions of 

2ui02 + 2ui202 + ''I = , 2mi20? + 2U202 + ^^2 = 0. (60) 

These equations determine the order parameters uniquely, except in the case when the determinant of the linear 
equations vanishes. At the point when 



and 



the order parameters satisfy the relations 



UiU2 = Mi2 (61) 



(62) 



y/u2 y/ul 



~ const, (63) 



but they are not individually determined. In fact, with the re-scaling 0^ = 0i/y^ and 02 — (fy^/ \/ui^ the free energy 
is exactly 50(5) symmetric with respect to the scaled variables, and Eq. (|63|l becomes identical to Eq. H35I) in the 
50(5) case. Since the free energy only depends on the combination -I- 0|, one order parameter can be smoothly 
rotated into the other without any energy cost. Equation H61|) is the most important condition for the enhanced 
symmetry. We shall discuss extensively in this paper whether this condition is satisfied microscopically or close to 
some multi-critical points in the HTSC cuprates. On the other hand, equation (|62|l can always be tuned. In the case 
of AF to SC transition, the chemical potential couples to the square of the SC order parameter, as we can see from 
Eq. H39I) . Therefore, ri can be tuned by the chemical potential, and equation H62|l defines the critical value of the 
chemical potential /ic at which the phase transition between AF and SC occurs. At this point, the chemical potential 
is held fixed, but the SC order parameter and the charge density can change continuously according to Eq. H63() . Since 
the free energy is independent of the density at this point, the energy, which differs from the grand canonical free 
energy by a chemical potential term can depend only linearly on the density. The linear dependence of the energy 
on doping is a very special, limiting case. Generally, the energy versus doping curve would either have a negative 
curvature, classified as "type 1," or a positive curvature, classified as "type 2" (see Fig. |^). The special limiting 
case of "type 1.5" with zero curvature is only realized at the 50(5) symmetric point. The linear dependence of the 
ground state energy of a uniform AF/SC mixed state on the density is a crucial test of the 50(5) symmetry which 
can be performed numerically, as we shall see in section FV. Bl and IV.CI The constancy of the chemical potential and 
the constancy of the length of the 5*0(5) superspin vector H63|l as a function of density can be tested experimentally 
as well, as we shall discuss in section IV.BI 

The constancy of the chemical potential as a function of the density in a uniform system is a very special situation 
which only follows from the enhanced symmetry at the phase transition point. In a system with phase separation. 
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the chemical potential is also independent of the total density, but the local density is non-uniform. The two phases 
are generally separated by a domain wall. The SO{b) symmetric case can be obtained from the phase separation 
case in the limit where the width of the domain wall goes to infinity and a uniform state is obtained. This situation 
can be studied analytically by solving Eq. KiOII . Defining the parameters that characterize the deviation from the 
symmetric point as w — ui2 — y/uiU2 and g — ("^757 ~ "^^J-^/^' obvious that the phase transition between the 
two forms of order is tuned by g, while w determines the nature of the phase transition. The phase diagram in the 
{g,w) plane is shown in Fig. 1^;. For w > 0, the two ordered phases are separated by a first order line. This type 
of transition is classified as "type 1." On the other hand, when w < 0, the two ordered phases are separated by 
two second order phase transition lines with an intermediate mixed phase where two orders coexist, i.e. ^ 

and {(f>2) 0. This type of transition is classified as "type 2." The limiting "type 1.5" behavior corresponds to the 
symmetric point w — 0. Approaching this point from w > 0, the first order transition becomes weaker and weaker and 
the latent heat associated with the first order transition becomes smaller and smaller. Therefore, the symmetric point 
can be viewed as the end point of a first order transition. On the other hand, approaching the symmetric point from 
w < 0, the width of the intermediate mix phase becomes smaller and smaller, until the two second order transition 
lines merge into a single transition at w = 0. From the above discussion, we learn an important lesson: the phase 
transition between two ordered phases can be either a direct first order transition or two second order transitions with 
an intermediate mixed phase. Furtherm ore, the symrnetric p oint realizes a lim iting behavior which sep arates these 
two scenarios. Balents, Fisher and NavakfB alents et aZ.l . ll998|) . Lee and Kivelson l)Lee and Kivelso'n !.'2003') pointed out 
that the "type 1" and "type 2" transitions of a Mott insulator induced by varying the chemical potential is analogous 
to the two types of superconductor to normal state transitions induced by a magnetic field. The magnetic field induces 
a direct first order transition from the SC state to the normal state in "type 1" superconductors, while it induces 
two second order transitions with an intermediate mix state in the "type 2" superconductors. Indeed, the limiting 
"type 1.5" behavior separating the "type 1" and the "type 2" superconductors also has a special symmetry, where 
the Bogomol'nyi's bound for the vortex is satisfied as an equality. We note that recent work of Senthil et al discussed 
an alternative scenario for a direct second order transition between two phases with different order parameters and 
without a higher symm etry at the transitio n point. This was achieved by having fractionalizated excitations at the 
quantum critical point ijSenthil et a/.ll200^ . 

Let us now turn to the finite temperature phase transitions. In D = 3, finite temperature phase transitions 
associated with continuous symmetry breaking are possible. Therefore, the order parameters (jfi and (f>2 can each 
have their own phase transition temperature, Tc and Tjv. The interesting question is how these two second order 
lines merge as one changes the parameter g or, equivalently, the chemical potential fi, which interchanges the relative 
stability of the two ordered phases. There are two generic possibilities. The "type 1" phase diagram is shown in Fig. 
llOb . where the two second order phase transition lines intersect at a bi- critical point, Ti,c, which is also the termination 
point of the first order transition line separating the two ordered phases. This type of phase diagram is realized for 
U12 > y/uiU2. The first order transition at separates the AF and SC states with different densities; therefore, the 
T versus 5 phase diagram shown in Fig. llOb contains a region of phase separation extending over the doping range 
Q < 5 < 5c- The "type 2" phase diagram is shown in Fig. llOb . where Tc and T^r intersect at a tetra- critical point, 
below which a uniform AF/SC mixed phase separates the two pure phases by two second order transition lines. This 
type of phase diagram is reahzed for U12 < ^uiU2. 

In contrast to the conventional superconductors with a long coherence length, the HTSC cuprates have a short co- 
herence length and a large Ginzburg region. Thus, they have the possibility of observing non-trivial critical behaviors. 
An interesting point concerns the symmetry at the multi-critical point where T/v and Tc (or, more generally, Ti and 
T2) intersect. At the multi-critical point defined by ri — r2 — 0, the critical fluctuations of the order parameters cou- 
ple to each other and renormalize the coefficients of the fourth order terms ui, U2 and U12. There are several possible 
fixed points. The symmetric fixed point, also known as the Heisenberg fixed point, is characterized by u\ = U2— u\2- 
The 0{Ni) X 0{N2) symmetry is enhanced at this point to the higher 0{Ni + N2) symmetry. Another fixed point, 
called the biconical tetra-critical point in the literature, has non- vanishing values of u\ , U2 and u\2 at the fixed point 
which deviates from the 0{Ni -\- N2) symmetry. The third possible fixed point is the decoupled fixed point, where 
U12 = and the two order parameters decouple from each other at the fixed point. The relative stability of these three 
fixed points can be studied analytically and numerically. The general picture is that there are two critical values, Nc 
and N^. For Ni + N2 < Nc, the symmetric bi-critical point is stable, for Nc < Ni + N2 < N'^, the biconical point is 
stable, while for A^i -I- N2 > Nl, the decouple d point becomes stable. The renormalization group calculations based 
on the 4 — e expansion ljKosterlitz et al] . ll976fl places the value of Nc close to 4 and the value of iV^ close to 11. The 
RG flow diagram is shown in Fig. ^2 for the case of A^i = 3 and N2 = 2. Initially, all RG trajectories flow towards 
the symmetric fixed point. The manner in which the trajectories diverge close to the symmetric point depends on 
the values of the initial parameters. The trajectories flow to the symmetric point when u\2 = uiU2, to the biconical 
point when < 'i*i'"2, and flow outside of the regime of weak coupling RG analysis when U12 > U1U2. In the case 
of competition between AF and SC, iV = A^i -f- iV2 = 5 is very close to Nc, leading to two important consequences. 
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First, the biconical point breaks the SO{5) symmetry weakly. The value of the interaction parameters at the biconical 
fixed point is given by (u^; W12) = 27r^e(0.0905; 0.0847; 0.0536). Extrapolating to e = 1 gives the root mean square 
deviation from the symmetric 5*0(5) point to about 26%, indicating weak SO{5) symmetry breaking. The second 
consequence is that the critical exponent associated with the flow away from the symmetric, 50(5) point is extremely 
slow. The first loop 4 — e expansion gives the value of 1/13 for the exponent associated with the flow away from the 
symmetric point. To get an estimate of the order of magnitude, we take the initial value of the scaling variable taking 
the flow away from the 5*0(5) fixed point to be 0.04. This value is obtained by considering the quantum corrections 
associated with a projected 50(5) model (Arrigoni and I lankc . 2000) . In this case, the significant deviation away 
from the symmetric point can only be observed when the reduced temperature is t = (T — Ti,c)/Ti,c ~ 10~^^, making 
the departure away from the SO{5) symmetric point practically unobservable. Indeed, numerical simulations of the 
5*0(5) models presented in section HV.CI are consistent with the 50(5) symmetric behavior in a wide range of tem- 
peratures and in very large systems. However, it should be noted that they do not prove the ultimate stability of the 
symmetric point. 

The qu estion of the stability of the 5*0(5) symmetric bicritical point has been raised and discussed extensively in 
litcraturc(AharonvV2002': Arrieoni and Hankc, 2000; Burgess et a/..'l998l: ICalabrese et 0111200 3': Hu and Zhang, 2000; 
IHu , 2001: Jostingmeier et al . 2003: Murakami and Nagaosa. ,2000). Because the possible flow away from the bicritical 
point is extremely slow, experimental and numerical observation of the 5*0(5) symmetric bi-critical behavior should 
be possible in a wide range of temperatures, if the starting microscopic parameters are already close to the symmetric 
point = The 5*0(5) symmetric bi-critical point has a distinct set of critical exponents, summarized in 

Ref. ijHu and Zhan g. 2000), which can be distinguished experimentally from the usual 50(3) and U{1) behavior. In 
this sense, the experimental observation of the bi-critical behavior would demonstrate that the microscopic model of 
the HTS C cuprates is close to the 50( 5) symmetry. In section IVII.DI we shall discuss the analysis of Murakami and 
Nagaosa ljMurakami and Nagaosal l2000|) showing bi-critical scaling behavior in the k-BEDT organic superconductors. 
If the microscopic paramete rs are far from the the symmetric point U12 = ^1*1^2, other critical behaviors could be 
observed. Aharonv l) Aharon v . 2002) proposed the decoupled tetra-critical fixed point with = 0. As previously 
discussed , this critical point can be observed in experiments only if the microscopic value of u^2 is already close 
to zero (due to the extremely slow flows of parameters). For the HTSC cuprates, the AF vortex core experiments 
discussed in section IVlI. Al clearly show that the AF and SC order parameters are strongly and repulsively coupled 
with w^2 > 0- Therefore, the decoupled fixed point is unlikely to be relevant for these materials. However, this 
behavior could be realized in som e heavy fermion syste ms where different b ands are responsible for the AF and SC 
phases separately. Kivelson et al l|Kivelson et all l200l|) and Calabrese et al ljCalabrese et ail l2003|) also considered 
the possibility of tri-critical points, where some of the quartic terms ui,U2,wi2 become negative and the sixth order 
terms become important. In this case, the phase diagram could have to pologies different from those listed here, and 
the readers are referred to the more extensive discussions given in Ref. ijKivelson et all\200]^ . especially Fig. Ic and 
Id of that reference. Negative values of the quartic coefficient in the free energy H59|l may come from the runaway 
flows shown in Fig. ^] A multi-critical point most closely related to t he bi-critical point is the b i conica l tetra-critical 
point. Its relevance to the HTSC cuprates has been discussed in Ref. ljZhaneL Il997t IZhang et a^J . l2002() . 

B. Phase diagram of the quantum model 

Having discussed the flnite temperature phase diagram of the classical model, we now present the global phase 
diagram of the quantum model at z ero temperatu re. The quantum phase transitions in the 5*0(5) model were discussed 
in Fig. (IC) and (ID) in re ference llZhand.ll99'<|). The quantum critical behavior of the 5*0(5) models have als o been 
studied extensively in Ref. ijKopec and Zalesk i'.'2001. 2003: Za lcski and Koped . 12000: Zal eski and Kope c. 2000). This 
section extends the original analysis to include quantum disordered states with inhomogeneous charge distributions. 
The analysis carried out in this section is based on the bosonic projected 50(5) mo del, which bears great similarities 
to the phase diagrams of the hardcore boson model studied extensively in Ref. (Bernardet et al, '200^;'Bruder et all 
[1993; Fisher et al, 1989; He bert_gi a l. 2002; Pich and Frey, 1998: Scalettar e t al, 1995; vanOtterlo e t aL, 1995). The 
iterative construction of the global phase diagram of the 5*0(5) model is also inspired b y the global phase diagram of 
the quantum Hall effect constructed by Kivelson, Lee and Zhang (iKivelson et J . 1199^ . 

The projected 50(5) model given in Eq. H50() contains the creation energy and the hopping process of the magnons 
and hole pairs. The variational wave function for this model has the general form given in (|55|) . with variational 
parameters 0{x), ma{x) and A(a;) = mi(x) -I- imz{x). The expectation value of the energy in this state is given by 



{^H\^) = E{e{x),ma{x)) 





a;a;';a— 2,3,4 



xx' 
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+As ^ siTi^ 0(x)ml{x) + Ac sin^ e{x)mf{x). (64) 

a:;a— 2,3,4 1,5 

The variational minimum is taken with respect to the normahzation condition (|56(l . In the regime when the quantum 
fluctuations are small, 0{x) can be taken to be fixed and uniform. In this case, the variational energy is nothing 
but the energy functional of a classical, generally anisotropic SO{5) rotor model, which has been studied extensively 
numericallv l)HuLl200l[l . At the point Jc = 2Js and Ac = in parameter space, this rotor model is 50(5) symmetric 
at the classical level. However, unlike the classical 50(5) rotor, the projected SO{^) model also contains quantum 
fluctuations and quantum disordered phases. T he phase diagram of the projected 50(5) model has been studied 
extensively by Quantum Monte Carlo simulations ljChen et a/J . l2003aHDorneich et a/] . l2002btlJostingmeier et aILl2003t 
iRieraL 1200231^ ). and the results will be reviewed in details in section lIV.(/l When the quantum fluctuations are not 
strong enough to destroy classical order, the general topology of the phase diagram is similar to that classifled in 
section llV.AI 

In this section, we discuss the regime when quantum fluctuations are non-negligible and focus on the global phase 
diagram when the classical order competes with the quantum disorder and uniform states compete with non-uniform 
states. In Fig. and table II, we see that the classically ordered states are obtained from the linear superpositions 
of quantum disordered states. The quantum disordered states are realized in the regime where the kinetic energy 
of the superspin As and Ac overwhelms the coupling energy of the superspin Js and Jc, and the superspin vector 
becomes disordered in the temporal domain. In this sense, the quantum description of the superspin goes far beyond 
the classical LG theory discussed in the previous section. 

By arranging the six elementary states from Fig. I^into a spatially non-uniform patterns, we have inflnitely many 
possibilities. In addition to the classically ordered AF and SC states, in Fig. ^] we illustrate some of the basic 
non-uniform states and their associated wave functions, expressed in terms of 0(x), mgjx) and A(cc ) . Stripe order was 
theoretically predicted and exper i mentally observed in the HTSC cuprates ljKivelson et alV Il998t iTranauada et alV 
Il995t IWhite and Scalapinol Il998t IZaanen and GunnarssonL Il989|) . In a typical stripe phase, a magnetic stripe of 
width 2a is separated by a charge stripe of width 2a, where a is the lattice constant. The stripe state come in two 
forms. For the in-phase stripes, both the charge and the spin periodicity is 4a in the direction transverse to the stripe 
direction. For the out-of-phase stripes, the charge periodicity is 4a, while the spin periodicity is 8a. The charge stripe 
can either be insulating, or superconducting. The SC stripes are defined by their phase angle; the two nearby SC 
stripes can be either in-phase or out-of-phase. The case when both the AF and the stripes are out of phase can be 
viewed as a superspin spiral^ in which the superspin direction rotates continuously along the direction transverse to 
the stripes. (See Fig. Et)- Both types of str ipes discussed her e have both AF and SC orders. Another possibihty 
is the checkerboard pair-density-wave rPDWl llChen et QiJ . l2002D . depicted in Fig. It can be obtained from the 

in-phase stripe by quantum disordering the hole pairs in the SC stripe. This state is insulating with AF and charge 
orders. We stress that all insulating states in the 50(5) theory are obtained from the quantum disordered states of 
the hole pairs; therefore, they are paired insulators, in contrast to ordinary band insulators or a Wigner crystal state 
of the electrons. 

Therefore, some of the inhomogeneous states observed in the HTSC cuprates can be described naturally in terms 
of the temporal and spatial ordering of the superspin. The key question is how they are energetically stabilized in 
the projected 50(5) model. These spatially non-uniform states are usually realized when extended interactions are 
considered. These extended interactions take the form 

)nh{x)nh{x') 

{xx') {{xx')) 

+ (F.(5t) E E ) E itH^)tH^'))sAtix)ti^'))s^ 

{xx') {{xx')) Sr=0,l,2 

E (4(a^)ia(a;')4(^')i/«(a^) + H.c.) + V^Y1 iMx)nt{x') + nh{x')nt{x)) + ... (65) 

{xx') {xx') 

Here {xx') and {{xx')) denote the summation over the nearest-neighbor and the next-nearest-neighbor on a square 
lattice. {t{x)t{x'))sT refers to the total spin St = 0, 1, 2 combinations of two magnons on sites {xx'). The Vc and Vc 
terms describe the interaction of the hole pairs, the Vg and V^ terms describe the interaction of the magnons, and 
the Jtt and terms describe the mutual interaction of the hole pair and the magnon. Since the projected 50(5) 
model is defined on a coarse-grained lattice, the density of the hole pairs, n^, is related to the hole do ping density 
by Uh = 2(5. The model Hamiltonian given hy H + Hcxt has been studied extensively by Chen et al ljChen et all 
l2nn3aD by using both quantum Monte Carlo methods and mean field theory. Here we summarize the basic qualitative 
results. In order to study the phase diagram of this model, we first focus on the charge sector. The charge sector 
of the projected 50(5) model is the same as the hard-core boson model introduced in Eq. (TJ of sectional This 
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model has been much studied in the context of superfluid-to-insulator transition(|Bern ardct et al'. '2002| iBruder et al\ . 
Il993t iHebert et all l2002t iPich and Frevl Il998t IScalettar etdl Il995t IvanOtterlo ~aL .1995.) . Without the extended 
interactions V{XtX') in Eq. (PJ, the phase diagram of the hard-core boson is given in Fig. Half-fiUing of the 
original electron systems in the cuprates corresponds to the vacuum state of the hole-pairs, or phase III in Fig. 
The chemical potential induces a transition into the SC state, labelled as phase II. Further increase of the chemical 
potential induces a transition into a checker-board ordered state, labelled as phase I. This is the "class B" transition 
shown in Fig. Phase I corresponds to nt = 1/2 of the hole pair bosons, or (5 = 1/4 of the original electrons. 
When extended interactions in V{x, x') are included, a new insulating phase develops near the overlapping region of 
phase I and phase III, with boso n density of tz^. = 1 /4 of the hole pair bosons, or 5 — 1/8 of the original electrons 
(see, for instance. Fig. 2 of Ref. ijBruder et a/.Lll993|) '). This insulating phase can have either stripe or checkerboard 
like charge order. Generally, str ipe type of insulating order is favor e d for y/ » V^, and the checkerboard-type order 
is favored in the opposite limit ljHebert et aTl 120021: iPich and Frevl Il998() . With even more extended interactions, 
additional phases develop at lower rational densities. These Mott insulating phases at various rational densities are 
shown in Fig. El The phase boundary between the insulating phases with charge order and the SC phases can be 
generally classified into "type 1, 1.5 and 2," according to the terminology developed in section HV. Al and Fig. IHl In 
the last two cases, a mixed phase, called the supersolid phase, develops near the phase boundary. After understanding 
the generic phase diagram of the hard-core lattice boson model, we are now in a position to discuss the full global 
phase diagram of the 5*0(5) model H + H^xt, depicted in Fig. E| Here J/V denote the typical ratio of Jc/Vc, but it 
can obviously be replaced by other similar parameters. The — phase corresponds to the AF state at half- filling, 
where magnons condense into the singlet ground state. For large values of J/V, a pure SC state is obtained where 
the hole pairs condense into the singlet ground state. However, besides these two robust, classically ordered phases, 
we also see new insulating phases at rih = 1/4, = 1/8 and Uh = 3/8, which correspond to S = 1/8, 6 = 1/16 
and S = 3/16 in the real system. These new insulating states are stabilized by the extended interactions and have 
both AF and PDW order (see example Fig. I12tl'l. As the chemical potential or the doping level is varied, a given 
system traces out a one dimensional slice in this phase diagram, with typical slices Bl, B2 and B3 depicted in Fig. 
I13f we expect the quantum parameter J/V to be independent of fi for a given family of materials). The nature of the 
phase transition Bl is similar to that of the classical model already discussed in section HV.AI In this case, the phase 
transition from the AF to SC state can be further classified into "types 1, 1.5 and 2," as discussed in section llV.AI 
with the two latter cases leading to an AF/SC mixed phase at the phase transition boundary. For lower values of 
J/V, the trace B3 encounters the 5 = 1/8 insulating phase. The key signature of this type of phase transition is that 
the SC Tc will display a pronounced minimum as the doping variation traces through the 6 — 1/8 insulating state. At 
the same time, the AF ordering (possibly at a wave vector shifted from (tt, tt)) will show reentrant behavior as doping 
is varied. The phase transition around the fractional insulating phases can again be classified into types "1, 1.5 and 
2," with possible AF/SC, AF/PDW, SC/PDW and AF/PDW/SC mix phases. 

So far we have classified all quantum phase transitions in the 5*0(5) models according to two broad classes. "Class 
A" describes transitions at a fixed chemical potential, typically at an effectively particle-hole symmetric point. "Class 
B" describes transitions in which the chemical potential or the density is varied. Each broad class is further classified 
into three "types, 1, 2 and 1.5," depending on whether the transition is a direct first order, two second order, or an 
intermediate symmetric point in between. The full quantum 5*0(5) symmetry can only be realized in the "class A 
type 1.5" quantum phase transition. The Heisenberg point in the hard-core boson problem discussed in section ITU is 
one such example. The g — point in the 50(5) quantum non-linear a model H34f) is another example. On the other 
hand, the static, or projected 50(5) symmetry can be realized in "class B type 1.5" transitions. We believe that the 
AF to SC transitions in the YBCO, BCCO and NCCO systems correspond to "class Bl" transition. These systems 
only have an AF to SC transition, which can be further classified as types "1, 1.5 and 2," but they do not encounter 
additional statically ordered fractional insulating phases. On the other hand, the phase transition in the LSCO system, 
where Tc displays a pronounced dip at 5 = 1/8, corresponds to the "class B3" transition (see Fig. I13f) . In the HTSC 
cuprates, the charge gap at half-filling is very large, of the order oiU ^ 6eV; it is not possible to induce the "class Al" 
transition from the AF to SC state by conventional means. However, the charge gap at the fractional insulating states 
is much smaller, of t he order of J, and it is possible to induce the "class A2" insulator to superconductor transition 
by applying pressure ijArumueam et aO . l2002t [Locauet et al ]. ll998HSato et ad 120001 iTakeshita eroll l2003j) . It would 
be interesting to determine if this transition point could have the full quantum 5*0(5) symmetry. 

Therefore, we see that the concept of the 50(5) superspin indeed gives a simple and unified organizational principle 
to understand the rich phase diagram of the cuprates and other related systems. This construction of the global phase 
diagram can obviously be iterated ad infinitum to give a beautiful fractal structure of self-similar phases and phase 
transitions. All of this complexity can be simply reduced to the five elementary quantum states of the superspin. 
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C. Numerical simulations of the classical and quantum models 

In this section, we revi ew essentially exact numerical studies of the classical S O(5) model and the quantum pro jected . 
50(5) model on a lat ticefC hen et aL 2003a';'Dorneic h et a^J .l2002b':'H?. ll999lj. 120011: iJostingmeier et a;J . l2003tlRieraL 
I22029 b). In section IV.DI we shall discuss the transformation from the microscopic models into the effective 5*0(5) 
models and determine the effective parameters. Once this is accomplished, the phase diagram of the model can be 
determined reliably by bosonic QMC simulations. These calculations can be carried out for system sizes up to two 
orders of magnitude larger than fermionic QMC simulations, wit h the latter being p lagued, in the physically interesting 
regime - i. e., close to half- filling - by the minus-sign problem (|der Lindenlll993) . The effective models can also be 
studied numerically in 3D; this is crucial, since there exists no AF ordered phase in 2D at finite temperature (nor 
long-range SC order). Thus, we are forced to study the 3D case in order to determine the phase diagram and to show 
that the scaling behavior is consistent with an 50(5) symmetric critical behavior within the parameter regime studied 
(temp erature and system size). This was possibl e due to a major step forward in the numerically accessible system 
sizes ijDorneich and Troveil I2OOII: ISandvikL Il997l Il999j) : in the boso nic projected 50(5) model ^10.000 sites were 
included, in contrast to just ~100 sites in fermionic QMC calculations fDagotto', '1994'; "Do pf et all 11992; ' Imada et all 
^98). The numerical results, obtained by the QMC technique of Stochastic Series Expansion (SSE) (Sandvik, 19971 
[1299) and reviewed here, show that the projected 50(5) model can give a realistic description of the global phase 
diagram of the HTSC cuprates and accounts for many of their physical properties. 

The form of the projected 50(5) Hamiltonian is given in Eq. (|50|l . The extended 50(5) model also include s the 
interactions expressed in Eq. We shall discuss the simple 50(5) model first. In Ref. l)Zhang et a/.lll999[l . this 

Hamiltonian was studied analytically within a mean-field approach. At the special point Jj, = 2Js = J and Ag — A^, 
the mean-field level of the ground-state energy of Hamiltonian ( 1501) depends on the AF and SC order parameters 
X ~ (tl.) and y = (fj^) only via their combination + y^, which reflects the 50(5) invariance of the mean-field 
approximation. In the full model, however, quantum fluctuations modify the zero-point energy of the bosons in 
Eq. f I50|l : therefo re, giving a co rrection to the ground-state energy which depends on x and y separately and destroys 
50(5) symmetry l)Zhang et al\ . \l99^ . Hence it is essential to study the full quantum-mechanical model (|5n|l . including 
all quantum fluctuations, which can only be done by means of numerical simulations. We then compare the properties 
of the projected 50(5) model first in two dimensions (2D) with a variety of salient features of the HTSC such as 
the global phase diagram and the neutron-scattering resonance. Finally, we review an extension of these studies to 
the 3D projected 50(5) model. In particular, we show that the scaling behavior near the multi-critical point, within 
the parameter regime studied (system size and temperature), is consistent with an 50(5) symmetrical behavior. The 
departure away from 50(5) symmetric scaling can only occur in a narrow parameter regime which is hardly accessible 
either experimentally or numerically. 

After numerically solving the projected 50(5) model, we obtain Fig. 114b . which gives the mean hole-pair 
and magnon densities as a function of the chemical potential for T/J = 0.03 and their T extrapolations 
Ijostingmcicr et aZ., 2003). Similar to the mean-field results, ajump in the densities can be clearly seen at ^c = — 0.175, 
with a shift in respect to the mea n-field value due to the stronger fluctuations of hole pairs, as seen in the Gaussian 
contributions ijZhang et all Il999() . The nature of the phase transition a,t fi — can be determined by studying 
histograms of the hole-pair distribution for fixed fj, = ^c- While in a homogeneous phase the density peaks at its mean 
value, at ^i = fj,c we obtain two peaks, which indicates a first-order transition with a phase separation between (almost) 
hole-free regions and regions with high hole-pair density. From Fig. 114b we see that the transition is of first order for 
T<rp = (0.20± 0.01)J at /i = /ip = (-0.168 ±0.002) J. Above Tp, the histograms show strongly fluctuating hole-pair 
densities, suggesting the presence of critical behavior. 

Based on these results, the phase diagram of the 2D projected 50(5) model is obtained in Fig. ^| Unlike the 
generic three dimensional phase diagrams presented in Fig. 1101 there can be no finite temperature Neel transition in 
D=2 because of the Mermin- Wagner theorem. On the other hand, a continuous transition of the Kosterlitz-Thoulcss 
(KT) type is possible for the SC to normal state transition at finite temperature. The 3D phase diagram shown in 
Fig. llOb takes the form of Fig. ^|in D=2, where the first order line separating the AF and SC phases merges into 
the continuous KT transition at a tricritical point P. The SC phase with finite superfluid density is identified by a 
power-law decay of the SC correlation function: 

C/.(r)-(4(r)+t,,(r))(4(0)+<;,(0)) . 

The KT transition line in Fig. I15h , separates power-law {Ch(r) cx r~") from rapid exponential decay {Ch(r) (x e"^*"). 
A reliable and accurate distinction between these two decay behaviors requires a finite-size scaling with large system 
sizes, as well as an efficient QMC estimator for the Green functions appearing in the correlation function. With its 
non-local update scheme and with our new estimators for arbitrary Green functions, SSE provides both (for details 
see Ref. CDorneicli and Trover,.2001^) ). An alternative method for detecting a KT transition exploits the fact that the 
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superfluid density jumps from zero to a finite value at the KT temperature Tr-t (* Nelson and Kos terlit^.'lQTT*). Within 
SSE the superfluid density can be measured quite easily by counting winding numbers l)Harada~nd KawashimaL .1997). 
Numerically, this criterion is preferable to the arduous process of direct determination of decay coefficients. Fig. 115b . 
plots the phase diagram obtained by applying both criteria independently. The figure shows that the projected 5*0(5) 
model indeed has a KT phase with quasi long-range order whose dome-like form in /i-T space looks like that of the 
HTSC cuprates. Both criteria produce the same cle arly pronounced phase se paration line. It is well known that a 
similar transition cannot occur for antiferromagnets ljChakravartv et all . 1198 8^ and that the finite-T AF correlation 
length ^ is always finite and behaves like ^ oc q'^^pJ^bT ^ with ps being the spin stiffness. This fact is confirmed by 
our numerical results. 

One condition required for an SO{b) symmetric point is that the formation energies of hole-pair bosons and of 
magnons are identical. This condition is fulfilled along the line from S to the tricritical point P in Fig. ^| Another 
necessary condition is that hole pairs and magnons behave in the same way at long distances. This condition is fulfilled 
on the dashed line in Fig. where the AF and SC correlation lengths ^ become equal. Interestingly, these two 
conditions are met (within error bar accuracy) at the tricritical point P. Of course, the correlation length is still finite 
here; however, we find relatively large ^ values of order 10 to 15 in the immediate vicinity of point P, demonstrating 
the importance of 5'0(5) critical fluctuations in this region. 

In addition, in realistic electron systems, the long-range part of the Coulomb repulsion between the doubly-charged 
hole pairs disfavors phase separation, while extended short ranged interactions described by Eq. (|65|l could lead to 
the formation of stripes and checkerboard types of states, as discussed in section HV.BI To study the effect of off-site 
Coulomb interaction, we have added additional nearest-neighbor and next-nearest-neighbor Coulomb repulsions Vc 
and V^ = 0.67Vc to the projected 5*0(5) model. Indeed, a relatively modest Coulomb repulsion of T4/J«0.2 is enough 
to completely destroy the phase separation. Thus, one interesting effect of Coulomb interaction in two dimensions is 
to push down the tricritical point into a quantum-critical point at r = 0. In section FlV. Al and in Fig. we showed 
that the 5*0(5) symmetric behavior is recovered at the special point when a direct first order transition changes into 
two second order transitions. Therefore, the extended Coulomb interaction plays the role of the w parameter in Fig. 

and could restore the 50(5) symmetry at the quantum critical point. 

When larger values of extended interaction parameters in Eq. Ht)5l) are considered, new in sulating phases are 
expected, following from the general discussions in scction llV. Bl and Fig. E| Indeed, Chen et al l)Chen et a/.l . l2003ai) 
have performed extensive QMC simulation of the 50(5) model and have determined its generic phase diagram, as 
shown in Fig. 1161 In addition to the AF and SC phases, there is an insulating pair-density-wave state around doping 
range of x = 1/8, where hole pairs form a checkerboard state in the AF ordered background, as depicted in Fig. 
I12H . Near the phase boundaries between the AF, PDW and SC phases, there are mixed phases with coexisting order. 
The topology of the phase diagram obtained from the QMC simulation agrees well with the mean field theory of the 
extended 50(5) model. One of the main features of the 50(5) theory is that it provides an e l egant explanation fo r the 
neutron resonance peak observed in some HTSC cuprates at a = (tt. Trl llDemler and Zha^ Il995t IZha^L Il99it . We 
refer the reader to the detailed discussion of the resonance mode in section fvTI Experiments show that the resonance 
energy Wrcs is an increasing function of Te, i.e . LOies increases as a function of doping in the underdoped region and 
decreases in the overdoped region lFong et aIl . l200Cl(l . Fig. 115b plots the resonance frequency determined from the 
spin correlation spectrum obtained for the projected 50(5) model. As illustrated in Fig. |S1 the spin- wave excitations 
are massless Goldstone modes in the AF phase at fj, < (and T = 0) and become massive when entering into the 
SC phase. Wrcs increases monotonically up to the optimal doping ^opt ~ 1- In the overdoped range of the simple 
5*0(5) model, however, uj^cs is increasing more, in contrast to what happens in the cuprat es. The resonance peak 
continuously loses weight as /i increases, which is consistent with experimental observations ijFong et all . \200dl} . 

A comparison of the critical temperature Tc obtained from Fig. 1151 and the resonance frequency Wres at optimal 
doping yields the ratio Tc/wresjopt = 0.23. Again, this is in qualitative accordance with the corresponding ratio for 
YBa2Cu306-i-x, for which the experimentally determined values Tc = 93 K (thus fcs?c = 8.02 meV) and Wrcs.opt = 41 
meV yield Tc/wi.cs,opt = 0.20. 

Now we turn to the numerical simulations of the 5*0(5) models in D=3. Two aims motivate our studies of the 
projected 50(5) model in three dimensions (3D). First, we expect to find an AF and SC phase with real long-range 
order. We need to determine which of the two types of phase diagrams introduced in section IIV.AI (see Fig. I1U|) 
is realized in the numerical simulations. Second, we would like to determine whether the projected 50(5) model 
has a certain multi-critical point at which the 50(5) symmetry is asymptotically restored. Since the cuprates have 
a pronounced 2D layer structure with relatively weak couplings between adjacent Cu02 planes, the 2D and the 
isotropic 3D model (discussed here) should be two extreme poles for the possible range of properties of real HTSC 
material s. Most numer i cal data reviewed here have been obtained by a finite- s ize sc aling with lattice sites up to 
10,000. l|Dorneich et a^J . l2002albl: iDorneich and Troveill200ltlJostingmeier et al\ . \200^ 

The pha se diagram and the scaling behavior of the classical 5*0(5) model has been studied in detail by Hu in Ref. 
|h1 . I2001D by means of classical Monte Carlo simulations (MC). Classical MC are by orders of magnitude easier to 
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perform and less resource demanding than QMC simulations; hence, very large system sizes can be simulated and 
highly accurate data can be obtained. The classical S0{5) model can be obtained directly from the quantum 5*0(5) 
model by taking the expectation value of the Hamiltonian in the variational state, as given by Eq. (j64|l and assuming 
a constant value of 9{x). It takes the form: 

= -J ^ ma{x)ma{x') + g'^niKx) +w'^ml{x)mf{x), (66) 

{x.x') X X 

where g = — is the quadratic symmetry breaking term, and w is an additional quartic symmetry breaking term. 
Hu established the T{g) phase diagram, which is of the type illustrated in Fig. llUb . The model has an AF and SC 
phase which meet at a bicritical point {Tf^c, gbc = ^)- The boundary lines between the disordered and AF phases, and 
between the disorder ed and SC phases merge tangentially at the bicritical point, which is an important characteristics 
of 5*0(5) symmetry ijHuL l200ll) . The following scahng properties were determined by Hu and will be used to study 
the restoration of 5*0(5) symmetry in the projected 50(5) model. 

For an analysis of the crossover phenomenon, an Ansatz for the behavi or of the helicity modulus T in the range 
T<Tc{g) and g>0 \s used, which is suggested by scahng theorv l|HiJ . l200lt) : 

T(T,5) « {g-gbcT''^ x f{{T/n, -l)/{g~g,,)^/^) . (67) 

Here, is the critical exponent for correlation length at ri = 5 and 4> the crossover exponent. Using H67I) . the values 
of and </> can be determined in two steps. First, performing a g scan of T(T — The, g) returns the ratio v^/(j): 

T(r,„g)/T(T,„.g') = {{g - g,,)/{g' - g,,))"''"^ . (68) 

Then, (j) is obtained from the slopes -^{T (T, g)/T{T, g')") via 

, , f 92-9bc \ L ( d TiT,g,) I d T(T,ff2) \ 



if gi, g[, 32, and are related by (51 - gtc)/ {g[ - gbc) = [92 - gbc)/{92 ~ 9bc) > 0. From the scahng plots presented 
in Figini Hu finds the values 2^5/0 = 0.523 ± 0.002 and 0= 1.387 ± 0.030. 

According to the scaling Ansatz in H67|) . the transition lines between the disordered and AF phases, and between 
the disordered and SC phases near the bicritical point should be of the form 

B,-ig-g,^r/^^^^-l and B, ■ {g,^ ~ g)^/^ = - 1 . (70) 

J 6c -Lbc 

The ratio B2/B3 should be given by the inverse ratio between the AF and SC degrees of freedom, i.e. 

82/33 = 3/2. (71) 

The values numerically determined by Hu indeed have the correct ratio: ,62 = 1/4 and B :^ = 1/6 . 

We now proceed to the phase diagram of the 3D quantum 50(5) model (Dorncic h et al\ . l2002albl) . Figure 1181 
shows the AF and SC phases, as expected. Furthermore, the two phase transition lines merge tangentially in to a 
multi-critical point (at Tbc = 0.960 ± 0.005 and gtc = —0.098 ± 0.001) just as in the classical 5*0(5) system llHul . 
|5001). The line of equal correlation decay of hole-pairs and triplet bosons also merges into this bicritical point P - 
a necessary con dition at t his point for the restoration of 50(5) symmetry. Unlike the corresponding phase in the 
classical model i HuL l200l|) . the SC phase only extends over a finite g range due to the hardcore constraint of the 
hole-pair bosons and agrees with experimentally determined phase diagrams of the cuprates. Obviously, the quantum 
mechanical 50(5) model is 'more physical' in this aspect than the classical 50(5) model. In real cuprates the ratio 
between the maximum temperatures and T/v is about 0.17 to 0.25, whereas in the projected 5*0(5) model we obtain 
the values Tc/J= 1.465 ± 0.008 at ^lopt/ J ~ 1-7 and T^/ J = 1.29 ± 0.01 at /i ^ co\ hence, Tc is slightly larger than 
Ttv. In order to obtain realistic values for the transition temperatures, it is necessary to include the and T^- terms 
in Eq. I|65|) . These terms represent the repulsion between the magnons and the hole pairs. If we take the expectation 
values of the hole pair operators, these terms effectively represent a doping dependent Js, which can produce a more 
realistic phase diagram. Such terms break the 5*0(5) symmetry of the static potential at T=0 (see Section llH. CI after 
equation (|48|l '). However, the static symmetry may still be recovered at the bicritical point, as discussed in Section 
IIV.AI At this point we are primarily concerned with the multi-critical behavior, so we stay with the simple 50(5) 
model. 
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A closer look at the phase transition line between the points S and P (see Fig. I18|) reveals that this line is slightly 
inclined, unlike the vertical line seen in the classical 5*0(5) model. This indicates that a finite latent heat is connected 
with the AF-SC phase transition. In addition, this means that fi is not a scaling variable for the bicritical point P, as 
it is in the classical model. The result in Fig 1 181 shows a phase separation regime at /i = /ic on the entire transition 
line from S to P. 

We now rev iew the results of a scaling an alysis for the 3D quantu m 5*0(5) model, similar to the one performed by 
Hu IHuI. 120011) in a classical 50(5) system (iDorneich et all l2002aj) . From this analysis we also find that the 50(5) 
symmetry is restored in the region around the bicritical point (Tf,c = 0.96, /i = — 0.098). 

We have determined the critical exponents for the onset of AF and SC orders for various chemical potentials as a 
functi on of temperature. Far into the SC range, at /i = 1.5, we find that the SC helicity modulus follows the scaling 
form ijFisher a/.l Il973|) 

T cx (1 - T/TcY with u = 0.66 ± 0.02 , 

which agrees with the values obtained by both the e-expansion and numerical analysis of a 3D XY model. On the AF 
side, error bars are larger. For fi = —2.25, 

CAFioo) cx (1 - T/T.f'^ with 133 = 0.35 ± 0.03, 

as expected for a 3D classical Heisenberg model. 

To determine v and 0, we use Eqs. (I68|) and (|69|1 . which express the scaling behavior in the crossover regime (cf. 
Ref. ljHl . l200ll) ). We obtain the ratio 

1/5/0 = 0.52 ±0.01, 

which matches the results of the 6-exDansion (|Hu and Zhan3. 1200(1 iKosterlitz et all ^2^. (/> is then obtained by 
using H69|l . The result is 

(/)= 1.43 ±0.05 

which also agrees with the e-expansion for an 50( 5) bicrit ical point and with the results of Ref. l|HuL 12001]) . 

Let us finally return to the comment by Aharony l) Aharonv, ,2002) . who, via a rigorous argument, demonstrated that 
the decoupled fixed point is stable, as opposed to the biconical and 50(5) fixed points. However, he also commented 
that the unstable flow is extremely slow for the 5*0(5) case due to the small crossover exponent. 

The scaling analysis of the 3D projected 50(5) model has produced a crossover exponent which matches the value 
obtained from a classical 50(5) model and from the e-expansion. This provides strong evidence that the static 
correlation functions at the 50(5) multicritical point are controlled by a fully 5*0(5) symmetric point, at least in a 
large transient region. However, the isotropic 50(5) and biconical fixed points have very similar critical exponents. 
Thus, given the statistical and finite-size errors, as well as the errors due to the extrapolation of the e-expansion value 
to e = 1, we cannot exclude the possibility that the multicritical point on the phase diagram is actually the biconical 
one. On the other hand, the biconical fixed point should be accompanied by a uniform AF/SC mixed region (as a 
function of chemical potential), which was not observed. The decoupled fixed point appears to be the least compatible 
with the nu merica l results presented above. Even if the bicritical point is fundamentally unstable, as suggested by 
Aharony in ijAharo nv. 2002^, one would have to come unr ealistically close to Ttc to observe this. For example, for the 
projected 50(5) models Ref. (jArrigoni and Hankelll999() estimated that deviations from the 5*0(5) behavior may be 
observed only when the reduced temperature becomes smaller than \T — Tbd/Tbc < 10"^"'^. On the other hand, the 
other scaling variables, although initially of the order of 1, rapidly scale to zero due to the large, negative exponents. 
Therefore, the 5*0(5) regime starts to become important as soon as the AF and SC correlation lengths become large 
and basically continues to affect the scaling behavior of the system in the whole accessible region. 

Summarizing, the accurate QMC calculations show that the projected 5*0(5) model which combines the idea of 
50(5) symmetry with a realistic treatment of the Hubbard gap, is characterized by an 5*0(5) symmetric bicritical 
point, at least within a large transient region. Possible flow away from this symmetric fix point occurs only within 
a narrow region in reduced temperature, making it impossible to observe either experimentally or numerically. This 
situation is com mon to many system s in co ndensed matter physics. For example, due to the well-known Kohn- 
Luttinger effect l)Kohn and Luttingerl Il965j) . the Fermi-liquid fixed point is always unstable towards a SC state. 
However, this effect is experimentally irrelevant for most metals since it only works at extremely low temperatures. 
Another example is the "ordinary" superconductor to normal-state tr ansition at T^. Strictly speaking, due to the 
coupling to the electromagnetic field this fixed point is always unstable l)HalDerin et aZllT974|) . However, this effect is 
experimentally irrelevant since the associated critical region is extremely small. Similarly, irrespective of the question 
of ultimate stability, the 50(5) fixed point is a robust one in a large transient regime, and it can control the physics 
near the AF and SC transitions. For all practical purposes, the multi-critical point is dominated by the initial flow 
toward the 50(5) symmetric behavior. 
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V. MICROSCOPIC ORIGIN OF THE S0(5) SYMMETRY 
A. Quantum lattice models with exact S0(5) symmetry 

Soon after the general .SO (5) t heory was proposed, a class of microscopic fermion models w i th exact SO(5) sy mmetry 
was constructed (|Burgess et all 119981: iHenlevi Il998l: lR,ahello et d[ Il998l: IScalapino eTdl Il998l: IWii et all H)03a). 
These models fall into three general classes. The first class contains models with two sites per unit cell, such as the 
ladder and the bi-layer models. In these models, a simple condition among the local interaction parameters ensures 
the full quantum 5*0(5) symmetry. The second class contains models with only one site in the unit cell but with longer 
ranged interactions. The third class contains higher spin fermion models, in particular the spin 3/2 Hubbard model. 
Remarkably, in this case the models are always SO{5) symmetric without any fine tuning of the local interaction 
parameters and doping level. 

The microscopic 50(5) symmetri c models in the lad der or bi-layer models were first constructed by 
Scalapino, Zhang and Hanke fSZH'llls'calapino all Il998^ and have been studied extensively both analyti- 
cally and numerically (Arrigoni and Hankc, 1999*: Bouwknegt and Schoutens, 1999; Duffy et al, 1998; Edcr et " 
[1999; Frahm and Stah lsmeier, 2001; Furusaki and Zhang, 1999; Hong and Salk, 1999; L in et ai, 1998; Schulz, 199 



Shelton and Senechaill998lh . In these models, there are two sites and 4^ = 16 states in the unit cell. In section llll.AI 
we already discussed the construction of 5*0(5) symmetry operators in terms of the fermion operators for two sites in 
the unit cell. Here we shall address the question of whether the microscopic Hamiltonian commutes with the 50(5) 
symmetry generators. Three interaction parameters, [/, V and J, fully characterize the most general local interactions 
on the two sites, which takes the form 

H{x) = U{nci - i)(n,i - 1) + (c ^ d) + V{n, - l){nd - 1) + JSja - /i(nc + Ud). (72) 

This Hamiltonian can be solved easily for the 16 states on two sites and the 6 energy levels are given in Fig. 1201 
Since the 50(5) symmetry generators can be expressed in terms of the microscopic fermion operators, we can easily 
determine the transformation properties of these states under the 5*0(5) group. There are three 50(5) singlet states, 
and two fermionic quartet states, which form the fundamental spinor representations of 50(5). We see that the four 
fermionic states in each group are always degenerate, without any fine tuning of the interaction parameters. The 
three spin triplet states at half filling and the two paired states away from half-filling form the five dimensional vector 
representation, but they arc only degenerate if we specify one condition, namely 

J^A(U + V). (73) 

This condition ensures the local 50(5) symmetry within the unit cell. Remarkably, under this condition, a global 
5*0(5) symmetry is also obtained for a bi-partite lattice including nearest-neighbor hopping. This is best demonstrated 
when we write the model in a manifestly 5*0(5) covariant manner. On a bi-partite lattice, we introduce the four- 
component spinor operator 

*„(x e even) = ''j^j^J^ ) ^^{x S odd) = g ) . (74) 
The microscopic Hamiltonian including intra-rung hopping t± and inter-rung hopping t|| is given by 

H = -2t|| (4(a^)c.(^') + dUx)da{x')) - 2t^ ^(4(x)d,(x) + h.c) + ^ H{x). (75) 

{x.x') X X 

Under condition (|73|l . this Hamiltonian can be expressed in a manifestly SO{5) invariant manner: 

H = 2t|| ^ {^^{x)R'^^^^ix') + h.c.) + t^i^^R'^f^^, + h.c.) + J2 {lUx) + (^ + ^)(*1*„ _ 2)^ (76) 

{x.x') X 

where the R matrix is defined in the Appendix. This model was originally constructed for the two-legged ladder 
system, but it works equally well for a two dimensional bi-layer system. 

The phase diagram of this 5*0(5) symmetric model has been studied extensively in the literature. This simple model 
has a rich and rather complex phase diagram, depending on the coupling strength and doping. However, because of 
the constraints imposed by the 5*0(5) symmetry, the phase diagram is much better understood compared to other 
related models. In the strong coupling limit, three phase boundary lines are determined from the level crossing of the 
bosonic states on two sites. At = — 2[/, the Eq state becomes degenerate with the states; aX V = —U, the Eq 
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state becomes degenerate with the Ei states; finally, at V = the Ei states become degenerate with the E^ states. 
The strong coupling phase diagram at half-filling is shown in Fig. 120b . 

In the strong coupling Eq phase, a robust ground state is obtained as a product of 5*0(5) singlets on the rungs. 
This type of insulating state does not break any lattice translational or internal rotational symmetry. Since there are 
two electrons per unit cell, this insulating state is also adiabatically connected to the band insulator state. This state 
is separated from the excited SO{5) quintet vector states by a finite energy gap, A = Ei — Eq = J. In this regime, 
we consider the low energy manifold consisting of six states, namely one Eq state \Q) and five Ei states na\^) per 
rung. The low energy effective Hamiltonian can be obtained easily by the second order strong coupling expansion 
and is exactly given by the 50(5) quantum non- linear a model Hamiltonian given in Eq. (|34|l . with — J ^^'^ 
P = J\\ = tfJ{U + J/2). The operators Lat and act on the six states in the following way: 



Since the quantum model is exactly S'0(5) symmetric, the anisotropy term V{n) vanishes identically. Therefore, 
we see that the 50(5) quantum non-linear a model, phenomenologically introduced in section Fill. Bl can indeed be 
rigorously derived from the microscopic SZH model defined on a ladder and on a bi-layer. 

In the Eq regime, the SZH model on the half-filled ladder has a 50(5) rung singlet ground state with a finite gap 
towards the 50(5) quintet excitations. A chemical potential term of the order of the gap induces a second order 
quantum phase transition into the SC phase. On the other hand, the SZH model on the bi-layer has a quantum phase 
transition even at half-filling, when J\\/J ~ 1. For J > Jy, the ground state is a Mott insulator without any symmetry 
breaking with a finite gap towards the quintet excitations. For J < Jy , the ground state is classically ordered and 
breaks the 50(5) symmetry spontaneously by aligning the superspin in a particular direction, which can be either AF 
or SC. Since the residual symmetry is 50(4), the Goldstone manifold of the a model is the four dimensional sphere 
50(5)/50(4) = 5''*. Away from half-filling, the 50(5) symmetry is broken by the chemical potential term. According 
to Table I, the tt operators carry charge ±2, and we have [H,ttI^] = 2/i7rjj,. However, although the Hamiltonian does 
not commute with all the 50(5) generators, it still commutes with the Casimir operator L^^. For this reason, all states 
are still classified by 50(5) quantum numbers and the 50(5) symmetry makes powerful predictions despite a broken 
symmetry away from half-filling. The phase diagram for the two dimensional SZH bi-layer model is shown in Fig. 
1211 For Jy 3> J, the ground state is classically ordered. The chemical potential induces a quantum phase transition 
from the 50(5) uniform mixed AF/SC state to the SC state at = 0. This transition is exactly the superspin fiop 
transition discussed in section UlI.BI For J|| <C J, the ground state is quantum disordered at half-filling. A second 
order quantum phase transition from the singlet Mott insulator state to the SC state is induced at finite = /ic. 
The exact 50(5) bi-layer model offers an ideal theoretical laboratory to study the collective modes, especially the tt 
resonance mode discussed in section Ivll since their sharpness is protected by the exact 50(5) symmetry. The Mott 
phase has five massive collective modes, a doublet of charge modes and a triplet of spin modes. The energy of the 
two charge modes splits at finite chemical potential, and the energy of one of the charge modes vanishes at the second 
order phase transition boundary. This charge mode continues into the SC phase as the phase Goldstone mode. The 
spin triplet mode of the Mott phase continues smoothly into the SC phase and becomes the pseudo-Goldstone mode, 
or the TT resonance mode of the SC phase. The ordered phase at half-filling has four Goldstone modes. The direction 
of the order parameter can be smoothly rotated from AF to SC at half-filling. When the order parameter points 
in the AF direction, the four Goldstone modes decompose into two spin wave modes and two charge modes. When 
the order parameter is rotated into the SC direction, the four Goldstone modes decompose into a spin triplet and a 
Goldstone phase mode. The energy of the triplet Goldstone mode (the massive tt mode) increases continuously with 
the chemical potential, while the phase Goldstone mode remains gapless. 

Having discussed the Eq regime at length, let us now turn to the Ei regime, where the 50(5) quintet state has 
the lowest energy. In this case, we can restrict ourselves to the low energy manifold of five states on each rung. 
The effective theory within this low energy manifold can again be obtained by the strong coupling second order 
perturbation theory, and is given by 



where K = iy/([//2 — J/4). This effective Hamiltonian is the 50(5) generalization of the AF spin 1 Heisenberg 
model. Here we must distinguish between the one dimensional ladder model and the two dimensional bi-layer model. 
In one dimensional models, the ground state is separated from the 50(5) vector excitation by a finite energy gap. 
In fact, an exact ground state can be constructed for the 50(5) vector model by generalizing the AKLT model for 
the spin 1 chain. Such a state also preserves the lattice translational and internal rotational symmetry. However, in 



Lab{x)\fl{x)) = , Lab{x)\ncix)) = iSbc\naix)) ~ iSac\nb{x)) 
na{x)\n{x)) ^ \naix)) , na{x)\nb{x)) ^ Sab\^{x)). 



(77) 




(78) 



30 



two dimensional bi-layer models, the effective exchange coupling between the 5*0(5) vectors will lead to a state with 
spontaneously broken SO{5) symmetry, with the 50(5) adjoint order parameter (Lab) 0. This order parameter is 
formed by the linear superposition of two S'0(5) vector states, and n;,. Without loss of generality, let us consider 
the case where (iis) ^ 0. In this case, the 50(5) generators L15, {^23,^24,^34} leaves the state invariant. These 
set of generators form a U{V) x SU{2) symmetry group. Therefore, the Goldstone manifold is the coset space 

S0{5)/{U{l)x SU{2))^CP3, (79) 

where OP3 is the six (real) dimensional complex projective space, which can be described by the complex coordinates 
(2:1,22,2:3,^4), satisfying |2ip + I22P + I23P + I24P = 1 and with the points related by a U{1) gauge transformation 
Zi e^°'Zi identified. Since the OP3 manifold is six dimensional, there are six Goldstone bosons in this case. Here 
we see that there is an important difference between the 5*0(5) symmetric SZH model and the 50(3) symmetric 
Heisenberg model. In the Heisenberg model, the vector representation is identical the adjoint representation, there is 
only one type of classically ordered AF state. In the 5*0(5) case, the symmetry breaking can occur either in the vector 
or the adjoint representations of the 50(5) group, which are inequivalent, and the resulting Goldstone manifolds are 
5*^ and OP3, respectively. The adjoint s ymmetry breaking patte rn has been used by Murakami, Nagaosa and Sigrist 



to unify p wave SC with ferromagnetism (|Murakami et a/. . Il999() 



In the weak coupling limit, powerful renormalization group (RG) analysis has been applied to study the 50(5) 
symmetry in ladder models. The main conclusions are similar to the strong cou pling analy s is; th erefore, we will 
only r eview the most remarkabl e and dist inct r esults . Lin, Balents and Fishcr(Lin et ai, 1998), Arrigoni and 
Hanke ljArrigoni and Hank3 . Il999() . Schulz (Schulj . Il998|) . Shelton and Senechal (Shclton jmd S cn echal . -1998) carried 
out detailed RG analysis and showed that RG transformation always scales the most generic ladder model towards 
an 50(5) symmetric ladder model. This is a remarkable result and showed that the quantum 50(5) symmetry does 
not need to be postulated at the microscopic level but co uld emerge as a result of scaling in the long wave length 
and low energy limit. More over, Lin, Balents and Fisher ijLin et all Il998|) showed that even the 5*0(8) symmetry 
could emerge at half-filling. Another interesting and remarkable result was obtained recently. In the transition region 
between the singlet Eq phase and the charge ordered E3 phase, the RG analysis of the weak coupling limit showed the 
existence of a new phase, called the staggered flux phase, or the DDW (d-density-w ave) phase, which has staggered 
circulating current on the plaauettes ljMarston et all l200a ISchollwoeck et al This phase has been prop osed 

to explain the pseudo-gap behavior in the HTSC cuprates i) Affleck and Marstonlll98dlChakravartv et a/ll200lti. 

Ex actly 50 (5) sy mmetric models can also be constructed for the single layer model l)Burgess'^r'ainT998t iHenlevl 
Il998t .Rabello et all [1998) . In this case, there is no natural way to group two sites to form a local, four-component 
5*0(5) spinor. However, one can introduce a 50(5) spinor in momentum space by defining 

= {cpT.Cpi,ff(p)clp+n,T'5(p)c^-p+n,i} ' (^0) 

where g{p) = sgn(cosp2: — cospy) = ±1 is the form factor introduced by Henlev l)Henlevl H 998|) . As discussed in section 
IIII.AI this factor is needed to ensure the closure of the 50(5) algebra. Indeed, with this choice, the ^E* spinors form 
the canonical commutation relation 

V/3} = '5«/3'5p,p'> (81) 

^Ul - {*P«, *p'/3} = -g(p)i?a/3 W,n. (82) 

If we restrict both p and p' to be inside the magnetic Brillouin zone, the right hand side of the second equation 
vanishes and the '^pa spinors commute in the same way as the Cpa spinors. Any Hamiltonian constructed by forming 
singlets of the the basic spinors would be manifestly 50(5) symmetric. 

Because of the non-analyticity associated with the function g{p), this class of 50(5) symmetric models contain 
long ranged interactions in real space. However, similar kinds of long ranged interactions are also present in the 
original BCS model due to the truncation of interactions in momentum space. Therefore, this class of 50(5) models 
can be best viewed as low energy effective models resulting from integrating out states far from the Fermi surface. 
These models may address an important issue in the field of HTSC, which concerns the nature of the quasi-particle 
spectrum at the d wave SC to AF transition. In the pure d wave SC state, the SC order parameter is described 
by the form factor d{p) = {cospx — cospy). When the system is rotated into a uniform mixed AF/SC state, the 
form factor of the resulting AF order parameter is given by g{p)d{p) = \ cospx — cospyj, which contains nodes at the 
same positions as in the pure d wave SC state. When doping is further reduced, a uniform component of the AF 
gap develops across the Fermi surface, filling the d wave nodes. This uniform AF gap gradually evolves in to the AF 
Mott insulating gap at half-filling. See Fig. |^ Based on this scena rio, Zacher et a l (Zacher et al\. l2000() explained 
the d-wave like dispersion of the quasi-particle in the insulating state ljRonning et a^J . il998.) . Filling the d nodes with 
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the uniform AF gap also naturally ex plains the "small gap" observed in the photoemission experiments in the lightly 
doped cuDrates ljShen and et aj E004j) . This theory of the q uasi-particle evolution is a l so similar to the scenario of 
quan tum disordering the nodal quasi-particles developed in l|Balents et all Il998t Il999t iFranz et all l2002hl : iHerbntL 
120021) . Recent studies have found that the generalized Hubbard model for spin 3/2 fer mions enjoys an e xact and 
generic 5*0(5) symmetry without any fine tuning of model parameters and filling factors ijWu et a^J . l2003a|) . Such a 
model can be accurately reali zed in systems of ultra-cold ato r ns on optical lattices, whe re the interaction is local and 
s wave scattering dominates l|Greineill2002HHofstetter et a/J . l2002HJaksch et a/.l . Il998|) . In the Hubbard model with 
spin 1/2 fermions, two fermions on the same site can only form a total spin St = state; the 5t = 1 state is forbidden 
by the Pauli principle. Therefore, only one local interaction parameter specifies the on-site interaction. By a similar 
argument, two spin 3/2 fermions on the same site can only form the total spin St = 0, 2 states; the St = 1,3 states 
are forbidden by the Pauli principle. Therefore, the generalized Hubbard model for spin 3/2 fermions is given by 

H = -tY^ {4aCja + h.C.} ~ 

+ UoJ2P^{{)Po{t) + U2 J2 4nW^2™(*), (83) 
i j,m=±2,±l,0 

where t is the hopping integral, ^ is the chemical potential, and Pq, p\,^-^ are the singlet [St — 0) and quintet {St — 2) 
pairing operators, defined as 

-y/2 2 2 2 2 

-^2 2(*) ~ '^3^1 , P^ = ctc^l, 

' 2 2 ' 2 2 

_ 1 (0 = ct cl 3 , Pl-2i^^cl, cl 3 . (84) 

2 2 2 2 

Remarkably, this generalized Hubbard model for spin 3/2 fermions is always 50(5) symmetric, without any fine tuning 
of parameters and filling factors. This can be seen easily from the energy level diagram of a single site, which contains 
16 states and 6 energy levels for spin 3/2 fermions, as depicted in Fig. ESI The fi'i^.e levels are non-degenerate, the 
degeneracy of the i?2.5 levels is four-fold, and the degeneracy of the level is five-fold. We see that without any fine 
tuning of interaction parameters, this pattern of degeneracy exactly matches the singlet, the quartet (fundamental 
spinor) and the quintet (fundamental vector) representations of the 50(5) group. It can also be easily verified that the 
hopping term also preserves the global 50(5) symmetry. In fact, it preserves an even larger symmetry group, namely 
50(8). The 50(8) symmetry is always broken by interactions; however, under special circumstances, its subgroups, 
50(7), 50(6) and 50(5) x 5t7(2) can be realized in addition to the generic 50(5) symmetry. In this article we 
mainly focus on application of the 50(5) theory to the AF/SC systems. However, from the above discussions, we see 
that ultra-cold atoms on optical lattices also provide a fertile ground for investigating higher symmetries in strongly 
correlated systems, because the higher spins of the atoms and the accuracy of local interaction approximation. In the 
case of the spin 3/2 systems, the generic 50(5) symmetry makes powerful predictions on the symmetries at quantum 
phase transition lines, spectrum degeneracies, topology of the ground state manifolds and low energy effective theories 
of the Goldstone bosons. With the emerging convergence between the atomic and condensed matter physics, we expect 
symmetry concepts and its multiple manifestations to play an ever increasing role in these fields. 

Ferm ions in exact 50(5) models have a beautiful non-abelian holonomy associated with them ljPemler and Zhand . 
11999a'). The four components of an 50(5) spinor represent four states but only two energy levels, each being doubly 
degenerate. As one varies some adiabatic parameters and returns to the same starting value, the states inside a doublet 
can be rotated into each other by a unitary transformation. This interesting mathematical pro perty has been used 
to pre dict 50(5) generalization of the Andreev effect and the non-abelian Aharonov-Bohm effect (jPemler and Zhand . 
ll999aD . 

B. Variational wave functions 

In this section we shall discuss a crucial test of the 50(5) symmetry by investigating the microscopic wave functions 
of the t — J model. In section llV.AI we showed that the transition from the AF state at half-filling to a pure d-wave 
SC state away from half-filling can generally be classified into three types. Within the general form of the static 
potential as given in Eq. H59|) . the "type 1" first order transition is realized for > miM2- For < uiU2, the "type 
2" transition involves two second order transitions with an intermediate mixed phase where the AF and the d-wave 
SC order coexist uniformly. Only for = ""1^2 is an intermediate "type 1.5" transition realized, where the potential 
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can be re-scaled to take an SO{5) symmetric form and a smooth rotation between the AF and the d-wave SC states 
is possible. If we only investigate states with uniform densities, these three possibilities can be distinguished easily by 
curvature in the plot of the ground state energy as a function of doping 5. The curvature would be negative (concave), 
positive (convex), or zero (flat) for these three possibilities, as shown in Fig. @. In the concave case, the uniform 
phase would be thermodynamically unstable, and a Maxwell's construction leads to a phase separated ground state, 
where each phase has a distinct density. 

This interesting prediction can be tested numerically in the t — J model. At this moment, exact diagonalization of 
the t— J model with large system size is not possible due to the exponential growth of the Hilbert space, and reliable 
Monte Carlo simulation cannot be carried out due to the fermion minus sig n problem. A successful method employs 
the va riational Qua ntum- Monte-Carlo (VMC) method (see, in particular, ijCalandra and Sorellal l2000t ICrosl Il989t 
iHimed a and Qgatal Il999() and references therein). Historicall y, the VMC met hod was first used to investigate the 
RVB type of variational wave functions proposed by Anderson ljAndersonL Il987j) . By investigating various variational 
wave functions, this method can address the issue of d-wave pairing in the ground state and the possibility of a uniform 
mixed phase with AF and d-wave SC order for the 2D t — J model. 

The first question is whether the uniform mixe d state h as a lower energy than the pure d-wa ye SC or AF state near 
half-fi lling. In earlier work by Zhang et al. ( Zhang et a^J . ll988,'l and by Yokoyama and Ogata ijYokovama and Qgatal 
ll99fiD . it was shown that the Gutzwiller approximation (GA) gives a reliable estimate for the variational energies for 
the pure d-wave SC state. However, if the AF order parameter is taken into account in the GA, there exists no region 
in the phase diagram where the AF state is stabilized. On the other hand, in Ref. (Himeda and Ogata. 1999) . it was 
shown that when the variational parameters A^, A^^^ and /i were determined from a VMC simulation, where the 
double occupancy prohibition is rigorously treated, then the Gutzwiller-projected trial wave function of the uniform 
mixed state has a lower energy than the pure d-wave SC state with Aaf = 0, in the doping ran ge < 6 < 10%. Using 
Gree ns-function Monte Carlo with stochastic reconfiguration (GFMCSR), Calandra and Sorella l|Calandra and Sorellal 
|200(]|) also concluded that the AF correlations coexists with SC and persists up to (5 = 10%. Himeda and Ogata used 
the following Gutzwiller projected trial wave function: 

I^A) =PG|V'o(Ad,AAF,/^)) , (85) 

where Ad, Aaf and ^ are the variational parameters relating to d-wave SC and AF order and /i is the chemical 
potential. Pg = Hi (1 ~ ^'^iT'^il) stands for the Gutzwiller projection operator. The wave function jV'o (Ad, Aaf, m)) 
is a mixed BCS/Spin-Density-Wave function, i. e. 

|Vo(Ad,Aaf,^)) = n {^k^+^k^difd%)\0) , (86) 

where the index s = {±} takes care of the electron operators acting on the A{B) sublattice i n the AF state. The 
Mfc's and Vk^s contain the variational parameters A^, Aaf and fj. and are defined in detail in ijHimeda and Ogatal 
Il999j) . Fig. 1^ is reproduced from this paper and plots the ground-state energy and the staggered magnetization as a 
function of doping d. 

We see that in the Himeda and Ogata variational QMC work the uniform mixed phase of AF and d-wave SC has a 
lower energy than the pure d-wave SC state up to a doping of about 10%. At half-filling, the energy was found to be 
close to the best estimated value in the Green's function MC method (-0.1994 to -0.20076), which provides support 
for the wave-function Ansatz Eq. H86|) . 

The second point of interest is that, according to the Himada and Ogata results in Fig.[^ the ground-state energy 
is a linear function of doping S in this region, with essentially zero curvature. This implies that the chemical potential 
fi is constant. Since the wave function of Himada and Ogata describes a mixed state with uniform density, the energy 
versus doping plot can generally have three distinct possibilities, as enumerated in Fig. @. Therefore, from the fact 
that the curvature is nearly flat we determine that the condition Mir, = u -\ is fulfilled, which places th e t— J model at 
J/t = 0.3 into the domain of attraction of the 50(5) fixed point ifArrigoni and Hankg . 120001 lMurakami~ and Nagaosal 
1200(1 . 

C. Exact diagonalization of the t-J and the Hubbard model 

In the previous section we discussed the test of the 5*0(5) symmetry through the variational wave functions in the 
t — J model. In this section, we shall describe numerical calculations of the dynamic correlation functions and the 
exact diagonalization of the spectrum, which also tests the 50(5) symmetry of the microscopic t — J and Hubbard 
models. A microscopic model has a symmetry if its generators G commute with the Hamiltonian H, i.e. [H, G] = 0. 
In the 50(5) theory, the tTq operators are the non-trivial generators of the symmetry. In models constructed in 
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section lV.AI the tTq operators indeed commute with the Hamiltonian. However, there are models where the symmetry 
generators do not commute with the Hamiltonian, but they satisfy a weaker condition, [H, G^] = ztAG^, where A is a 
c number eigenvalue (see e.g. Eq. H13() in Section^. These operators are called eigen-operators of the Hamiltonian. 
In this case, from one eigenstate of the Hamiltonian, one can still generate a multiplet of eigenstates by the repeated 
actions of G*. However, the eigenstates within a multiplet are not degenerate, but their energies are equally spaced 
by A. A classic example is the precession of a spin in a magnetic field, where 

H^uoS,; [H,S±]^±u;oS±. (87) 

and llIq is the Lamor frequency of the spin precession. Although in this case the spin-rotational symmetry is broken 
explicitly by the magnetic field in the z-direction and the eigenstates within the multiplets are no longer degenerate, 
the multiplet structure of the symmetry is still visible in the spectrum and can be sampled by the ladder operators. 
If one calculates the dynamical response function of S±, only a single (5-peak is present at w = wq. 

The tTq operators defined in equation (|28|l do not commute with the Hubbard or t — J model Hamiltonian, but 
analytical and numerical calculations show that they are approximate eigen-operators of these model, in the sense 
that 

[iJ,7rt]«c..^t (88) 

is satisfied in the low energy sector. This means that the dynamic auto-correlation function of the tTq operators 
contains a sharp pole at Wtt, with broad spectral weight possibly istributed at higher energies. Using a T- matrix 
approximation, Demler and Zhang ijDemler and Zhang. 1995) verified this approximate equation with = J(l 
n)/2- — 
in a microscopic : 

exact diagonalization technique. Analysis presented in this paper showed that the dynamical correlation function of 
the TT-operator 



'j^i.Ai.J.ClJ iJi-WAJ. ^ VjAJ.J.i.V_(i Clji-AVJ. /_JAJ.C(rXJ.^ £J J.i.(.(jJ.i.(i, , _L I V AAAV^VJ. UiAJ-O CAr^^i VJj^AAAAClJ iJ VjVJ^ LlCAr Li A W AA VV A iJ AA l*^ t/ 1 -A 

-2/x. This calculation will be reviewed in section lVl.DI The first numerical tes t for a low-energy 5* 0(5) syrnmetry 
ricroscopic model has been performed by Meixner et al. jMeixner et aZ.lll997j) using the Lanczos jLanczo£lll95Cl|) 



niicj) = -llm«k„ „J^N+2^. (89) 

(with H being the standard Hubbard Hamiltonian, | ^E*^ ) its ground state with N electrons and the corresponding 
ground state energy) yielded a single sharp excitation peak at low energy w^, accompanied by an incoherent background 
at higher energies. The large separation between the peak and the continuum and the large relative spectral weight of 
the peak demonstrated that indeed the 7r-operator is an approximate eigenoperator of the Hamiltonian (see Fig. I25|l . 
Also in accordance with the perturbative result of Ref. l)Demler and Zhangl Il995(l . the "precession frequency" 
decreases for decreased doping. Furthermore, a comparison with a bubble approximation for this correlation function 
showed that the sharp peak near originated solely from vertex corrections (i.e. collective behavior). 

Not only can the dynamic correlation function of the tTq operators (|89|l be measured numerically for microscopic 
models, thus providing a test of the approximate 5*0(5) symmetry, but they can also be directly measured in neutron 
scattering experiments in the SC state. We shall discuss these experiments in section IvTl 

Exact numerical diagonalization of the t — J and Hubbard models gives eigenstates and eigenvalues on a finite size 
cluster, whose degeneracy pattern can be used directly to test the 50(5) symmetry. In order to explain the main 
idea, let us first examine the variational wave function of the projected 50(5) model given in Eq. fK5|l . This wave 
function describes a broken symmetry state formed by a linear superposition of states with different spin or charge 
quantum numbers. This type of state can only be realized in infinite systems. On a finite size system, all eigenstates 
must have definite spin and charge quantum numbers. Denoting t\x) = ma{x)tl^{x) + A{x)tj^{x) , we can expand the 
coherent state described by Eq. H55|l as 

I *) = {cos 61^ + cos e^-i sin 6* ^ {x) + cos O'^-^ sin 9^ ^ (a;)t^ (y) + cos 9'^-^ sin 6^ ^ {x)t'' {y)t'' {z) + ...}\Q). 

(90) 

For A(x) = 0, we see that the AF ordered state can be expressed as a linear superposition of states with different 
numbers of magnons, forming states with different total spins. While states with different total spins are separated 
by finite energy gaps in a finite size system, these energy gaps could vanish in the thermodynamic limit, allowing 
magnons to "condense" into the ground state. For ma{x) = 0, we see that the SC state can be expressed as a linear 
superposition of states with different numbers of hole pairs, with different total charge. A smooth rotation from the 
AF state to the SC state becomes possible if one can freely substitute each magnon by a hole pair without energy 
cost. This places a powerful requirement on the spectrum. The '^^t^ {x)\^) term in (|90|l contains a single magnon 
state with (5 = 1, Q = 0) or a single hole pair state with {S — 0, Q = —1). 50(5) symmetry requires them to be 



34 



degenerate. This can be easily achieved by tuning the chemical potential, which changes the energy of the hole pair 
state without changing the energy of the magnon state. Once the chemical potential is fixed, there are no additional 
tuning parameters. The J2x^y'^H^)tHy)\^) term in contains a two-magnon state with (S = 2, Q — 0), a one- 
magnon-one-hole-pair state with {S = 1, Q = —1) and a two-hole-pair state with {S — 0, Q = —2). 50(5) symmetry 
again requires them to be degenerate, which is a highly non-trivial test. We can easily perform this analysis for states 
with different numbers of magnons and hole pairs. 

This pattern of the e nergy levels has b een tested directly in the exact diagonalization of the t — J model by 
Eder, Hanke and Zhang ijEder et J . 1 1 998^1 . The t—J model, because of its more limited Hilbert space (no double 
occupancies), allows the exact diagonalization of larger systems (18, 20 and more sites). Since the J model explicitly 
projects out the states in the upper Hubbard band, some of the questions about the compatibility between the Mott- 
Hubbard gap and SO{b) symmetry can also be answered explicitly. In the exact diagonalization studies, total energy, 
momentum, angular momentum, spin and the charge quantum numbers of the low energy states can be determined 
explicitly. These quantum numbers are summarized in Fig. H2t)|) a. 

Eigenstates obtained from the t — J ov Hubbard Hamiltonian can always be interpreted as multi-particles states 
of the underlying electron. However, it would be highly non-trivial if the low energy states could also be interpreted 
as multi-particle states formed from the collectiv e degrees of freedom such as the magnons and the hole pairs. The 
first non-trivial finding of Rcf. (jEder et all Il998|) is that t his is indeed the c ase. Fig. 1271 shows the first four {v = Q 
io V — sets of low-lying states of an 18-site t—J model l)Eder et a?J . ll998|) . We see that the lowest energy state in 
the S = 1, Q = sector indeed has s wave like rotational symmetry and total momentum (tt, tt), as expected from a 
magnon; the lowest energy state in the S — Q,Q — ~1 sector indeed has a d-wave like rotational symmetry and total 
momentum 0. Similarly, states with higher S and Q have quantum numbers expected from multiple magnons and 
hole pairs. This finding confirms the basic assumption of the 5*0(5) theory, that the low energy collective degrees of 
freedom can be described by the superspin alone. 

At the next level, the pattern of symmetry can itself be tested. The level z/ of a given multiplet indicates the total 
number of magnons and hole pairs. If 50(5) symmetry is realized at a given chemical potential /ic, we would expect 
the free energy to depend only on v, the total number of magnons and hole pairs, but not on the difference between the 
number of magnons and hole pairs. As shown in Fig. 1261 the energy can depend on Q with three generic possibilities, 
similar to the discussions we presented in section HV. Al and Fig. O Only when the energy depends linearly on Q can 
the free energy be independent of Q at a given critical value of the chemical potential. From (|27ll we see that the 
energy levels indeed have this remarkable structure: states whose total charge differ by AQ = — 1 have nearly the 
same difference in energy. Therefore, the energy is approximately a linear function of Q or doping, similar to the 
situation discussed in section FV.BI To be more precise, the mean-level spacing within each multiplet (up to Q = ~2) 
is —2.9886 with a standard deviation of 0.0769. This standard deviation is much smaller (^ J/8) than the natural 
energy scale J of the t—J model and comparable to or even smaller than the average SC gap. Therefore, if one now 
adds the chemical potential term = —2^Q, and chooses ji = jic equal to the mean-level spacing, the superspin 
multiplets are nearly degenerate. Kt ^ — fic, magnons can be smoothly converted into hole pairs without free energy 
cost. This means that in each term of the expansion in (|90() one can freely substitute t\, or tj^ for t^ , and the direction 
of the superspin vector can be freely rotated from the AF to the SC direction. The smallness of the standard deviation 
indicates the flatness in the energy versus doping plot discussed in the previous section. If the standard deviation is 
significantly different from zero, this would indicate significant curvature in the energy versus doping plot. Therefore, 
the smallness of the standard deviation obtained by the exact diagonalization is consistent with the flatness of the 
energy versus doping plot obtained from the variational wave function discussed in th e previous sect ion. 

Another important aspect of the 50(5) symmetry is the Wigner-Eckart theorem ijGeorgil 119821) . This theorem 
provides a selection rule for the matrix elements of the operators based on the 50(5) symmetry of the system. It 
implies, for example, that the tt operators (see equation (|28|) 'l can only move us within a given multiplet, since they 
are symmetry generators. On the other hand, AF and d-wave SC order parameters (see equations (12 611 ) should mov e 
us between different multiplets. Both features have been verified in the numerical calculat ions in (Ed er et g/.Ul998(l . 

We conclude this subsection with a general remark. Exact diagonalizations (e.g. ()Dagottclll994|) ') commonly study 
ground-state correlations, but their spatial decay is often inconclusive as a test of order due to small system size. 
Discussions in this section show that it is possible that the (excited) eigenstates reveal a well-defined structure 
characteristic of a particular symmetry. Our strategy is to use the finite size calculations as input for effective models 
describing the collective degrees of freedom such as the superspin, or the magnons and the hole pairs. Since quantum 
Monte Carlo calculations can be performed for these models in a large size systems, the question of long range order 
and their competition can be firmly established. 
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D. Transformation from the microscopic model to effective S0(5) models 

From the two previous sections we have learned that both the variational wave function and the exact diagonal- 
ization of the t — J model show that the ground state and low lying excited states in the low doping range can be 
completel y described in terms of the sup erspin degree of freedom, with an approximate 50(5) symmetry. Altman and 
Auerbach ljAltman and Auerbachl l2002|) pioneered a systematic procedure in which they derived the effective bosonic 
50(5) model directly from the microscopic Hubbard and t — J models through a renormalization group transfor- 
mation called the Contractor Renormalization ("CORE") method(^'Iorningstar and Weinstein, 1996). This mapping 
has several distinct advantages. First, this approach helps to visualize clearly which processes and which excitations 
dominate the low energy physics of the system. Second, their work directly determines the parameters of the effective 
models defined in Eq. H50() and H65|l in terms of the microscopic parameters. The bosonic systems are often much 
easier to analyze numerically, as one does not have to worry about Pauli principles, Slater determinants, and the 
infamous sign problem in the quantum Monte Carlo algorithms. In this section, we shall describe their work. 

Since we want to construct bosonic quasiparticles, we have to divide the lattice into effective sites containing an 
even number of elementary sites (with one electron per site). In order to conserve the symmetry between the x and 
y-direction in the system, the original projected 50(5) model was formulated on a plaquette of 2x2 elementary sites. 
First we begin with the low energy eigenstates of the Heisenberg plaquette, which are determined easily. We find the 
nondegenerate ground state jfi) (see Fig. |Slfor a real-space representation in terms of the microscopic states on a 
plaquette) with energy Eq — —2 J and total spin 5 = 0. This singlet state will be the vacuum state of the effective 
bosonic projected 50(5) model. The next energy eigenstates are three triplet states with energy Et = —J and 

spin quantum numbers of 5 = 1. All other energy eigenstates of the Heisenberg plaquette have higher energies and can 
be neglected in the low energy effective model. It should be noted that the quasiparticles ta, which carry spin 1 and 
charge 0, are hardcore bosons because one cannot create more than one of them simultaneously on a sin gle plaquette. 

In their CORE study of the 2D Hubbard model, Altman and Auerbach (lAltman and Auerbachl l2002i) started from 
the spectrum of lowest-energy eigenstates of the 2x2 plaquette for 0, 1 and 2 holes, respectively. The corresponding 
lowest spectrum of the triplet (tj^), pair boson (tj^) and fermionic excitations is presented in Fig.|2Hl The ground state 
of two holes, also depicted in Fig. |H1 is described by 



where dji is +1 (-1) on vertical (horizontal) bonds within a plaquette and . . . stands for higher-order (U/t)-operators. 
Zf, is the wavefunction normalization. We note that ij^ creates a "Cooper" -like hole pair with internal d-wave symmetry 
with respect to the vacuum. The crucial point here is that while there is no hole pair binding for the Hubbard model 
on a dimer, there is bin ding in the range oi U/t G (0, 5) for a plaquette, a rather well-known fact (see, for instance, 
Rcf.(Hirs ch et aZ.l Il988|) '). However, this does not guarantee the integrity of the pair binding on the infinity lattice, 
documented by the fact that the hopping energy t is much larger than the pair binding energy, nor does it guarantee 
long-range SC order. To get more insight into these questions, one has to construct H^ff via a CORE procedure. 

In order to understand how the triplets and pair bosons behave on the infinite lattice, we must determine the 
boson hopping energies and the corresponding effective Hamiltonian. A suitable approach for this has been suggested 
by Mo rningstar and Weinstein on t he basis of the CORE technique, which has been shown for the 2D Hubbard 
model man and Auerbach, "200^, t ~ J ladders ^ Capponi and Poilbland |2002|) and earlier for Heisenberg chains 
and ladders (Mor ningstar a nd Weinstein, 1996) to be extremely accurate. For example, Morningstar and Weinstein 
obtained a very accurate ID Heisenberg model ground-state energy. This is even more impressive considering the 
latter model has long-range, power-law decaying spin correlations. 

In order to implement the CORE technique, the lattice is decomposed in small block units, as shown in Fig. 1301 
where Hq is the intra-block Hamiltonian and V is the part describing the coupling between the two neighboring blocks. 
The M low-energy states {|Q!°)}f^ are kept in each block i (here M = 4 in the 2x2 plaquette i) to define a reduced 
Hilbert space. The full Hamiltonian is then diagonalized on N connected units (in our example in Fig. 1201 = 2), 
i. e. for the (superblock) Hamiltonian Hs- The (in our case, — 4^) lowest energy states |5'„) with energy 
e„, n = 1 . . .M^ are retained. These true eigenstates of the N{= 2) block problem, {j^I^n)}, are then projected onto 
the reduced Hilbert space spanned by the tensorial product |a" . . . a'^) of the Af(= 4) states of each block, i. e. 




(91) 




(92) 



ai...Qjv 



and Gram-Schmidt orthonormalized, finally yielding the states {|^'n)} 



Then, the new superblock (renormalized) 
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Hamiltonian is defined as 

= |^«)(«'„| . (93) 

n 

By construction Hs, has the same eigenvalues e„ as Hs for n = 1, . . . , M^. Having constructed the new superblock 
or renormalized Hamitonian Hs, one can write (in our N = 2 example) 

Hs = Ho<»l + l<»Ho + V, (94) 

where Hq is simply the projected block Hamiltonian: 

M 

71=1 

The above equation H94|l gives the new renormalized interblock coupling V , restricted to the reduced Hilbert space. 
In the next step, one repeats the above procedure, replacing Hq and V in the original superblock Hamiltonian Hs by 
Hq and V , and so on. 

The projection onto the original plaquctte product basis in the Eq. H92|l expresses, of course, the above-discussed 
proliferation and possibly spatial decay of the block excitations. More generally, this is incorporated within the CORE 
method, in a superblock consisting of N blocks and a corresponding Hamiltonian containing A^-body interactions. 
The construction to obtain V (Eq. (jHH)) is different and, obviously, one also obtains y-terms, connecting N clusters 
instead of just N = 2 (called range- A^-approximation). It has been shown in the above-cited various applications that 
the above range-2 approximation {N — 2) and at most A'^ = 3 interactions already yield very accurate results. Thus, 
with a pr oper and physically motivated choice of the truncated basis, range-N interactions decay rapidly with N . 

In Ref . (|Ait man and Auerbachil2?ffll . the CORE calculation was limited to range-2 boson (triplets, hole pair boson) 
interactions leaving out the fermion state. From the above Fig. 1201 i- e. it is clear that this amounts to diagonalizing 
two coupled (2 x 2) plaquettes, for instance, an 8-site Hubbard cluster, which is very straightforward by the Lanczos 
technique. The resulting effective Hamiltonian for this range-2 four-boson model is exactly the projected 5*0(5) model 
defined in Eq. (|5()|l plus more extended interactions defined in Eq. H65(l . Following Alt man and Auerbach, we compar e 
in Fig. 13 II the magnitudes of the magnon hopping (denoted as Ji/2 ~ Jtt/2 in Ref. ijAltman and AuerbachLl200l 'l 
and the hole pair hopping Jc (denoted as Jb in Ref. (Altman and Auerbach, 2 00^ ') for a range of (t//t)-values. 

The first observation is that Jt ^ Jtt ~ 0.6 J; therefore, the triplet terms have a s i milar magnitude as those 
previously (see also our simple pedagogical Heisenberg example) obtained ijGopalan et alV 119941: ISachdev and BhatU 
[1990). 

The second finding is crucial. The region of equal Jt {Jtt) and Jh, equal magnon- and pair-boson hopping, occurs 
very close to U/t = 8. Thus, the value of the projected SO{5) model with Jt = Jb occurs in the physically relevant 
regime: It is known from a large body of numerically essentially exact (for example QMC) evaluations of the 2D 
Hubbard model t hat it reproduce s salient features of the HTSC cuprates precisely in this regime (see for example 
(Dagotto, 1994: Imada et al\ . \l99^ . This gives yet another piece of evidence, in addition to those discussed in sections 
IV. HI and IV. CI that realistic microscopic models can be described effectively by the projected 50(5) model close to 
the symmetric point. 

\ • : ■ - and Auerbach l)Altman and Auerbachl |2002|) . Capponi and Poilblanc ljCapponi and PoilblancL I2002D also 
calculated the coefficient and terms on in Eq. (|65|l . which contains triplet-triplet, pair-pair and pair-triplet 

interactions. These interaction terms were found to be small compared to H^ and /f*, but their influence has yet 
to be studied in detail. They also estimated the truncation error of discarding range-3 terms which, for physically 
relevant [/-values, was found to be very small (1%). 

An issue still left open is the role of fermions. Altman and Auerbach have extended the above 4-boson model 
to a boson-fermion model by augmenting the bosons with single-hole fermions "by hand." This is certainly a first 
step. However, a consistent low-energy theory has to treat bosons and fermions within the CORE procedure on equal 
footings. It should, however, be noted that the short-range effects of the fermions on the effective boson couplings 
were included in the above range-2 calculation. Altman and Auerbach estimated the fermion-boson interaction by 
including the hole fermions dispersion "adhoc" , i. e. using the single-hole band-structure extracted previously by 
various groups for large clusters (Dagotto, 1994J. 

In summary, the application of the CORE algorithm to the Hubbard model has demonstrated two features which 
are of immediate relevance for the 50(5) theory: the d-wave hole pairs already present in the 2D Hubbard model on 
a single 2x2 plaquette maintain their integrity in the "infinite" square lattice. The low energy degrees of freedom are 
indeed described solely by the superspin. Secondly, the hole-pair and magnon (triplet) hopping fulfills the projected 
50(5) condition in the physically relevant (U/t)-range. 
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VI. PHYSICS OF THE tt RESONANCE AND THE MICROSCOPIC MECHANISM OF SUPERCONDUCTIVITY 

A key experimental manifestation of a higher symmetry is the emergence of new particles or new collective modes. 
Historically, this line of reasoning has led to important discoveries in particle physics. For example, Gell-Mann used 
the SU{3) symmetry of the strong interaction to predict the fl~ resonance. Similarly, the electro-weak unification 
based on the SU{2) x U{1) symmetry has led to the prediction of the and the Z bosons. In a strongly interacting 
system, whether in particle physics or in condensed matter physics, typical excitations have short life times and broad 
line shapes. However, if higher symmetries are present, the selection rules associated with the symmetry prevents 
the excitation from decaying. The SO{5) symmetry of antiferromagnetism and superconductivity naturally predicts 
a new class of collective excitations, called the tt resonance or tt mode for short, which are the (pseudo-) Goldstone 
modes of the spontaneous symm etry breaking. The tt resonance can be identified naturally with the neutron resonance 
observed in the HTSC cuprates l)Fong et a/.l.ll995tlMo'ok adll993tlRossat-Mignod e^"a?] . ll991a(l . In this section we 
will review basic experimental facts about such resonances and discuss a theoretical scenario in which they originate 
from the pseudo Goldstone modes associated with the SO {5) symmetry. The operator of the 7r-mode is a symmetry 
generator of the 5*0(5) symmetry, so the appearance of the low lying resonance tells us about a small energy difference 
between the c?-wave SC and AF ground states of the doped cuprates. The idea of the near degeneracy of the d-wave 
SC and AF states lies at the heart of the 50(5) approach, which assumes that fluctuations between these two states 
exhaust the low energy sector of the system. Hence, experimental observation of the low lying resonances provide a 
key foundation to the 50(5) approach to competing AF and SC in the cuprates. In this section we provide several 
perspectives on the 7r-excitations. First, we use the 5*0(5) non-linear sigma model to describe them as pseudo- 
Goldstone modes of the approximate 5*0(5) symmetry of the system. Second, we show that the Fermi liquid analysis 
of the weakly interacting electron gas in a two dimensional tight binding lattice produces the 7r-mode as a sharp 
collective mode and gives a simple picture of this excitation as an anti-bound state of two electrons with the total 
spin 5=1 and with the center of mass momentum H = (tt, tt). Such excitation contributes to the spin fluctuation 
spectrum, measured by the inelastic neutron scattering, only in the SC state when there is a condensate of Cooper 
pairs. Finally, we discuss an important role that the 7r-resonance plays in stabilizing the SC state. 



A. Key experimental facts 

Resonant peak in the SC state of the cuprates was fir st observed in optimally doped Y Ba^Cu^^Or 
l|Rossat-Mignod et ad ri991alll992lll991bj) . Further experiments l|Fong et ad Il995t iMook et alUl99^ established 
that this is a magnetic resonance (spin 5* = 1) at the AF wavevector H = (tt, tt) which appears in the SC state. It has 
a constant energy luq = 41 meV at all temperatures and intensity that is strongly t emperatur e dependent and vanishes 
at Tc- Simi lar resonances have then been found in underdoped YBaoCu-iOn ^^ l |Dai et all. Il996l Il998t Fong et al\ . 
HOOO, 19961 iMook et all Il998t Istock et all l2003ft and in Bi2Sr2CaCu208+s l|Fong et all ll999HHe et all l200ll) and 
TloBa-yCuOf,^^ (|He et a/.ll2002l) . 

An important feature of magnetic scattering in underdoped YBa2Cu30e+x (|Dai et all Il996l Il998t iFong et aUi 
l2000t iMook et all I 1998f) is that the resonance precursors are detectable above Tc- Magnetic correlations, however, 
are strongly enh anced in the SC st ate, and there is a cusp in the temperature dependence of the resonant scattering 
intensity at Tc ijFong et all |2000() . Doping dependence of the resonance energy and intensity indicate a strong 
enhancement of magnetic fluctuations as we approach half-fillin g: for underdoped Y Ba2Cu^O(,j^x the resonance energy 
decreases with decreasing doping, and t he intensity increases (IFong et For overdoped Bi2Sr2CaCu20s+s 

it was found that the energy decrease d liHe et a/l l20q2l l200l[l . which led to a suggestion that the resonance energy 
follows the SC transition temperature ijHe et a^.l . l200l|) . 

The presence of the magnetic resonance in the SC state of many cup rates suggests th at it is closely related to the 
SC pairing. This idea was reenforced by the experiments of Dai et.al ijPai et all l200Cl|) . in which the SC coherence 
in Y Ba2Cu30e,e was suppressed by applying a magnetic field. It was found that the resonance intensity decreased 
without any noticeable change in the resonance energy. Finally, in (Dai et al, 1999) it was demonstrated that the 
exchange energy associated with the resonance has the right magnitude, the temperature and doping dependences to 
describe the SC condensation energy of YBa2Cu30e+x materials. 



B. Contribution of the tt resonance to the spin correlation function 

Resonance that appears in the SC state suggests that what gets scattered is Cooper pairs which are only present 
below Tc- Based on this idea it was proposed (Demlcr and Zhang, 1995) that the resonance observed in inelastic 
neutron scattering experiments is due to the presence of the 7r-mode, a sharp collective mode in the particle-particle 
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channel at momentum 11 = (tt, tt) with spin S = 1. In the normal state such excitation does not contribute to 
the magnetic spectrum, since the latter is determined by fluctuations in the particle- hole c hannel. Belo w T^, on 
the other hand, condense d Cooper pairs couple the particle-hole and particle-particle channels l)Demler et a l. 1998^ 
iDemler and Zhand . Il995|) and cause the tt excitation to appear as a sharp resonance in the magnetic spectrum with 
intensity set by the strength of mixing of the two channels, |A(T)p, where A(r) is the amplitude of the SC order 
parameter. Such a scenario provides a natural explanation for the key properties of the observe d resonance: its energ y 
is essentially the energy of the tt mode in the normal state and is temperature independent ijDemler et a/.l ri998hfl . 
whereas the intensity of the resonance is set by |A(r)p and vanishes at Tc. 

Coupling of the particle-particle 7r-channel and the particle-hole AF channel may be understood using the com- 
mutation relations between the operators tTq, and given in equation (|18|) . In the SC state the d-wave SC order 
parameter that enters the right hand side of equation (|18|) can be replaced by its expectation value in the ground 
state. Hence, the commutator of tt and N becomes a c-number, and the two fields become conjugate variables, just as 
coordinate and momentum are conjugate to each other in elementary quantum mechanics. The result of such coupling 
is that the 7r-mode appears as a sharp resonance in the spin fluctuation spectrum. To demonstrate this we consider 
the spin-spin correlation function at wavevector 11 

X(n,c.) = -z / e-^-HTN^m^mdt^J2\(^^\N^\^)f( 1^ _ ^ \, (96) 

J ^ \^UJ - En + l() LO + En - ] 

Here |0) is the ground state and n-summation goes over all excited states of the system. One of the excited states is 
created by the 7r-operator defined in equation (|28|l 

ka> - ^^im, (97) 

where M is the normalization factor. 

It is useful to realize that if tt''' acting on the ground state creates a n excited state, then tt sh ould annihilate it 
(otherwise it would create a state of lower energy than the ground state ijPines and Nozierej . 11966) ). Then we have 

1 = (tt^Itt^) = -i^(0|7r„^t |o) ^ ^(0| [7r„,7rt] |0) « ^i^, (98) 

where n is the filling fraction [n — 1 corresponds to half-filling) . In writing the last equality we assumed {{g{p)Y) — 1 
when averaged around the Fermi surface. 

If we separate the contribution of the tt state to x(H, lo) we have 

x(H,cj) = |(7r„|7V„|0)|2- ^ + part regular at 60^. (99) 

The resonant contribution to x(n, a;) can be expressed as 

x'^'=^(H,c.) = -L|(0K7V„|0)p^ 1 |(0|[,r,iV]|0)r 



mm? 1 

(1 — n) (tj — cJtt + «0) 



(100) 



The expectation value in the numerator of the last expression is simply the amplitude of the superonducting d-wave 
order parameter. We emphasize that Eq. (|100|l does not rely on the details of the microscopic model but only on the 
com mutation relations between the tt, N , and A given by the equation H18|l (this is somewhat analogous to the /-sum 
rule ijPines and Nozieresl Il96^ ). To relate the order parameter to what one typically measu res in experiments we 
use BCS type arguments t o connect the order parameter to the quasiparticle gap (see however (jEmerv and Kivelsonl 
119951 lUemura et alV IT981 ) (0|A|0) = CAo/Vscs- Here Ao is the maximal gap for Bogoliubov quasiparticles at the 
antinodal point, Vbcs is the interaction strength that we expect to be of the order of the nearest-neighbor exchange 
coupling J, and C is a dimensionless constant of the order of unity. Therefore, we find 

X-^(H,c.) = J^i^ L^. (101) 

J^(l — n) LJ — LJ^ -I- zO 

As we go to the underdoped regime, Aq remains constant or increases slightly, and the factor 1 — n decreases. Eq. 
(|101|l pr edicts that t he int ensity of the resonance should increase; this increase has been observed in the experiments 
inRef. l|Fong et adl200(]|) . 
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It is useful to note that contributions from modes other than the tt excitation do not spoil the result in H100() . If 
most of the 7r-spectrum is accommodated in an in terval (oj^ — v^uJt^ + v '), one can use the Cauchy-Schwarz inequality 
to prove a rigorous and model-independent result l|Demler et all . ll998h|) that 

-/ dtjlmx''°'(n,cj) > -L-^. (102) 

The left hand side of this equation represents the contribution of the tt mode to the spin excitation spectrum and 
the right hand side gives its lower bound. Exact equality holds when tt operator is an exact eigenoperator and hence 
there is only one energy eigenstate which satisfies (0|7ra|?^) 7^ 0. 

Thus, a simple picture of the resonant neutron scattering is as follows: when an incoming neutron is scattered off 
one of the electrons in a Cooper pair, it transfers a momentum of (tt, tt) to this electron and flips its spin. At the end 
of the scattering process the Cooper pair has quantum numbers of the 7r-mode, spin S = 1 and momentum 11. If the 
energy transfer matches the energy of the tt excitation, we have a resonance. In the next two sections we build upon 
this simple argument to establish a more detailed picture of the 7r-resonance in the two cases - the strong coupling 
limit described by the 5*0(5) non-linear a model, and the weak coupling limit where the Fermi-liquid type analysis 
can be applied. 



C. TT-resonance in the strong coupling: the S0(5) non-linear a model and the projected S0(5) model 

In this section we review how the resonant peak observed in the inelastic neutron scattering experiments appears in 
the 5*0(5) non-linear a model, signalling competition between the AF and SC ground states. We use the Hamiltonian 
of this model (see equation H34|) ') to write the operator equations of motion (O = i[7i, O]) for the order parameters, 
Ua, and symmetry generators. Lab, with a,b = {1, . . . , 5}. For > fic ~ \\/ qIx the system is in the SC ground state, 
which we can take to be along the n\ direction. Linearizing equations of motion around n\ we obtain 

XdW = P^lns (103) 
Xd^no, = pVln^ - [x{2^if - g]n^. (104) 

The first equation describes the Goldstone mode of the spontaneously broken charge U(l) symmetry (Bogoliubov- 
Anderson mode) and the second equation describes a triplet massive excitation of the superspin in the direction of 
the AF state (see Fig. 1^ . 

In a model with exact 50(5) symmetry superspin ordering reduces the symmetry from 50(5) to 50(4) and should 
be accompanied by the appearance of four Goldstone modes (50(5) and 50(4) have ten and six symmetry generators 
respectively). In the case of approximate 50(5) symmetry that we discuss here, explicit symmetry breaking turns 
some of the Goldstone modes into pseudo-Goldstone excitations, i.e. they acquire a finite energy. This is similar to 
a chiral symmetry breaking in quantum chroniodynamics, where a smal l mass of the qu arks leads to a finite mass of 
pions, which are the Goldstone bosons of the chiral symmetry breaking l)Weinberelll995(l . but it does not change the 
fundamental nature of the latter. 

The doping dependence of the resonance energy follows immediately from the equation (|104|l 

uj^ = 2V/i'-Mc- (105) 

The resonance energy is zero at the 5*0(5) symmetric point fj, = Hc and increases with doping. Vanishing of the 
resonance energy at Hc is a special property of the 5*0(5) symmetric point, and for a generic first order transition 
between the AF and SC phases (see Fig. llOb ) the resonance energy would remain finite at the transition point. 
When there is an intermidiate uniform mixed AF/SC phase ("type 2" transition shown in Fig. El(c)), the doping 
dependence of the resonance energy also obeys ()105f) with fi^ corresponding to the boundary between the SC and 
AF/SC phases (/ic2 in Fig. HOD. Softening of the tt mode in this case demonstrates a continuous transition into a 
state with magnetic order ( Demler et all I2OOII : ISachdev and M.Voita. 200Q') . 

I ■ ■ dispersion of the it resonance mode is model dependent. Hu and Zhang ljHu and ZhaneL l200l[l studied the 
dispersion of the tt resonance mode in the projected 5*0(5) model using the strong coupling expansion, and concluded 
that the tt mode can have a downward dispersion away from the 11 point, reaching a minimum at some incommen- 
surate wave vector. This model could possibly give a unified description of the neutron resonance mode and the 
incommensurate magnetic fluctuations in the HTSC cuprates. 

In Section [111. CI we discussed the projected 50(5) model that forbids double occupancy of the Cooper pairs by 
introducing chirality into SC rotations. As was pointed out before, such a projection does not affect small fluctuations 
around the SC state (see Fig. [S]) and does not change the relation H105|) . 
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D. TT-resonance in weak coupling: the Fermi liquid analysis 

In this section we consider a weakly interacting electron gas in a two dimensional square lattice and show that the 
Fermi liquid analysis of this system gives rise to the 7r-mode that is very similar to the collective mode we discussed 
earlier in the strong coupling limit. Using perturbative Fermi liquid analysis to describe strongly interacting electron 
systems, such as cuprates, may cause reasonable objections from some readers. We remind the reader, however, that 
the goal of this exercise is to complement strong coupling discussion presented in the earlier sections. The benefit of 
the weak coupling discussion is that it provides a simple picture of the 7r-mode as an anti-bound state of two electrons 
in the spin triplet state having a center of mass momentum H and sitting on the neighboring lattice sites. 

Our starting point is the t-J type model on a two-dimensional lattice 

{ij)a i (ij) 

Note that we do not impose a no-double occupancy constraint but include the on-site Hubbard repulsion. Within 
the Hartree-Fock discussion presented here, the Hubbard U only renormalizes the band structure, but it does not 
affect collective excitations of the order of J. Therefore, in the rest of the paper we will disregard the U term in the 
Hamiltonian H106|l and assume that we work with the renormalized parameters. 

To begin, we consider adding two non-interacting electrons into an empty two dimensional lattice with the condition 
that the center of mass of the pair has momentum q. For a general q the energy of such a pair, given by eq-k + 
Cfe, depends on the relative momentum of the two electrons. Therefore, we have a continuum of particle-particle 
excitations. When the center of mass momentum is H = (tt, tt) the whole particle-particle continuum collapses to a 
point. This can be verified by taking the tight binding dispersion — —2t{coskx -\-cosky) and is shown schematically 
in Fig. 1331 The collapse of the continuum makes it easier to create resonant states by adding interaction between the 
electrons. For example, the J term in the Hamiltonian H10t)|) introduces an energy cost of J/4 for electrons sitting 
on the nearest neighbor sites when their spins point in the same direction. Thus, if we make a two electron pair 
in such a way that the two electrons form a triplet pair on the nearest neighbor sites and have a center of mass 
momentum (7r,7r), we get an anti-bound state separated from the continuum by energy J/4. The argument above 
can be generalized to the case of adding two electron s on top of the filled Fermi sea. We recall that collective modes 
correspond to poles of the vertex functions lAbrikoso v et al\ . ^993) . In th e case of the 7r-resonance, we are interested 
in the particle-particle vertex, which we describe by the Dyson's equation f Pemler and Zhaii3 . Il995|) after separating 
the spin triplet component of the interaction at the center of mass momentum H with the d-wave symmetry of the 
electron pair 

"Hj = E '^!3^fec|,+nQ('^2CT)a/3clp^C_fcT.((TCr2)75Cfc+n5 + ■ ■ • (107) 

From the equation presented in Fig. 1341 we find the triplet particle-particle vertex 

T(p,p', H, Lo) = (108) 

and observe that it has a pole at energy 

= -2^+^(1 -n). (109) 

The first term in H109I) originates from the kinetic energy of the tight binding Hamiltonian €p + ep+n = and 
the second part of H109|) describes the nearest-neighbor exchange interaction of the triplet pair of electrons in the 
presence of a filled Fermi sea. The (1 — n) factor describes the blocking of the states below the Fermi energy from the 
phase space available for two particle scattering. In the Hartree-Fock theory the chemical potential is proportional to 
doping; hence, we find that the resonance energy in Eq. H109|l scales with x. It is useful to point out that including 
the near- neighbor density interaction ^y^(;j \ njiij in th e Hamiltonian (|106|l will not change our discussion as long as 
the system remains in the d-wave SC state (jPemlerTii a/.. . .1998b: .Meix ne r et alV Il997|) . Such an interaction affects 
equally the tt mode and Cooper pairs that constitute the ground state. 

One can also ask how to use the perturbative approach to demonstrate the appearance of the 7r-resonance in the 
spin-fluctuation spectrum below Tc. In Fig. I35l we show that when we compute the spin-spin correlation function in the 
SC state, we need to include scattering of spin fluctuations at momentum (tt, tt) into the tt pair, which corresponds to 
mixing the particle-particle ladder of diagrams into the particle-hole bubble. This contribution requires two anomalous 
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Green's functions and is therefore proportional to | Ap. Detailed calcula tions based on generalized random phase 
approximation for the model 1)1 U9|) were presented in ijDemler et a?J . IT998bf) : in Fig. 1201 we only show a representative 
plot of spin-spin correlation function x(g,ti)) computed with an account of the tt channel. 

In summary, we used a Fermi liquid analysis to establish a simple picture of the 7r-resonance as a triplet pair 
of electrons sitting on the nearest neighbor sites with the d-wave function of the pair and with the center of mass 
momentum IT. 



E. Resonance precursors in the underdoped regime 

In the unde rdoped cuprates t he resonance does not d i sappear above Tc b ut remains as a broad feature at higher 
temperatures l)Dai et all . Il996l Il998t iFong et all Il996t iMook et all Il998ll) . with only a cusp in t he te mperature 
dependence of the i ntensity signa l ling t he onset of the long range d-wave SC order l)Fong et all |200(]|) . In Ref. 
l)Demler and ZhaneL IT999bt IZhand. [l99l , it was pointed out that the most likely origin of these resonance precursors 
is the existence of strong c?-wave SC fluctuations in the pseudogap regime of the underdoped cuprates. Precursor of 
the TT-resonance in the spin-spin correlation function can be identified with a process in which a 7r-pair and a preformed 
Cooper pair propagate in opposite directions, as shown in Fig. |37| Because uncondensed Cooper pairs have a finite 
energy, we expect precursors to appear at a slightly higher energy than the resonance itself and have a width of the 
order of temperature (^Dcmlcr and Zhang> .1999b,'l . 



F. Implications for experiments and comparison to other theories 

In Section IVl. CI we discussed the 7r-resonance as a pseudo-Goldstone mode of the SO{5) non-linear cr-model, and 
in Section fVI.DI we gave a simple microscopic picture of the 7r-mode as a sharp collective mode in the particle-particle 
channel with spin S = 1 and momentum 11 = (tt, tt). From Eq. (|100|l . we see that the tt resonance intensity due to 
the contribution from the particle-particle channel scales with the square of the SC order parameter, namely 

/(n) = J dw/TOx''°'(n,cj) cx \{A{x,B,T))\^. (110) 

Here we have explicitly exhibited the dependence of the SC order parameter A(x,B,T) on doping x, magnetic 
field B and temperature T. Therefore, this simple scaling relation makes powerful predictions on the resonance 
intensity and has been tested in a number of experiments. Our analysis explains several puzzling features of the 
resonance observed in experiments. The first is the striking contrast betw een its temperature dependent intensity and 
temperature independent energy. Taking the Bardasis-Schrieffer exciton l*Bard asis and Schrieffeilll96l that appears 
as a bound state below the quasiparticle gap for s-wave superconductors, both energy and intensity of the exciton 
will be determined by the SC gap; hence, as the temperature is increased in the SC state, both the resonance energy 
and its intensity decrease. In the case of the 7r-mode, on the other hand, different behavior of the resonance intensity 
and energy are expected. The energy is essentially given by the energy of the 7r-mode in the normal state and does 
not change with temperature. The resonance intensity is set by the d-wave SC order parameter, as given in Eq. 
(|110|l . and decreases with increasing temperature and vanishes at Tc- Eq. (|110|l also predicted that the suppression 
of the SC coherence by a magnetic field should lead to a rapid decrease in the resonance intensity withou t changing 
the r esonance energy. This prediction was confirmed experimentally in a striking experiment by Dai et.al l)Dai et all 
|200(]|) . reproduced here in Fig. |2H1 The S0{5) theory predicted that with decreasing dopi ng the resonance intensity 
should increase (see Eq. (IIOII) ) and its energy should decrease (see Eqs. (|1U5|) and (|109|) ') llDemler and Zhan3.ll995t 
IZhanglll99'A . both of which have been observed in experiments, as we show in Fig. (JHl Note that for small values of 
the chemical potential there is a small difference in the precise vs relation obtained from the non-linear a model 
and the Fermi liquid analysis. We expect the strong coupling expression (|105|l to be more reliable close to the AF/SC 
transition where fi ^ fic and suggest that comparison of the doping dependence of the resonance energy l)Fong et al\ . 
I2OOO) and the chemical potential (Fujimori et ai, 1998 ; Ino et at, 1997) should be an important test of the SO{5) 
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ton that is overdamped in the normal state but becomes sharp in the d-wave SC state when a gap opens up for 
single particle excitations. In the d-wave SC state, the particle-particle channel and the particle-hole channels are 
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mixed into each other and there are, strictly speaking, no rigorous distinctions among these different theories. How- 
ever, important quantitative predictions differ in details. Near the transition, the tt resonance theory predicts a 
sharp onset of the magnetic resonance due to the coupling to the particle-particle channel, whose contribution to 
the magnetic scattering can be rigorously established via the Cauchy-Schwarz inequality, as shown in Eq. (|102|l . 
Some of these alternative theories expect a gradual broadening of the resonance rather than a sharp reduction of the 
intensity as Tc is approached from below. The tt resonance theory predicts that the energy of the magnetic resonance 
mode is independent of the temperature near Tc, while some of the alternative theories predict that the mode energy 
should vanish as the SC gap. In section I VII. El we shall discuss a rigorous distinction between the 7r-mode in the 
particle-particle channel and the magnetic exciton in the particle-hole channel in the normal state, and discuss an 
experimental proposal where this distinction can be tested. 

Se veral proposals h ave been made regarding implications of the resonance peak for various properties of the cuprates 
(see llKee et a^J . l2002|) for a critical review). S cattering of quasipart icles on the 7r-mode was argued to be responsible 
for the "kink" in the quasiparticle dispersion llJohnson et aZ.ll200l|) . "peak-dip- hump" structure m easured in ARPES 
(lAbanov e,t al\ . l200lt lEschrig and NormanL l200(]l) . and the pseu dogap seen in optica l conductivity llSchachinger et all . 
2001 ) . SC pairing mediated b y the resonance was suggested in l|Abanov et oi.l . Eofill ICarbotte et oi.Lll999HOrensteinL 
199S ; IZasadzin ski et a O-bOOSl). and relation between the resonance intensity and the condensate fraction was pointed 
out in Ref. ifchakraTOrtT^Lnd Keel 12000(1 . We do not discuss these proposals here, but in the next section we will 
review an important role that the resonance plays in thermodynamics of the SC state. We will argue that the SC 
condensation energ y may be accounted for b y lowering of the spin exchange energy due to the appearance of the 
resonance below Tc l|Demler and ZhangL[l99^ . 



G. Microscopic mechanism and the condensation energy 

The central question in the field of HTSC concerns the microscopic mechanism of superco nductivity. In conventi onal 
superconductors, the pairing interaction is mediated by the phonon interactions (see Ref . l|Maksimov et a/lll99? l for 
a review). Within the weak coupling BCS theory, the vertex corrections are suppressed by a small parameter, namely 
the ratio of the electron mass to the nuclei mass. Thus, the interaction which mediates the pairing of electrons 
can be unambiguously determined. In the case of HTSC, the dominant interaction is the Coulomb interaction and 
the AF exchange interaction. In such a strongly coupled system, the traditional approach based on the Feymann 
diagram expansion does not work, and the nature of the pairing interaction is not easily revealed by studying low order 
diagrams. However, the mechanism of superconductivity can still be addressed by identifying the interaction terms in 
the Hamiltonian which is lowered in the SC state. By comparing the magnitude of the energy saving associated with 
a particular interaction term with the actual experimental measurement of the condensation energy, the mechanism 
of superconductivity can be unambiguously identified. In our discussion in the previous section we showed that the tt- 
mode co ntributes to the spin fluctu ation spectrum below Tc and, therefore, enhances AF correlations in the SC state. 
In Ref. l|Demler and Zhaii3 . Il998|) . it was shown that the 7r-resonance can be promoted from being a consequence 
of superconductivity to being the real driving force behind the electron pairing. By analyzing the neutron scattering 
data, Demler and Zhang demonstrated that lowering of the AF exchange energy in the supercoducting state due to 
the appearance of the 7r-resonance can be sufficient to stabilize superconductivity in the first place. In this section 
we provide the details of this condensation energy argument focusing on the microscopic t-J model and discuss its 
relevance to the condensation energy of Y Ba2Cu^0Q+x materials. We also demonstrate that this scenario can be 
formulated as an additional contribution to the BCS coupling constant in weak coupling. 



1. The TT Resonance Contribution to the Condensation Energy 

The SC condensation energy is defined as the energy difference between the SC and the normal states at T = 
ISchricffcr, 1964; Tinkham, 1995). In type I superconductors it can be obtained directly by measuring the critical 
value of the magnetic field, at the first order transition between the normal and SC states. At the transition point, 
the energies of the two phases are equal (note that at T = the free energy is equal to the energy) and, assuming 
that the normal state is not affected by the magnetic field, we obtain the condensation energy per unit cell 

Ec^En-Es=^^, (111) 



where Vq = a x 5 x c is the volume of the unit cell. For type II superconductors including the HTSC, such a simple 
argument is not available. However, one can use LG theory to relate the condensation energy to Hci and Hc2, or 
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alternatively to the SC coherence length and London penetration depth A ijTinkhamL Il995j) : 

where $o = hc/2e is the SC flux quantum. An alternative approach to measuring the condensation energy is to 
integrate the difference between the SC and the normal state specific heat from T = to Tc, wh ere the norma l state 
speci fic heat below Tc is defined as extrapolation from temperatures above the transition point ijLoram et all 1 19941 
To be more precise, let us consider the condensation energy of the optimally doped YBa2Cu30T. Taking 
the characteristic values = 12 — 20A and A = 1300 — 1500A, with a = b = 3.85A and c = 11.63A, we find the 
condensation energy of Ec = 3.5 — 12 K per unit cell. The d etermination of the Ec of this material using specific 
heat measurements by Loram et.al. (|Loram et a^J . Il994l ll990') gave Ec ~ Q K per unit cell. 

Ideally, one would like to start with a microscopic model that has kinetic energy of electrons and ions, and the 
Coulomb energy of all particles, and calculate the condensation energy from first principles. Although possible in 
principle, in practice this approach is very hard to accomplish because of large scales involved in both the kinetic and 
the Coulomb energies. A method that is easier to pursue in practice is to start with an effective model defined on a 
much smaller energy scale and try to calcu late the condensation energy within this effective model. This approach has 
been undertaken by Scalapino and White ijScalaoino and Whitel IT99ql within the t-J model. In the t-J Hamiltonian 
in equation J^J, we have two terms: the kinetic energy of electrons (with the Gutzwiller projection operator) and the 
exchange energy of electrons. Analogous to conventional superconductors, we expect that the transition into the SC 
state is driven primarily by lowering the interaction part of the Hamiltonian , i.e. the exchange term (in conventional 
superconductors the relevant interaction is electron-ion Coulomb interaction). Is it possible then to find the change 
in the exchange energy between the normal and SC states? Scalapino and White made the insightful observation 
that the value of the J term in equation Q is directly related to the dynamic spin structure factor ^("{q^uj)^ the 
quantity that is being measured directly in neutron scattering experiments. And the change in the exchange energy, 
/S.Ej = Ej — Ej, can be directly expressed as a frequency and momentum integral of the difference in dynamic spin 
structure factors Xjv(<Z, w) — x'silj'^) with a form factor coming from the interaction being near neighbor 

AEj = 3J(^)2 f d'q ^^{x'l^{q^co)^xUQ,^Wo.siq,a)+cosiqya)). (113) 
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This equation applies to the quasi- two-dimensional systems, and q — {qx,qy) is a two dimensional in-plane mo- 
mentum. The generalizat ion of (|113|l to the bilayer systems, the case relevant for YBa2Cu3 06.35, is given in Rcf. 
l|Demler and Zhand. [l993) . 

The quantity 'x'^{q^uj) in equation (|113|) is not the normal state spin structure above Tc but rather an extrapolated 
normal state quantity at T = 0. Experimentally, one has to carefully identify features in ^("{q^uj) which change 
abruptly at Tc. From inelastic neutron scattering experiments we know that the most drastic change between the SC 
and the normal state spin structure factors is the appearance of the 41 meV resonance. Even for underdoped materials, 
which have many more AF fluctu ations in the norm al state, the main change between the normal and SC states is 
the appearance of the resonance ijFong et a^J . l2000|) . It is then reasonable to take formula (|113|) for AEj, calculate 
the contribution of the 7r- resonance a.nd arg ue that this will be the dominant contribution. For optimally doped 
YBa2Cu30Q,35 Fong et.al ijFong et a/.Ul996j ^ measured the absolute intensity of the resonance d{hLj)xs(^T^) to 
be 0.52 at T = 10 K. This resonance has a Gaussian profile centered at 11 with a width K2d = 0.23A~^, so the two 
dimensional integral can be easily estimated, and 

AEj = Ini^n^oA — = 0.016J. (114) 

Z Z Z TT 

Taking J = lOOmeV we find that the change in the exchange energy between the normal and SC states is approximately 
18 K per unit cell. This remarkable number tells us that the resonance alone can account for the SC condensation 
energy. 

Regarding our estimate of AEj in equatio n 11141) . a comment must be made. The dynamic spin structure factor 
S{q, w) satisfies the sum-rule llScalaoino et all il998il 

3(f)' d\ r ^x"(9, c.) = (1 - x)S{S + 1). (115) 

271" J-7r/a Jo 

Therefore, the spectral weight for the resonance needs to come from other regions in q-u) space. In obtaining H114|l 
we made an additional assumption that this weight was spread uniformly in q in x"((7,u;), and it did not contribute 
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to (|113|l . since any uniform component in is cancelled by the {cos{qxa) + cos{qya)) factor in Eq. ( I113II . It 

is also useful to point out that the weight of the resonance is less than 1% of the total sum rule l^Demler and Zhanj . 
ll998tlKee et oi.Ll2002^ . which, when multiplied by the AF exchange energy J, gives the correct order of magnitude 
for the condensation energy. 

The condensation energy argument can be generalized to finite temperatures. In this case the resonant peak 
intensity at temperature T should be related to the free energy difference between the SC and normal states, which in 
turn is given by the integral of the specific heat difference above T. This hypothesis has been analyzed by Dai et.al. 
fOai et aL, 1999), who showed that the temperature derivative of the resonant peak intensity follows very closely the 
specific heat anomaly for different dopings of YBa2Cu30e+x- We show this comparison in Fig. ^U] For optimal 
doping there is a BCS type anomaly in the specific heat at T^, which corresponds to the resonance appearing abruptly 
in the SC state. For underdoped samples the specific heat anomaly is broadened, which agrees with the resonance 
precursors appearing above Tc- This highly non-trivial experimental test establishes the connection between the tt 
resonance and its contribution to the condensation energy. 

Therefore, we see that the tt resonance mode naturally accounts for the condensation energy in the HTSC. The 
AF exchange interaction is lowered in the SC state, and this energy saving can drive the transition from the nor- 
mal state to the SC state. Within this scenario, the AF exchange energy is decreased, while the kinetic energy 
is increased below the SC transition. On the other hand, a number of theories argue that the domin a.nt driving 
mechanism of HT SC is the saving of the kinetic energy, either along the c axis, or in the Cu02 plane l) Andersonl 
[1997; Chakravar tv et all 11999: Hirsch and Marsiglio, 200fl}. The c axis kinetic energy saving mechanism has been 
definitively ruled out by experiments lMoler et al\. Il99^. The experimental measurement of the ab plane kinetic 
energy has not yielded conclusive results l)Keimeill2004HMolegraaf et all . l2002(l . The tt resonance based AF exchange 
energy saving is an experimentally established mechanism which can account for the condensation energy in the 
HTSC cu prates. Recent exper iments indicate that phonon mediated attraction also plays a role in the mechanism 
of HTSC l|Lanzara et all l200l[l . It is possible that various mechanisms contribute constructively to the condensa- 
tion energy in the HTSC. In this case, it is important to quantitatively measure the relative magnitudes of various 
contributions and identify the leading contribution to the condensation energy. 

2. Microscopic Discussions and Relation to the BCS Pairing 

In the theory of Ref. l|Demler and Zha'ii3 . ll99^ . the saving of the AF exchange energy arises from the coupling of 
the AF order parameter N to the tt operator in the SC state. This coupling leads to the additional spectral weight, 
proportional to |Ap, in the AF spin correlation function, thus lowering the AF exchange energy. This argument is 
generally valid, in both strong and weak coupling limits. However, it is also useful to connect this theory to the 
conventional BCS pairing theory in the weak coupling limit. In the limit of weakly interacting electron gas, we can 
formulate this scenario as a contribution to the BCS coupling in the d-wave channel. In Fig. ^2 we show a schematic 
representation of such a contribution: a Cooper pair splits into two virtual excitations - a magnon (iV) and a 7r-particle 
(tt) - which then recombine into a Cooper pair. One can easily verify that the quantum numbers are matched in this 
process: quantum numbers of the combination of the 7r-mode (charge 2, momentum H, spin S = 1) and the magnon 
(charge 0, momentum H, spin S = 1) sum to exactly the quantum numbers of the Cooper pair (charge 0, momentum 
q = 0, spin S* = 0). This may also be formulated using electron Greens functions, as shown in Fig. 021 We start with 
a Cooper pair formed by the electrons {p ]) and [~p [). After the latter electron emits a magnon, shown as an upper 
particle-hole ladder with total momentum H and spin Sz = —1, we have two electrons with momentum H and spin 
Sz = ^- These are exactly the quantum numbers of the 7r-mode that we describe by the lower particle-particle ladder 
in Fig. E21 

VII. KEY EXPERIMENTAL PREDICTIONS 

A. The antiferromagnetic vortex state 

A fundamental prediction of the 50(5) theory is the smooth rotation from the AF state to the SC state as the 
doping density is varied. As shown in sections IV. Bl and IV. CI this prediction has been tested numerically within the 
t — J model, with good agreement. However, testing this prediction directly in experiments would be much harder, 
since the doping level of most cuprates cannot be controlled w ell in the regime where the transitio n from the AF to 
SC state is expected to occur. Therefore, Zhang |)Zhangl [199^ and Arovas et al l)Arovas et aZ.lfl997j) proposed testing 
this prediction in the vortex state of underdoped cuprates. Around the center of the vortex core, the phase of the SC 
order parameter winds by ±27r, and the amplitude of the SC order parameter is constrained to vanish at the center 
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for topological reasons. In conventional BCS superconductors, the metallic Fermi liquid ground state is realized inside 
the vortex core. In the S0{5) theory, the SC order parameter is embedded as components of a higher dimensional 
order parameter, namely the superspin. When the amplitude of the SC order parameter vanishes in the vortex core, 
the amplitude of the superspin order parameter can still remain constant, provided that the superspin vector slowly 
rotates from the SC direction into the AF direction as the vortex core is approached. The superspin configuration 
near the vortex core is shown in Fig. 1431 This type of topological field configuration is known as the meron solution, 
meaning half of a Skyrmion( Raiaranianl^82). Fig. 1431 shows the rotation of the superspin in the vicinity of a vortex 
core. The AF order, which develops around the center of the vortex core, can be measured directly in experiments 
and can provide a quantitative test of the SO{5) symmetry. 

The key idea behind this prediction is more general. When the SC order is destroyed in the vortex core, the 
closest competing order develops in the vortex state. Aside from the commensurate or incommensurate magnetic and 
charge o rder, a number of novel correlation st ates have been proposed, including, for exarnple, the circulating orbital 
curre nts l)Chakravartv et all . l200ll: iLeel l2002|) and the fractionalized excitations ljSachdevl Il992t ISenthil and Fisheil 
Theref ore, the vortex core state c an provide a key test for various forms of the competing orders that have 
been proposed fSa chdev and Zheingl 12002'). 

Magnetic field provides a clean tuning parameter that can be used to investigate quantum transitions between the 
SC and AF phases. By solving both the 5*0(5) no n-linear sigma model and the LG model of competing AF and SC 
order parameters, Arovas et al ljArovas et al\ . \l99lf\ predicted the existence of the AF vortex state in the underdoped 
cuprates and further suggested a systematic experimental search for the AF vortex state in neutron scattering and 
muon spin rotation experiments. These authors also predicted that the magnetic field induced AF moment should 
increase linearly with the applied magnetic field, or the number of vortices in the system, when the applied magnetic 
field is small compared to the upper critical field Bc2- While the original analysis of Arovas et al focused on the 
regime where the transition between AF and SC is a direct first or der transition (corresponding to Fig. llOb of the 
phase diagram), Demler et a/ l)Demler et a^J . l200lHZhang et a^J . l2002(l analyzed the case in which there are two second 
order phase transitions with an intervening uniform AF/SC mixed phase, corresponding to Fi g. llOh of the phase 
diagram. In this case the AF order extends far beyond the vortex core region. Analysis in Refs. ijDemler et all l200lt 
IZhang et al. . 2002 ) demonstrated that the suppression of the SC order in this regime is dominated by the circulating 
super-currents and leads to a logarithmic correction to the linear dependence of the field induced moment. Recently, a 
number of experiments have been performed to test the prediction of AF order in the vortex state. Neutron scattering 
under a magnetic field can directly measure the field induced AF moment. Katano et al tKatano et a l.. 2000) measured 
enhanced magnetic scattering in the La2-xSrxCu04 crystal at a; = 12% doping. The i ntensity of elastic magnetic 
peaks around the (tt, tt) point increases at B = lOT by as much as 50%. Lake et al l|Lake et all l200lj) observed 
enhanced dynamic AF spin fluctuations in optimally doped La2-xSrxCu04 crystal at a; = 16% doping in an applied 
magnetic field. Without an applied field, the SC state has a spin gap of about 6meV. An applied field of B — 7T 
introduces a spectral weight in the energy range of 3 ~ AmeV. The mixed AF/SC phase has been also investigated 
in both the underdoped La2-xSrxCu04 crystal at x = 10% doping and i n the LaiCuO^+i, crystal . ^ In both material s 



the applied magnetic field strongly enhances the quasi-static AF ordering l)Khavkovich et a^J . 12002': 'Lake_et_£l 
The field dependence of the induced AF scattering is approximately linea r, as predicted in Re f. (Arovas et al 
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and it agrees quantitatively with the B\og{B/ Bc2) form proposed in Ref. ijPemler et a/.ll2boij) . with the correct value 
of Bc2. Another method to measure the AF order is the nuclear magnetic resonance (NMR). In the vortex state, 
the magnetic field is distributed inhomogeneously over the sample, with the maxima centered at the vortex cores. 
Therefore, the NMR frequency correlates directly with the location of the nucleus in the vortex lattice. Using NMR 
on the ^''O nucleus of YBa2Cu307 under a magnetic field as high as 40T, Mitrovic et al ljMitrovic et a^J . l200lll200^ 
detected a sharp increase o f the l/TjT r ate ne ar the vortex core as the temperature i s lower ed, indicating enhanced 
AF ordering (see also Ref. ljCurro et all l2000() ). Kakuyanagi et al ljKakuvanagi "all . l2002(l performed Tl NMR in 
the Tl2Ba2CuOQ+s sample. Tl NMR provides a more direct test of the AF ordering, since ^°^TI nucleus is located 
directly above the Cu spins. The temperature dependence of the l/TiT rate shows that the AF spin correlation is 
significantly enhanced inside the vortex core, compared with regions outside. The last class of magnetic experiments 
we discuss is the muon spin rotation (/isR) experiments. When muons are stopped inside a solid, their spin precesses 
around the local magnetic field. Since the muon decays predominantly along the direction of its spi n, the spatial decay 
pattern yields direct information about the local magnetic field distribution in a soHd. Miller et al ljMiller et allhOO'j) 
performed a /isR experiment in the underdoped Y Ba2Cu30G,5 system under a magnetic field of i? = 4T. They found 
that the local magnetic field distribution has a staggered pattern, superimposed on a uniform decay away from the 
vortex core. The staggered magnetic field detected at the muon site is about ISGauss. All the experiments discussed 
above were carried out at fields far below the upper critical field Bc2 , which in hole doped materials typically exceeds 
60r. In order to establish the nature of the competing state, one has to perform experiments close to i?c2- This was 
achieved in recent neutron scattering experiments on the electron dop ed Ndi ^^Ceo.i^CuOi crystal in magnetic fields 
up to 14T, far above the upper critical field, Bc2 l|Kang et all\2003^ . Kang et al found field induced AF scattering 
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at (7r,7r,0) and observed that the AF moment scales approximately linearly with the applied field up to Bc2- The 
AF moment decreases with the magnetic field in the range between 3^2 an d lAT. Their expe rimental data and the 
theoretical fit are shown in Fig| 44l The experimental findings of Kang et al l)Kang et all l2fl0,'^ have been confirmed 
by Fujita et al ljFuiita et a/.l . l2003 h in a related, electron doped material Pri^xLaCexCuOi. While Ndi,^^Cei^,i]^CuOi 
material contains the magnetic Nd moment, the Pri-xLaCcxCuOi material studied by Fujita et al does not contain 
such magnetic ions, th us confirming that t he field induced AF moment cannot be due to any spurious effects associated 
with the Nd moments (iMang et a/l l2(¥)^ . As we shall see below, the wide field range of the neutron data enables 
quantitative comparisons with theoretical models. 

Since the original theoretical predict ion of the AF vortex state, tremend o us th e oretical prog r ess has been made 
on the subject of AF vortex latticesllAlam a et al l Il999t Andersen et all l200(]t iBruus all. Il999t IChen et al 



2002HChen and Tind. I200I iDemler et a(.l. l200lllFranz ei ad l2002aHFranz and Tesanovicll200lHGhosal all\ 200: 



Han a nd Lcc"."2000V'Hu and Zhang. '200^ lHul . ll999aHJuneau et a?J . l2002HKivelson et g/l 12002: 'Mortcnscn et al.. 200C ; 
Ogatatil99& .Zhang et a/.„.2002 ). Based on the variational solution of the t — J model, Ogata [Ogata..l999.) concluded 
that the vortex core has AF with an ordered moment about 10% of the full moment. This calculation established 
the microscopic basis of the AF vortex core. The initial AF vortex solutions were based on the static mean field 
theory. In the weak field regime where the vortex cores are separated far from each other, the enhanced AF order can 
be viewed either as dynamic fluctuations of the AF order parameter due to the finite size of the vortex core, or as 
the bulk AF fluctuation pulled below the spin gap and spatially bound near the vortex cores. This dynamic picture 
was developed in Refs.(Bruus et al., 1999; Dcmlcr et al., 2001; Hu and Zhang, 2002} and could apply to experiments 
by Lake et al in optimally dop ed LSCO. Classical Monte Carlo calculations of the 5*0(5) model also show the exis- 
tence of the AF vortex lattice l)Hxl Il999a h. While the original theory of the AF vortex state was developed for the 
commensurate AF order, it can also be generaliz ed to the case where the AF ordering w ave vector deviates from the 
(tt, tt) point, as in the case of the LSCO svsteml Hu and Zhangll2b02HZhang et a^J . l2002() . The AF ordering inside the 
vortex core has a profound effect on the electronic structures of the vortex, since it opens up an insulating-like energy 
gap inside the vortex core where the conventional SC gap vanishes. The conventional theory for d-wave vortices 
based on Bogoliubov-de Gennes (BdG) mean-field theory predicts a large and broad peak at the Fermi energy in the 
local density of states (LDOS), the so-called zero-energy peak (ZEP), and at the vortex core (| Wang and MacdonaldL 
[199I. However, scanning tunneling spectroscopy (STS) spectrum in BSC CO, giving the LDOS around the vortex 
core directly, shows only a small-double peak structure at energies 7 meVI Pan et a/.l l2000() . A similar situation was 
observed in YBCO compounds! Maggio aprile e t al., 1995). The suppression of the local density of s tates due to the 
AF ordering inside the vortex core could naturany explain this phenomenon ll Andersen et all 120001: |Ch en and Tingl 
l2002t lOgatal '19991^^ Howeve r, other forms of order, or the smallness of the core size could also offer alternative 
explanations! Tsuc hiur a etall |2003). 

While the experimental observation of the AF vortex state confirms a major prediction of the 50(5) theory, most of 
these experiments have not directly tested the symmetry between AF and SC in the strictest sense. In the following, we 
shall discuss two aspects of the AF vortex state which directly pertain to the 50(5) symmetry. The spatial variation 
of the AF and SC order parameters around vortex core lead to a region of space where both order parameters coexist. 
In this region, the tt order parameter, whose magnitude can be quantitatively predi cted by the 50(5) orthogonality 
relation in Eq. (|41|l . also develops. Ghosal, Kallin and Berlinskv ljChosal et aUl2002() have quantitatively verified this 
relationship from their numerical solution of the t — J model around the vortex core. It would be desirable to find a 
way to measure the tt order parameter and test this relation experimentally. 

The detailed experimental data now available up to Bc2 in electron doped cuprates allows for a quantitative test of 
the 50(5) symmetry. As discussed in section llV.AI within models of competing AF and SC order, a crucial test for 
the 50(5) symmetry is the relation = U1U2 for the quartic term in Eq. (|59|l . Deviation from the 50(5) relation 
determines the curvature of the ground state energy versus doping plot, which can be used to determin e the nature 
of the transition between the AF and SC states. Recently, Chen, Wu and Zhang ljChen et a/l l2003b(l numerically 
solved the LG model with competing AF and SC order in the vortex state and found that the deviation from the 
50(5) relation = also determines the curvature of the field induced AF moment versus the ma gnetic field 
plot for magnetic fields up to Bc2. The neutron scattering data obtained in the NCCO superconductors ( Kang et all 
12003) can be fitted by U12/U1U2 = 0.95, showing that this system only has a 5% deviation from the 50(5) symmetry. 
When the magnetic field exceeds i?c2, it causes canting of the spin moments, thereby reducing the AF moment 
while increasing the ferromagnetic moment. Thus the 50(5) theory quantitatively expl ains the experiment al data 
in the entire magnetic field range below 14T. The experimental results of Fujita et al ljFuiita et all l2003(l in the 
Pri-xLaCexCuOi material are quantitatively similar. We note that the mean-field analysis of the GL free energy 
does not include quantum fluctuations of the AF order (the first term in equation H34|l 'l. The latter should be important 
when the AF moments are strongly localized inside the vortex cores. We expect that proximity effect type coupling 
between neighboring AF vortices should be sufficient to suppress such fluctuations. 

In the above discussions we focused on the AF moments of static votices in the SC state. The 50(5) model has 



47 



also been extended to study thermally activated phase slips in one-dimensional wires llSheehv and Goldbard ITool . 
One can also construct a du al effect to the AF vortices: Goldbart and Sheehy proposed AF hedgehogs with SC cores 
in Ref. l|Goldbart and Sheehv . 1998). 



B. The pair density wave state 

In the quantum disordered phase of the SO{5) model, the hole pair bosons become localized, forming a pair 
density wave. Since the superfluid density is low and pairing is strong in the underdoped regime of HTSC cuprates, 
the pair density wave state competes with the d wave SC state. In the global phase diagram shown in Fig. ^1 
aside from the half-filled AF insulator, there are several possible pair density waver states surrounded by the SC 
phases. In contrast to the supeconducting state, which can be realized for any charge density, each pair density wave 
state has a preferred charge density, the dominant one being at doping level a; = 1/8. Since the projected 5*0(5) 
model is formulated on the plaquettes of the original lattice, the pair density wave naturally forms a checkerboard 
pattern, as depicted in Fig. 1451 This state has a rotationally symmetric charge periodicity of 4a x 4a near doping 
level X — 1/8. However, connecting period of charge modulation to hole density in realistic systems is not always 
straightforward. In most cases we find states that have both superconductivity and periodic density modulation. 
Hence, they may be best described as supersolids. Supersolid phases are compressible and can accomodate extra 
charge without changing the period. Expressed differently, the excess charge can always be taken by the superfluid 
part of the Cooper pair density without affecting the localize d part. The pair density wave state differs from the 
stripe state l)Emerv et all . ll999HZaanen and GunnarssonL Il989() . since it does not brea k the sym metry of 7r/2 lattice 
rotations. It is also distinct from the Wigner crystal of individual holes proposed in Ref. l)Fu et a^J . 12004). which should 
have a cha rge periodic i ty of VSa x \/8a at the same doping level. The pair density wave state was first proposed by 
Chen et al ljChen et a^J . 12002;) in the context of the 5* 0(5) theory of the vortex sta te. It also arises na turally from the 
plaqu ette boson approach of Altman and Auerbach ljAltman and AuerbachL |2002() . Podolsky et al. ijPodolskv et all 
|200J?) discussed how unconventional states with translational symmetry breaking, including the pair density wave 
sta te, can be detected in STM experiments. Relevance of this state to tunneling experiments has also been considered 
in ijAndersen et a^J . l2003t IChen et aZ.Ll2003:IVoitaLl200^ . 

As we see in the global phase diagram shown in Fig. 1131 the pair density wave state can be stabilized near doping 
of a; = 1/8, when the superfluid density (or the kinetic energy of the hole pairs) is small compared to interaction 
energy. This situation can be realized in the vortex core , near the impurities, in the underdoped cuprates or in the 
pseudo-gap phase. The STM experiments (|Hoffman et al] . l2002ft measuring the local density of sta tes near t he vortex 
core demonstrated a 4a x 4a checkerboard pattern, consistent with the hole pair checkerboard state ljChen et at, 200^ 
shown in Fig. 1451 The vortex core can be either positively or negatively charged, d epending on whether the bulk density 
is greater or smaller than that of the nearby pair density wave state l)Wu et g^J . 12002,) . For example, if the chemical 
potential is such that the bulk SC state is on the left (right) side of the i5 = 1/8 insulator, w e expect the vortex core 
to have more (less) hole density. The STM experiment of Howald et al fHowald et aZ.','2002^ sees a similar real space 
modulation wit hout the applied magn etic field, possibly induced by impurities ( McElroy et aL , 2003^) . More recently, 
Vershinin et al ijVershinin et all . l2004|) discovered a real space modulation of the density-of-states in the pseudo-gap 
phase above Tc- Enhancement of the translational symmetry breaking in the pseudogap regime of the cuprates has 
been proposed theoreti cally in Ref . (Sachdcv and Demler, 2003). The microscopic pictu re of thi s phenomenon has 
been studied in Ref. l|Chen et all i2"004|l using an extensio n of the formalism in Ref. l|Podolskv et al\. I200:| for 
the pseudogap regime. Analysis of Ref. ijChen et all l2004|) shows that the experimentally observed modulation is 
inconsistent with an ordinary site centered charge density wave and the corresponding modulation of the Hartree-Fock 
potential. However, the pair density wave state provides good agreement with the experimental data. 



C. Uniform mixed phase of antiferromagnetism and superconductivity 

The phase diagram obtained from the classical competition between the AF and SC states is shown in Fig. ^| We 
classified the phase transition broadly into three different types. The "type 1" transition involves a direct first order 
phase transition between the AF and the SC phases, and the "type 2" transition involves two second order phase 
transitions, with an intermediate phase which is a uniform mixture of AF and SC. The marginal "type 1.5" transition 
describes the special 5*0(5) symmetric case where the chemical potential remains constant in the entire uniform mix 
phase. Therefore, both "type 2" and "type 1.5" transitions predict a uniform mixed phase of AF and SC. 

Evidence for the AF/SC mixed phase exists in the excess oxygen doped La2CuOiJ^y material. Neutron scattering 
meas urement detects the onset of the AF or spin-density- wave orders at the same temperature as the SC T.^ l)Lee et alV . 
I1999D . This remarkable coincidence is the hallmark of a multi-critical point, which we shall return to later. Because 
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the La2Cu04+y system has an ordering wave vector similar to that of the La2-xSrxCuOi^ system, it should also 
be classified as the "class B3" trace in the global phase diagram of Fig. 1131 passing through the "1/8" Mott lobe. 
However, in this case, the Mott phase boundary likely belongs to "type 1.5 or 2," where the AF and SC order can 
coexist. 

For Y Ba2Cu'iOQ+x materials, static magneti c ordering extendin g to a; « 0.5 has been observed recently using 
muon spin rotatio n/rela xation measuremen ts in ijMiller eit oil l2f)(l^ . Preliminary neutron scattering experiments in 
l)Sidis et al 1 1200 ID and llMook et olV I2002D also reported magnetic ordering with a wavevector (tTjTt). Thus, in this 
case we have AF coexisting with SC without any additional charge order. However, it is unclear if two phases coexist 
uniformly in these materials. Assuming that future experiments verify the existence of a homogenious phase with AF 
and SC orders, we conclude that the phase diagram for Y Ba2Cu^OQ+x may be understood as moving along the Bl 
line line in Figure^] when the system avoids all the PDW lobes but only has AF and SC orders either separately or 
in a uniform mixed phase. 

Evidence for the mixed phase of superconductivity and antiferromagnetism has also been obtained recently in 
the five-layered HTSC cuprate HgBa2Ca4,CuzOy. In this system, the three inner layers are predominantly anti- 
ferromagnetic , while the two outer layers are predominantly superconducting. In a Cu-NMR study, Kotegawa et 
al (i,Kotegawa et a/.L l2004|) obtained firm evidence that the AF inner layers induce a small magnetic moment in the 
outer layers, establishing the case of an AF/SC uniform mixed phase in t his system. However , this type of AF/SC 
proximity effect was not observed in the artificially grown layer structures ijBozovic et a?J . l2003() . 

Above discussions show that there is evidence of a uniform mixed phase of AF and SC in the HTSC cuprates. On 
the other hand, microscopic probes such as STM l)Pan et alV l200lj) reveal electronic inhomogeneities characteristic 
of the "type 1" direct first order transition between AF and SC. Therefore, depending on material details, some 
HTSC compounds show AF/SC mixed phase, characteristic of the "type 2" behavior, while others show microscopic 
separation between these two phases, a characteristic more consistent with the "type 1" behavior. It is quite remarkable 
that such different physical effects can be obtained in materials that are so similar. A reasonable explanation is that 
the system is actually very close to the 50(5) symmetric point exhibiting "type 1.5" behavior. Only in this case can 
a slight variation tip the balance towards either the "type 1" or "type 2" behaviors. 

A genuine uniform mixed phase of AF and SC has been observed in several heavy- fermion systems in some regions 
of the pressure (P) versus temperature (T) phase diagram jKitaoka et a^J . l200lL 12003) . Recently, such coexistence was 
observed through NMR and NQR spectrum measurements in CeCu2{Sii-xGex)2 with a small concentration x = 0.01 
of Ge. In CeCu2Si2, SC coexists with slowly fiuctuating magnetic waves. However, for AF CeCu2Ge2, which has the 
same lattice and electronic structure as CeCu2-S'«2, it was found that a SC phase can be reached at a critical pressure 
Pc ~ 7.6GPa. Since CeCu2Si2 behaves at P = like CeCu2Ge2 at Pc, it is argued that SC in CeGu2Si2 occurs 
close to an AF phase at P = corresponding to a critical lattice density D — Dc- This appears to be the reason for 
the strong AF fluctuations at P = 0. A small concentration of Ge expands the unit-cell volume reducing D below 
Dc and is thus sufficient to pin the magnetic fluctuations and to produce AF long-range order within the SC phase. 
Noting that D = Dsi[l{VGe — Vsi)x/VGe] for Ge doping and that D inc reases with pressure, one can draw a combined 
phase diag ram as a funct i on of lattice density D ijKitaoka et a/.ll2002(l . 

In Ref. ijKitaoka et all l2002(l . it was shown that the phase diagram of Fig. 1461 could be understood in terms of 
an SO{5) superspin picture. This could suggest that SC in GeGu2Si2 could be mediated by the same magnetic 
interactions leading to the AF state in CeCu2{Sii-xGex)2- 

D. Global phase diagram and multi-critical points 

The 5*0(5) theory makes the key prediction of the existence of a multi-critical point where T/v and Tc intersect 
(see Fig. I10|l and the general topology of the global phase diagram in the space of quantum parameters (see Fig. I13|l . 
The goal of this section is to establish the connection between the theoretical quantum phase diagram proposed in 
Section IIV.BI and the experimental phase diagrams of various families of cuprates. The underlying assumption for 
making such a connection is that most of material specific properties can be absorbed into parameters of the effective 
Hamiltonian given in Eq. H5U|I and in Eq. 1651) . 

One of the most studied phase diagrams of the HTSC is for La2-xSrxCu04. The stri pe order's presence in these 
mate rials has been well documented by neutron scattering experiments (Wakim oto et alL 2000, 2001; Yamada et al\ . 
Il998|) . For less than 5% doping the system is in the insulating regime with diagonal stripes and for higher dopings the 
system is superconducting with coUinear stripes (see Fig. I47|l . It is natural to relate this family of cuprates to the B3 
trajectory on the J/V-fj, phase diagram shown in Fig. E| with increasing fj, the system goes through a hierarchy of 
states at fractional filling factors that correspond to insulating pair density wave states. Near these magic filling factors, 
the SC Tc drops dramatically, while magnetic ordering increases substantially. This is indeed the behavior observed in 
Fig. 07| As discussed in Section llV.BI the two possible patterns of charge ordering are checkerboard and stripes. In 
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the case of La2-xSrxCu04, stripe ordering may be stabilized by tilting the CuOq octahedron toward [100] tetragonal 
direction (parallel to the Cu-0 bonds). The phase diagram in Fig. E| predicts that the ordering wave- vectors take 
discrete values that correspond to different Mott insulating PDW lobes. For long range interactions PDW phases 
are very densely packed, so experimentally w e may observe a n almost continu ous dependence o f incommensuration 
on doping, such as the one discussed in Refs. ijYamada et "all, [l998,) and (,Wakimoto et aZ.Ll2000|) . However, different 
states in the hierarchy are not equivalent. For example, at 1/8 doping we have a very strong insulating phase which 
corresponds to insulating stripes or a simple checkerboard pattern of Cooper pairs (see Fig. I12|l . This may explain 
a famous "1/8 anomaly" in the vs doping relation for the La2-xSrxCuOi family of cuprates. Another strong 
PDW phase is for 1/16 doping, which may explain why supreconductivity disappears close to this filling (see Fig. HTjl . 
A "staircase" of ordering wavevectors for underdoped cuprates has also been discussed in the context of doping the 
spin-Peierls insulating phase in Refs. l|Sachdevl . l2002btfVoita a nd Sachdev, 1999). 

By adding another external parameter we can tune our system continuous ly between Bl and B3 trajectories. This 
was done in recent high pressure experiment s on LaiA8NdoASro.i2CuOi ijArumugam et all. I2OO2I: iLocauet et al\ . 
ll998l:ISato et oi.L 1200(1 l^keshita et oi.l . l200'^ . where the pressure of the order of 0.1 GPa was sufficient to suppress 
stripe ordering at 1/8 doping and stabilize the high temperature SC phase. Such pressure experiments correspond to 
moving up along the A2 path in Fig. 1131 Applying pressure along this path can directly induce a superconductor to 
insulator transition. 

In contrast to the LSCO family of HTSC cuprates, when one varies the carrier density in the YBCO or BSCO 
cuprates, there is no evidence for th e static charge order. In these materials, charge ordered PDW states can only be 
realized ar ound vortex c oreslHoffman e t aZ.U2002). when the effective Coop er pair kinetic energy is reduced, or near 
impurities l)Howald et ali 12002; .McElrov et all l2003t IVershini n eTd] . \2004} . Therefore, we identify these materials 
with the Bl trajectories in the global phase diagram of Fig. ^1 In this case, the AF/SC boundary can be either be 
"type 1" or "type 2". Given the evidence discussed in section lVII.CI these systems seem to be close to the "type 1.5" 
marginal case in between these two types of phase transitions, which means that they should have the approximate 
50(5) symmetry. 

Within the class of materials exhibiting the "i?l" type of trajectory in the global phase diagram, the 50(5) theory 
makes a distinct prediction of the finite temperature multi-critical point where Tc and intersect. An interesting issue 
discussed in Sees. IIV.'UI and HVl is the possibility of analyzing the critical properties of systems (such as many HTSC 
cuprates) showing a direct transition between an AF and a SC phase. In particular, measuring the critical exponent 
associated with various physical quantit ies near the bicritical AF-SC point can give information about the dimension 
of the symmetry group at the transition l)Hu and ZhaneLr2000(l . Unfortunately, in the HTSC cuprates, sample qualities 
are not high enough to enable a reliable measurement of the critical behavior near the multi-critical points discussed 
above. On the other hand, encouraging experimental evidence for an 50(5) bicritical point does exist in a class of 2D 
organic superconductors called BEDT salt. These material share most common physical properties with the cuprates, 
and the AF to S C transition can be induced by pressure. In particular, recent experiments on fc— [BEDT — TTF)2X 
l|KanodaL 1199 re vealed an interesting phase diagram in which and intersect each other at a bicritical point. 
Kanoda ("Ka nodal Il997j) measured the NMR relaxation rate 1 /Ti both in the AF and in the SC region near the 
bicritical point. Below a characteristic temperature T*, 1/Ti diverges towa rd the AF transition temperat ure, while 
it exhibits a spin-gap-like behavior on the SC side. Murakami and Nagaosa l)Murakami and Naeaosal l2000|) analyzed 
this experimental data in terms of a generalized LG model including both AF and SC fluctuations near the bicritical 
point. Their study concentrated on the dynamic critical phenomena, in particular the relaxation rate 1/Ti around the 
bicritical point. A detailed analysis of the data allowed the extraction of the corresponding critical exponent x. Before 
discussing the NMR line width, we would like to caution the r eaders that ther e is also a first order metal-insulator 
transition in addition to the AF/SC transition discussed here l)Lefebvre 6^*0? 1 12000). The presence of the critical 
endpoint of the metal insulator transition line may lead to some additional complications in the analysis. 

On the AF side of the phase diagram, the NMR line width is proportional to (T — Tm)~^ when approaching Tm 
from the normal state. For systems far away from the bicritical point, the dynamical critical behavior is governed by 
the 50(3) Heisenberg model, whose exponent x = 0:3 « 0.315. On the other hand, when the 50(5) bicritical point 
governs the critical dynamics, the exponent x should change to the 50(5) one x = x ^ ~ 0.584, as obtained from th e 
e-expansion. In Fig. HSl we present a lo g-log plot of 1/Ti vs (T - Tc)/T (from Ref. l|Murakami and Nagaosal l2000|) , 
data from Ref. l|Kawamoto et a/.Ul995|) ') for (A) /t-(BEDT-TTF)2Cu[N(CN)2]Cl (solid squares), and (B) deuterated 
K-(BEDT-TTF)2Cu[N(CN)2]Br (open squares). System (A) is located in the AF region away from the bicritical point 
and system (B) is nearly at the bicritical point. As one can see from the figure, the critical exponent, x, is 0.30 ±0.04 
for system (A) and 0.56±0.04 for system (B). These values of x are in reasonably good agreement with the theoretical 
ones, and, in particular, support the fact that the AF/SC bicritical point is governed by the 50(5) symmetric fixed 
point. This is the first experiment which directly measures the dimension of the symmetry group close to the AF/SC 
bicritical point and determines n to be close to 5. More extensive study near the critical region is certainly desired. 

A central issue of the HTSC cuprates concerns the phase boundary between the AF and SC phases. It is also in this 
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region that the 5*0(5) theory makes the most direct and distinct predictions. Current experiments discussed above 
seem consistent with the zero temperature and finite temperature phase diagrams presented in Fig. 1131 and Fig. IIUI 
However, detailed quantitative comparison is still lacking. As the material properties of the HTSC cuprates improve, 
direct quantitative tests of the 5*0(5) theory, such as those performed in the organic superconductors, may become 
possible. 

E. The particle-particle resonance mode in the normal state 

In this paper we discussed the scenario in which the resonance peak in the inela s tic neu tron scattering experiments 
l|Fong et all Il999t Il995l: iHe et a i. 2001; iMook et aA \im?t lR,ossat-Mignod et all Il991h() originates from the triplet 
TT mode in the particle-particle channel. This mode does not disappear above T^, but it ceases to contribute to the 
spin fluctuation spectrum, since the particle-particle and particle-hole channels are decoupled from each other in the 
normal state. An important question to ask is whether one can couple to the 7r-channel directly and establish the 
existence of the resonance already in the normal state. This cannot be done using conventional electro-magnetic 
probes, which all couple to the particle-hole channels only, but i t is possible using tunnelling experiments. Before we 
discuss the specific proposal of Bazaliy e^.aL l|Bazaliv et aJl\l991^ for detecting the n excitations, it is useful to remind 
the readers about earlier work on measuring pairing fluctuations in conventional superconductors above their transition 
temperature (Anderson and Goldman, 1970). As originally proposed by Scalapino (Scalapinol 119701) . the latter can 
be measured in a sandwich system of two superconductor SCi and SC2 with different transition temperatures in the 
regime Tc2 < T < Td- Resonant coupling between Cooper pairs from the superconductor SCi and the fluctuating 
pairing amplitude in SC2 leads to the peaks in the IV characteristics at voltages that correspond to half the energy of 
the "preformed" Cooper pairs in SC2 . The generalizati on of these t unnelling experiments for detecting the tt mode in 
the normal state of the cuprate has been suggested in l)Bazaliv et at . 1997) and is shown in Fig. 0^ In place of the 
SC2 region we now have some cuprate material that shows a resonance in the SC state, e.g. an underdoped YBCO, 
(electrode C in Fig I49|) . and in place of the SCi materials we have a different cuprate superconductor (electrode A 
in Fig I49|l with a higher transition temperature than material C. The sy stem should be i n the temperature regime 
< T < T/. The main difference with the set-up suggested in Ref. l|Scalapinol Il970|) is the presence of a thin 
layer of AF insulator between the A and C electrodes. The reason for this modiflcation is straightforward: we need to 
probe the tt channel in the C material that corresponds to the particle-particle mode with spin 5=1 and momentum 
n = (tt, tt), whereas the SC electrode A provides Cooper pairs with 5 = and momentum q = 0. If the two 
materials are connected as shown in Fig. 1491 a Cooper pair travelling across an AF layer B can emit a magnon, which 
converts this Cooper pair into a 7r-pair and allows resonant coupling between superconductor A and the 7r-channel 
of the "normal" electrode C. One expects to find a resonance in the IV characteristics of the junction with a peak 
in the tunnelling current at a voltage which is exactly half the energy of the 7r-resonance in the C electrode (note 
that this peak only appears when electrons are injected from A to C, so it appears on one side of the IV curve). A 
simple qualitative picture described above can be made more precise by considering a tunnelling Hamiltonian between 
materials A and C 



Here V is the applied voltage, the apa and Cka operators refer to the electronic operators in A and C with momenta 
p and k. The ratio of the spin flip matrix element T^j, to the direct matrix element T^j, is on the order of Asdw /U , 
where Asdw is the spin-density-wave gap of the AF insulating material B. The diagram responsible for the resonant 
contribution to the tunnelling current is shown in Fig. 1501 The triplet vertex F takes into account interactions needed 
to creat e a sharp 7r-resonanc e in the A electrode. The magnitude of the peak in the tunnelling current was estimated 
in Ref. l|Bazaliv et a^J . 0*997(1 to be 10 hA^jlY for a system of area 10~^ cm^. As argued in section IVTl it is not easy to 
distinguish the particle-hole and the particle-particle origin of the tt resonance below Tc since these two channels are 
mixed. Direct experimental detection of the triplet particle-particle mode in the normal state would give unambiguous 
evidence of the particle-particle nature of the tt resonance mode. 

F. Josephson effect in the SC/AF/SC junction 

When discussing the relationship between c?-wave SC and AF in the HTSC cuprates, one often flnds signatures of 
the nearby magnetic phase in experiments performed on the SC materials. An important question to ask is whether 
the AF insulating phase shows any signatures of the nearby SC state. An intriguing set of experiments that possibly 
provides such a demonstration is the long range proximity effect observed in insulating samples of Y Ba2Cuy,0Q-^x 




(116) 



pk(7 
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based materials coupled in the a-b plane directions l)Barner et all Il99ll: iDecca et al !. '200Cf; 'Hashimoto et al\. Il992t 
ISuzuki et oA\1994) . The AF/SC proximity effect was also observed by Kotee awa et al( | K otegawa a/.Ll200^ . how- 
ever, it seems to be absent in the case of artificially grown c-axis coupled lavers ljBozovic et flrn200^ . The appearance 
of the long range proximity effect is very natural from the point of view of the 5*0(5) theory, in which low energy 
degrees of freedom correspond to the order parameter rotation between the AF and SC configuration s. A theory of 
the lon g range proximity effect within the 5*0(5) non-linear sigma model has been developed in Eef . l|Demler et aJ\ . 
ll998aD . Let us consider the SC/AF/SC junction shown in Fig. If we set Re A = cos0cos(/), /to A = cos sin 0, 
and N3 = sin9, then according to our discussion in Section UlI.BI fsee Eqs. (jS^ - iP^ ). the junction can be described 
by the effective Lagrangian density 

C{9, 0) = ^{{d^ef + cos^ 9 ~ g sin^ 9. (117) 

The anisotropy term is given by 5^ > inside the A region, so that the AF phase would be established in the 
bulk. In the SC S regions on both sides o f the junction \ ye have gs < 0, and we should impose boundary conditions 
6* ^ as a; ^ ±00. As discussed in Ref. ljPemler et a/] . Il998a() . a simplified case corresponds to taking a "strong" 
superconductor limit for which 9(x — 0,d) = 0. The current phase relation can now be obtained by writing the 
Euler-Lagrange equations for the functional 1)117(1 at a fixed current. The maximal value of 9 reached aX x = d/2, 9o, 
is determined by the equation 

= 'jpl^^ (118) 

W2 + C0S2 6*0 

where K{k) is the the complete elliptic integral of the first kind, dimensionless current Us = IS,a, with / being the 
actual current through the junction and the characteristic length 



U = VpNm- (119) 

On the other hand, the equation for the phase difference across the junction, A^, is given by 

A$ = 2c.,^^^l£^=ni(-sin2 0o,fc). (120) 

y'LUf + COS^ 6*0 

Here Hi(n, k) is a complete elliptic integral of the third kind. Immediately, one can see that equation (|120|l describes 
two different kinds of behavior for d larger or smaller than dco = t^S.a- When d > dco we have a conventional proximity 
effect with /(A"I>) = /o(d)smA$ and Io{d) cx exp{—d/ £^a)- We observe, however, that the SC correlation length, S^a 
may be very long if the system is close to the 50(5) symmetric point {gA — + in equation (|119|l ). which corresponds 
to the long range proximity effect. When d < dco we get more intriguing behavior in H12(J|) . where for small currents 
the A region is uniformly superconducting, i.e. 9o = (proximity to a strong superconductor completely suppresses 
the AF order inside the A region), but when the current exceeds some critical value, the system goes into a state 
that has both d-wave SC and AF orders, i.e. < 9o < Tr/2. The resulting nontrival /(A$) are shown in in Fig. |^ 
We note that the analysis presented above does not take into account the long range part of the Coulomb interaction 
between electrons. This may become important for systems with sufficiently wide AF layers and lead to suppression 
of the proximity induced SC order in the AF layer. 

Several conseq uences of the non-sinuso idal behavior of the current-phase relation of the SAS junctions have been 
explored in Ref. ijden Hertog et a/1 Il999tl . including current-voltage characteristics in the presense of thermal fluctua- 
tions, Shapiro steps, and the Fraunhofer pattern. Decca et.al. (Decca et al.. 2000.) used near-field scanning tunnelling 
microscopy to photo-generate Josephson junctions in underdoped thin films of Y Ba2CuzOQj^x- They have verified 
a long range proximity effect through insulating layers but observed a conventional Fraunhofer pattern rather than 
the one pre dicted in (den Hertog et al., 199 9i). The geometry of their samples, however, is different from the system 
studied in l)Demler et al\ . Il998at Iden Hertog et all . Il99^ : the intermediate AF layer in their case is connected to 
large AF regions on both sides of the junctions, which suppresses rotat i on of the superspin into the SC direction. 

In a related context, Auerbach and Altman jAuerbach and AltmanL |2000|) applied the project 5*0(5) theory to 
predict multiple Andreev resonance peaks in the SC/AF/SC junctions. 



VIII. CONCLUSIONS 



In a large class of materials including the HTSC cuprates, the organic superconductors and the heavy fermion 
compounds, the AF and SC phases occur in close proximity to each other. The 5*0(5) theory is developed based on 
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the assumption that these two phases share a common microscopic origin and should be treated on equal footing. The 
S0{5) theory gives a coherent description of the rich global phase diagram of the HTSC cuprates and its low energy 
dynamics through a simple symmetry principle and a unified effective model based on a single quantum Hamiltonian. 
A number of theoretical predictions, including intensity dependence of the neutron resonance mode, the AF vortex 
state, the pair-density-wave state and the mixed phase of AF and SC have been verified experimentally. The theory 
also sheds light on the microscopic mechanism of superconductivity and quantitatively correlates the AF exchange 
energy with the condensation energy of superconductivity. However, the theory is still incomplete in many ways and 
lacks full quantitative predictive power. While the role of fermions is well understood within the exact 5*0(5) models, 
their roles in the projected 5*0(5) models are still not fully worked out. As a result, the theory has not made many 
predictions concerning the transport properties of these materials. 

Throughout the history of our quest for the basic laws of nature, symmetry principles have always been the faithful 
guiding light which time and again led us out of darkness. The enigma of HTSC poses an unprecedented challenge 
in condensed matter physics. Reflecting upon the historical developments of physical theories, it seems worthwhile 
to carry out the symmetry approach to the HTSC problem to its full logical conclusion. The basic idea of unifying 
seemingly different phases by a common symmetry principle may also prove to be useful for other strongly correlated 
systems. 



IX. NOTATIONS AND CONVENTIONS 
A. Index convention 

r" denote Pauli matrices. 

a, (3 = x,y, z denote 50(3) vector spin indices. 

cr, cr' = 1, 2 denote 5*0(3) spinor indices. 

a,b,c— 1, 2, 3, 4, 5 denote 50(5) superspin vector indices. 

IM.v — 1, 2, 3, 4 denote 50(5) spinor indices. 

i,j = 1,5 denote J7(l) vector indices for superconductivity. 

X, x' denote site indices. 



B. Dirac F matrices 

The general method introduced by Rabello et. al ljRabello et alV Il998|) to construct 50(5) symmetric models uses 
the five Dirac F matrices Fq (a = 1, ..,5) which satisfy the Clifford algebra, 

|pap&|^25'^^ (121) 

Rabello et. al introduced the following explicit representation which is naturally adapted for discussing the unification 
of AF and d-wave SC order parameters, 

Here r — (t^, Ty,Tz) are the usual Pauli matrices and *t denotes their transposition. These five F^ matrices form the 
5 dimensional vector irreps of 5*0(5). Their commutators, 

pah ^ _^ ^pa^pbj ^ (123) 

define the 10 dimensional antisymmetric tensor irreps of 50(5). In the above representation, the 10 F"'"s are given 
explicitly by 



F 
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These F matrices satisfy the following commutation relations: 

[pab^pc] ^ 2i{Sacr^ - SbcT'') (124) 

[pab^pcdj ^ 2i{6,j:^'^ + 6Mr'''-Sadr^'-6bcT'''^). (125) 

An important property of the 5*0(5) Lie algebra is the pseudo-reality of its spinor representation. This means that 
there exists a matrix R with the following properties: 

i?2 = = =tR = -R^ (126) 

i?r"j? = -*r", i?r<'''i? = *r"''. (127) 

The relations RV^^R^^ = — (r"**)* indicate that the spinor representation is real, and the antisymmetric nature of 
the matrix R indicates that it is pseudo-real. The R matrix plays a role similar to that of Caf} in 5*0(3). In our 
representation, the R matrix takes the form 

R=[\l). (128) 
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FIG. 1 Phase diagram of the NCCO and the YBCO superconductors. 
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FIG. 2 (a) The spin-flop transition of the XXZ Heisenberg modeL (b) the Mott insulator to superfluid transition of the hard- 
core boson model or the U < Hubbard model, (c) Both can be described as the spin or the pseudospin flop transition in the 
50(3) non-linear a model, induced either by the magnetic field or the chemical potential. 
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FIG. 3 (a) Zero temperature phase diagram of the XXZ Heisenberg model, the hard-core boson model or the negative U 
Hubbard model. Phase I is the Ising or the CDW phase, Phase II is the XY or the superfluid phase and phase III is the fully 
polarized or the normal phase. "Class A" transition is induced by the anisotropy parameter g — J/V , while the "Class B" 
transition is induced by the chemical potential or the magnetic field, (b) Finite temperature phase diagram for the "class A" 
transition in D = 2. Because of the 5*0(3) symmetry at J = V point, the transition temperature vanishes according to the 
Mermin- Wagner theorem. The dashed line denotes the mean field temperature, (c) Finite temperature phase diagram for the 
"class B" transition in D = 3. Ti,c denotes the SO{3) symmetric bicritical point. 
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FIG. 4 Illustration of hopping processes of the magnons and the hole pairs on a ladder. The cross denotes the center of a 
plaquette. An ellipse enclosing two sites denotes a spin singlet, (a) Js describes the magnon hopping, (b) Ji, describes the 
spontaneous creation and annihilation of a magnon pair, (c) Jc describes the hopping of a hole pair, (d) describes the 
spontaneous creation and annihilation of a hole pair and a particle pair. In the full SO{5) model, Js = Ji, and Jc = J'c- In the 
projected 5*0(5) model, the particle pair states are removed and Jc = 0. 




FIG. 5 The chiral 50(5) sphere has an 50(5) symmetric shape but allows only one sense of the rotation in the SC plane 
(ni,n5). Small oscillations around the equator, or the tt triplet resonance, are unaffected by the chiral projection. However, 
small oscillations around the north pole, or the tt doublet mode, are strongly affected: only one of the two such modes is 
retained after the projection. 
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FIG. 6 (a),(b) and (c) express the five bosonic states of the projected SO{5) model in terms of the microscopic states on 
a plaquette. (d), (e) and (f) represent states with well-defined superspin directions, which can be obtained from the linear 
combinations of (a), (b) and (c). These states are analytically defined in Eq. 1551 and Table II. 




FIG. 7 Phase diagram of the 
5c — Ac/4J . Variation of the 
dotted line, x — sin^O and y = 



projected 50(5) model l)5()|l (for the as a function of Ss = As/4 J and 

chemical potential changes Ac and traces out a one- dimensional trajectory as shown on the 
cos^4>. 
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FIG. 8 Spectra of the collective excitations of the projected 50(5) model as a function of density. The region < p < po 
corresponds to the uniform mixed phase of SC and AF. Region p > pc corresponds to the SC phase. The left panel shows the 
spectra of the spin excitations. For p < pc, there axe two gapless spin wave modes and one gapped spin amplitude mode. For 
p > Pc there is a spin triplet tt resonance mode. The right panel shows the spectra of the gapless charge excitations (in the 
absence of long range interactions). For p < pc the charge mode has quadratic dispersion. The dispersion relation changes 
from u; oc fe^ to w oc fc for the p > Pc regime. 




FIG. 9 The energy (a) and the free energy (b) can depend on the density of a uniform AF/SC mixed state with a negative 
curvature when ui2 > ^/urm, (classified as "type 1") or a positive curvature when wi2 < ^miW2 (classified as "type 2"). The 
5(9(5) symmetric limiting case of zero curvature, classified as "type 1.5," is realized when M12 = \/wiW2- (c) The "type 1" 

phase transition from the AF to SC state is a direct first order transition. There are two second order transitions from the AF 
to SC state in the "type 2" case. 50(5) symmetry is realized at the intermediate case of "type 1.5." 
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FIG. 10 The finite temperature piiase diagram in D=3 for tfie class "Bl" transition shown in Fig. 1131 (a) and (b), correspond 
to a direct first order phase transition between AF and SC, as a function of the chemical potential and doping, respectively. 
This type of transition is classified as the "type 1" transition, (c) corresponds to two second order phase transitions with a 
uniform AF/SC mix phase in between. This type of transition is classified as "type 2" transition. The AF and SC transition 
temperatures Tm and Tc merge into either a bi-critical Tbc or a tetra-critical point Ttc- 




FIG. 11 Renormalization group fiow in the (ui/iii2, W2/M12) plane. (In this figure, u = ui, v = U2 and w — ui2.) The 
renormalization group fiow is initially attracted towards the symmetric Heisenberg poi nt labelled by H. The RG traje ctories 
diverge near the Heisenberg model, with a very small exponent. Reproduced from Ref. ijMurakami and Nagaosal l200ff) . 
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FIG. 12 Some possible ground states of the projected 50(5) model (see also Fig. |^. The cross depicts an RVB like spin 
singlet state on a plaquette, the arrow denotes the direction of the superspin, and the open circles depict hole pairs, (a) The 
plaquette RVB state is described by 6{x) = on every plaquette. (b) The in-phase SC stripe with a{x) = 0, n/2, 0, n/2 on each 
stripe, (c) The superspin spiral with a{x) = 0, tt/2, tt, 3tt/2 on each stripe, (d) The hole pair checkerboard state with q(x) — 
everywhere, except on the hole pair plaquettes, where 9 — ■k/2 and a — -k 12. 




FIG. 13 A typical global phase diagram of the extended 50(5) model in the parameter space of chemical potential and the 
ratio of boson hopping energy over interaction energy (see Ref. fphcn et ai, 2003a) for details). This phase diagram shows 
self-similarity among the insulating states at half-filling and other rational filling fractions. There are two types of superfluid- 
insulator transition. The quantum phase transition of "class A" can be approached by varying the hopping energy, for example, 
by applying a pressure and magnetic field at constant doping. The quantum phase transition of "class B" can be realized by 
changing the chemical potential or doping. This theoretical phase diagram can be compared with the global phase diagram of 
the HTSC cuprates. Different families of cuprates correspond to different traces of "class B." For example, we believe YBCO is 
Bl-like, BSCO may be close to i32-like and LSCO is i33-like. The vertical dash-dot line denotes a boundary in the overdoped 
region beyond which our pure bosonic model becomes less accurate. All the phase boundaries in this figure can be classified 
into direct first order (type 1), two second order (type 2), or a marginal case with enhanced symmetry (type 1.5). Type 2 
transitions between CDW lobes and the superconducting state lead to intermediate supersolid phases. 
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FIG. 14 (a) Hole concentration 6 = p/2 = ^{t]^tf^) (filled squares) and magnon density ^{tata) (circles) as a function of 
the chemical potential fi at T/J — 0.03. The small inlay shows a detailed view to the /i region in which the hole-pair density 
jumps to a finite value, (b) Hole densities of the coexisting phases on the first order transition line from (almost) zero to finite 
hole density at /j. — /ic as a function of temperature. 
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FIG. 15 (a) Phase diagram of the projected SO{5) model (see Eq. 1)50^ with Jc = Js) in D—2: The squares between S and 
the tricritical point P trace the first-order line of phase separation. The solid line from P to the right edge of the plot traces 
the Kosterlitz-Thouless transition between the SC and the normal state. The dashed line separating Nt (^triplet dominated 
region) and Nh (=hole pair dominated region) describes the line of equal AF and SC correlation lengths. The small inlay shows 
the same phase diagram on a larger fi scale, covering the whole KT phase. The tricritical point P appears as a result of the 
Mermin- Wagner theorem, which does not allow spin ordering in D=2 at finite temperature, (b) Energy of a single magnon 
excitations in the projected SO (5) model as a function of the chemical potential. This corresponds to the resonance energy of 
the (7r,7r) peak of the spin correlations in the fermionic model (magnons are defined to carry the momentum of the AF order). 
The numbe rs in parentheses indicate the peak weights, i.e. the area under the peak. (20x20 lattice at temperature T/J = 0.1). 
(from Ref. llJostingmeier et all\200^ ) 
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FIG. 16 The phase diagram of the extended SO{5) model obtained by the QMC simulation. The parameters used in simulation 
are As = 4.8, Vc ~ 4.1010, Vc ~ 3.6329 and = V-n ~ 0. The lines are guides to the eye only. The overall topology of the 
phase diagram agrees well with the global phase diagram presented in Fig. 1131 
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FIG. 17 Scaling of the superfluid density near the SO{5) bi-critical point obtained by the classical Monte Carlo simulations. 
The critical behavior of the superfluid densi ty for various g fit into a single scaling curve, from which 5*0(5) scaling exponents 
were obtained. Reproduced from Ref. fiHull200ll) . 
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(a) (b) 

FIG. 18 (a) The r(/i) phase diagram of the three-dimensional projected 50(5) model with Js = Jc/2 and As = Ac — J. N;, 
is the hole-pair dominated part, Nt the triplet dominated part of the high-temperature phase without long-range order. The 
separation line between Nh and Nt is the line of equal spatial correlation decay of hole-pairs and bosons, (b) The T{5) phase 
diagram of the 3D projected 5*0(5) model as a function of hole doping 5 = nh/2. The first order transition line from S to P in 
the T{fi) diagram becomes a "forbidden region" due to phase separation. These two phase diagrams are consistent with those 
presented in Fig. llUb -b based on general arguments. 



1 

0. 99 

T 




-0.13 -0.12 -0.11 -0.1 -0.095 -0.09 -0.085 -0.08 -0.075 

/J'/J ^'/J 



FIG. 19 Scaling of Tn and Tc near the SO{5) bi-critical point. Both Tjv and Tc merge into the bi-critical point tangentially, 
with the crossover exponent of (jj — 1.43 ± 0.05. 
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FIG. 20 (a) Under the condition specified by Eq. 1731 . the 16 states on a rung are classified into 6 groups, each transforming 
irreducibly under the 50(5) group. is an 50(5) singlet state; na\fl) desribes five states that transform as SO{5) vectors; 
^/alfJ) are four states that form an SO{5) spinor; four states also correspond to a spinor; Ra/s'if c,'^ f3\iTj and Rd^'^a'^^plO,) 

are two 5*0(5) singlet states. The figure also gives energies of all multiplets for the 50(5) symmetric ladder model described 
by equations 17211 and (I73II . (b) Strong coupling phase diagram of the 50(5) symmetric ladder model in the {U, V) space. The 
Eo, El and E3 phases are regions in parameter space where the respective states have the lowest energy. 
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FIG. 21 Phase diagram of the bi-layer 50(5) model plotted as J\\/J versus fi. The entire phase transition line from the Mott 
phase into any of the ordered phases is a second order quantum phase transition. The Mott insulating state has 5 massive 
collective modes. The 50(5) symmetric AF/SC uniform mixed state at half-filling has 4 gapless collective modes. The SC 
state has a spin triplet tt resonance mode and one massless charge Goldstone mode. 
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FIG. 22 Evolution of the quasi-particle states when doping is reduced, (a) pure d wave SC gap with nodal quasi-particles. (b) 
the pure d wave SC gap is rotated into an AF gap of the form | cosp^ — cos py\. (c) A large uniform component of the AF/Mott 
insulating gap is developed on top of the | cospx — cospyj gap when doping is reduced close to zero. 
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FIG. 23 Eigenstates the spin 3/2 problem on a single site. The longer (shorter) arrows denote — Sz — = 3/2 (1/2) and the up 
(down) direction denote the +' (-) sign. The El, 4, 6 (singlet), E2,5 (quartet), and E3 (quintet) sets can also be classified as 
50(5) singlet, spinor, and vector representations. 
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FIG. 24 Doping dependence of the ground-state energy (two upper curves) and staggered magnetization (lower curve) for the 
t~ J model with J/t = 0.3. The state with umform AF and d-wave SC order has lower energy compared with the pure d-wave 
SC state for < S < 10%, furthermore, the energy of the uniform AF/SC mixed st ate depends linearly on 5, fitting into the 
SO{5) symmetric "type 1.5" transition classified in Fig. (|5J. Reproduced from Ref. iHimeda and Oeatal . Il999l) . 
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FIG. 25 Exact diagonalization results for the dynamical cor relation function of the vr-operator on a 10-site Hubbard system 
with U = 8t reproduced from reference iMeix ner et aLlligOTT) . A single 5- function-like peak with pronounced weight is visible 
near lu — for the 7r-operator, proving the eigenoperator relation 1881 in the low-energy regime. This "precession frequency" 
lJtt decreases with decreasing doping. An alternatively constructed "s-wave 7r-operator" , with g{p) in Eq. 1281 given by 
g{p) = cospx -\- cospy, shown in the bottom graph exhibits only incoherent behavior and hardly any weight (note the difference 
in the j/-scale). Here {n} denotes average electron density, with (n) = 1 being at half-filling. 
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FIG. 26 After expanding out the coherent state Ij55|l . we obtain the magnon and hole pair states at level v, which is the 
total number of magnons and hole pairs. These states are classified by their {Sz,Q) quantum numbers in (a). The energy is 
independent of the Sz quantum number because of the SO{3) spin rotation symmetry. The energy can depend on Q with three 
generic possibilities, as depicted in (b). (Compare with Fig. |UJ. If the energy depends linearly on Q, there is no free energy 
cost to rotate magnons and hole pairs into each other, and the potential energy is SO{5) symmetric. This multiplet structure 
was tested in the t — J model and shown in Fig. 1271 
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FIG. 27 The low energy states within each total spin and charge sector {Sz,Q) of the 18-site cluster t — J model with 
J/t — 0.5. The states are grouped into different multiplets and are labelled by the spin, charge, point group symmetry, and 
total momentum. Ai denotes the totally symmetric, Bi the d^2_y2-like representation of the C4v symmetry group. The 
quantum numbers of these states match that of the magnon and hole pair states shown in Fig. 1261 Furthermore, the energy 
depends approximately linearly on Q, demonstrating the 50(5) symmetry of the interaction potential among the magnons and 
hole pairs. 
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FIG. 28 Low energy spectrum of the Hubbard model on a plaquette. Eigenstates by total spin S and plaquette momentum 
QxjQy ~ 0,7r. Truncated high energy states are shaded. The vacuum is defined as | Jl), and quantized operators connect the 
vacuum to the lowest eigenstates as shown. (In this figure, denotes the magnon creation operator tj,, and denotes the hole 
pair creation operator t]^.) 




FIG. 29 Illustration of the basic idea of the CORE method. To implement the CORE method, first decompose the original 
lattice in plaquettes, and then truncate the spectrum of a given plaquette to five lowest states, i.e, singlet, hole-pair and three 
magnon states. An effective Hamiltonian for these bosons can then be calculated using the CORE method. Left: local bosons 
in the original lattice. Gray rectangle denotes the singlet RVB vacua, circles denote holes and the set of two parallel vertical 
arrows denote the magnon. Right: local bosons on the lattice of plaquette. Leaf-like pattern denotes a local d-wave hole-pair 
on a plaquette. Canted arrow denotes local magnon on a plaquette. The singlet RVB vacuum is denoted by an empty site. 
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FIG. 30 This figure illustrates the construction of the "superblock" and its Hamiltonian Hs out of two neighboring blocks, 
with intrablock Hamiltonian Ho and interblock coupling V (in the block basis: (Ho)^ „i = (ctn 1 H | a°i) = e^(5„^„' and 
(Vn) / / ~ (ce°a° I V 1 a°,a°i)). 
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FIG. 31 Boson hopping energies versus Hubbard U . The intersection region near [/ = 8 is close to the projected 5'0(5) 
symmetry point. All energies are in units of t. 



( Neel ) 

^^{2,3,4} 




FIG. 32 The order parameter space of the 50(5) theory, tt operator performs a rotation between the AF and the d-wave SC 
states. This small fluctuation is the new Goldstone mode of the 50(5) theory. 



74 





(0,Ji) 


q 








(it.O) 

























7t-5 7C 



FIG. 33 The two particle continuum and the tt excitation for the tight-biding model. Note, that the continuum of two particle 
states collapses to a point when the center of mass momentum is 11 = (tTjTt). The 7r-mode emerges as an anti-bound state 
above the continuum. 
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FIG. 34 Dyson's equation for the 7r-resonance. Function d{p) is defined in equation 12611 . 
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FIG. 36 Spin susceptibility in the SC state for the model II1U6II iPemler et aZ.l IT998bt) . The wavevector is along the to (tt, tt) 
direction. Susceptibility was computed using the self-consistent linear response formalism in Fig. 1351 The peak at (tt, tt) comes 
from the vr-resonance. 
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FIG. 37 Feynmann diagram for the vr-resonance below Tc contrasted with the diagram above Tc. The cross denotes the 
anomalous scattering in the SC state which converts a particle into a hole, and vice versa. 
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FIG. 38 Suppression of the resonance intensity by the magnetic field. Reproduced from Ref. llDai et aLll200i 
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FIG. 39 Doping dependenc e of the resonance energy and intensity measured in neutron scattering experiments. Reproduced 
from Ref. jFong 0^120001) . 
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FIG. 40 Temperature dependence of the resonance intensity compared to the specific heat. Reproduced from Ref. llDai et al\. 
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FIG. 41 SC pairing between electrons mediated by exciting a virtual magnon-vr-mode pair. 




FIG. 42 Diagrammatic representation of the SC pairing mediated by exciting a virtual magnon-7r-mode pair. Solid lines describe 
electron propogators and dashed lines describe interactions. The upper particle-hole ladder corresponds to the magnon and 
the lower particle-particle ladder corresponds to the 7r-mode. 
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FIG. 43 SC vortex with AF core. Far from the center of the vortex core, the superspin vector lies in the SC plane and winds 
around the vortex core by 2ti. The superspin vector lifts up to the AF direction as it approaches the center of the vortex core. 
The arrows represent the direction of the superspin and the color scale represents the magnitude of the AF order parameter. 




FIG. 4 4 The plot of field dependence of the AF moment for different parameters of the LG theory (defined in Ref. llChen et alV 
l2003lJl 'l. The parameters are pi = p2 = , ri = —1, r2 = —0.85, ui = «2 = 1 and x ~ 42.4. Here the parameters are chosen 
such that the maximum SC order is 1 and the SC coherence length at zero field equals the lattice constant a of the lattice 
model, (a) Field dependence for d ifferent values of ui2. The curvature strongly depends on ui2. (b) Fit to the neutron 
scattering results iKang et 0/1120031) of the Ndi.s^Ceo.isCuOi crystal with ui2 = 0.95. Bc2 is about 6.2T in this sample. The 
experimental data is obtained by subtracting the magnetic field response along the c axis by the magnetic field response in the 
ab plane, so that the response from the Nd moment can be removed. 
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FIG. 45 Illustration of the d wave pair density wave state at s = 1/8. In this state, the d wave hole pairs occupy every four 
non-overlapping plaquettes on the original lattice. The charge unit cell is 4o x 4a. The 50(5) model is defined on the center of 
the non-overlapping plaquettes. Such a state could be realized around the vortex core, whose center is depicted by the cross. 
In the actual realization of this state, the hole pair can be much more extended, and the AF ordering could be much reduced 
from the classical value. 




FIG. 46 The combined phase diagram as a function of lattice density D in CeCu2{Sii-xGex)2 {D < Dc ) and in CeCu2Si2 
{Da < D) under pressure P. Note t hat D oc where V is the unit-cell volume, and D = Dsi[l{VGe — Vsijx/Voe] in the 

former case. Reproduced from Ref. iKitaoka et aZ.l 1200211 . 
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FIG. 47 Doping dependence of the SC transition temperature and magnetic moment for La2-xSrxCu04. Reproduced from 
Ref. (jWakimoto et al... .20011 ') 
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(a) Phase diagram the two dimensional organic superconductor BEDT salt. Reproduced from Ref. llMcKenziel 
(b) Log-log plot of T^^ vs (r- Tc)/r for (A) k-( BEDT-TTF)2Cu[N(CN) 2]C1 (solid squares), and (B) deuterated 



K-(BEDT-TTF)2Cu[N(CN)2]Br (open squares). (Data from jKawamoto ef aZ.lli99 




FIG. 49 Setting of the tunnelling experiment for detecting the triplet particle-particle tt mode in the normal state. 
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FIG. 50 Second order tunnelling diagram that gives rise to the resonant coupling of Cooper pairs and tt excitations in the 
junction shown in Fig. 1491 
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FIG. 51 The Superconductor-Antiferromagnet-Superconductor (SC/AF/SC) junction 




FIG. 52 Predicted current-phase characteristics of a SC/AF/SC junction with different d/dco- 



